Metal-binding loop length and not sequence dictates structure

DeepBindBC: a practical deep learning method for identifying native-like protein-ligand complexes in virtual screening

Identifying native-like protein-ligand complexes (PLCs) from an abundance of docking decoys is critical for large-scale virtual drug screening in early-stage drug discovery lead searching efforts. Providing reliable prediction is still a challenge for most current affinity predicting models because of a lack of non-binding data during model training, lost critical physical-chemical features, and difficulties in learning abstract information with limited neural layers. In this work, we proposed a deep learning model, DeepBindBC, for classifying putative ligands as binding or non-binding. Our model incorporates information on non-binding interactions, making it more suitable for real applications. ResNet model architecture and more detailed atom type representation guarantee implicit features can be learned more accurately. Here, we show that DeepBindBC outperforms Autodock Vina, Pafnucy, and DLSCORE for three DUD.E testing sets. Moreover, DeepBindBC identified a novel human pancreatic α-amylase binder validated by a fluorescence spectral experiment (K_a= 1.0×10⁵ M). Furthermore, DeepBindBC can be used as a core component of a hybrid virtual screening pipeline that incorporating many other complementary methods, such as DFCNN, Autodock Vina docking, and pocket molecular dynamics simulation. Additionally, an online web server based on the model is available at http://cbblab.siat.ac.cn/DeepBindBC/index.php for the user's convenience. Our model and the web server provide alternative tools in the early steps of drug discovery by providing accurate identification of native-like PLCs.

A new density-modification procedure extending the application of the recent |ρ|-based phasing algorithm to larger crystal structures

The incorporation of the new peakness-enhancing fast Fourier transform compatible ipp procedure (ipp = inner-pixel preservation) into the recently published SM algorithm based on |ρ| [Rius (2020). Acta Cryst A76, 489-493] improves its phasing efficiency for larger crystal structures with atomic resolution data. Its effectiveness is clearly demonstrated via a collection of test crystal structures (taken from the Protein Data Bank) either starting from random phase values or by using the randomly shifted modulus function (a Patterson-type synthesis) as initial ρ estimate. It has been found that in the presence of medium scatterers (e.g. S or Cl atoms) crystal structures with 1500 × c atoms in the unit cell (c = number of centerings) can be routinely solved. In the presence of strong scatterers like Fe, Cu or Zn atoms this number increases to around 5000 × c atoms. The implementation of this strengthened SM algorithm is simple, since it only includes a few easy-to-adjust parameters.

CuA-based chimeric T1 copper sites allow for independent modulation of reorganization energy and reduction potential

Attaining rational modulation of thermodynamic and kinetic redox parameters of metalloproteins is a key milestone towards the (re)design of proteins with new or improved redox functions. Here we report that implantation of ligand loops from natural T1 proteins into the scaffold of a CuA protein leads to a series of distorted T1-like sites that allow for independent modulation of reduction potentials (E°') and electron transfer reorganization energies (λ). On the one hand E°' values could be fine-tuned over 120 mV without affecting λ. On the other, λ values could be modulated by more than a factor of two while affecting E°' only by a few millivolts. These results are in sharp contrast to previous studies that used T1 cupredoxin folds, thus highlighting the importance of the protein scaffold in determining such parameters.

An NMR Confirmation for Increased Folded State Entropy Following Loop Truncation

Previous studies conducted on flexible loop regions in proteins revealed that the energetic consequences of changing loop length predominantly arise from the entropic cost of ordering a loop during folding. However, in an earlier study of human acylphosphatase (hmAcP) using experimental and computational approaches, we showed that thermodynamic stabilization upon loop truncation can be attributed mainly to the increased entropy of the folded state. Here, using 15N NMR spectroscopy, we studied the effect of loop truncation on hmAcP backbone dynamics on the picosecond-nanosecond timescale with the aim of confirming the effect of folded state entropy on protein stability. NMR-relaxation-derived N-H squared generalized order parameters reveal that loop truncation results in a significant increase in protein conformational flexibility. Comparison of these results with previously acquired all-atom molecular dynamics simulation, analyzed here in terms of squared generalized NMR order parameters, demonstrates general agreement between the two methods. The NMR study not only provides direct evidence for the enhanced conformational entropy of the folded state of hmAcP upon loop truncation but also gives a quantitative measure of the observed effects.

Engineering a bifunctional copper site in the cupredoxin fold by loop-directed mutagenesis

Copper sites in proteins are designed to perform either electron transfer or redox catalysis. Type 1 and CuA sites are electron transfer hubs bound to a rigid protein fold that prevents binding of exogenous ligands and side reactions. Here we report the engineering of two Type 1 sites by loop-directed mutagenesis within a CuA scaffold with unique electronic structures and functional features. A copper-thioether axial bond shorter than the copper-thiolate bond is responsible for the electronic structure features, in contrast to all other natural or chimeric sites where the copper thiolate bond is short. These sites display highly unusual features, such as: (1) a high reduction potential despite a strong interaction with the axial ligand, which we attribute to changes in the hydrogen bond network and (2) the ability to bind exogenous ligands such as imidazole and azide. This strategy widens the possibility of using natural protein scaffolds with functional features not present in nature.

Engineering the metal-binding loop at a type 1 copper center by circular permutation

Circular permutation of the cupredoxin azurin creates a break on the metal binding loop, highlighting the loop's flexibility.

Asn47 and Phe114 modulate the inner sphere reorganization energies of type zero copper proteins

The geometric structures and electron transfer properties of type 1 Cu proteins are reasonably understood at the molecular level (E. I. Solomon and R. G. Hadt, Coord. Chem. Rev., 2011, 255, 774-789, J. J. Warren, K. M. Lancaster, J. H. Richards and H. B. Gray, J. Inorg. Biochem., 2012, 115, 119-126). Much understanding of type 1 copper electron transfer reactivity has come from site directed mutagenesis studies. For example, artificial "type zero" Cu-centres constructed in cupredoxin-azurin have showcased the capacity of outer-sphere hydrogen bonding networks to enhance Cu II/I electron transfer reactivity. In this paper, we have elaborated on earlier kinetics and electronic structural studies of type zero Cu by calculating the inner sphere reorganization energies of type 1, type 2, and type zero Cu proteins using density functional theory (DFT). Although the choice of density functionals for copper systems is not straightforward, we have benchmarked the density functionals against the recently reported ESI-PES data for two synthetic copper models (S. Niu, D.-L. Huang, P. D. Dau, H.-T. Liu, L.-S. Wang and T. J. Ichiye, Chem. Theory Comput., 2014, 10, 1283). For the Cu proteins, our calculations predict that changes in the coordination number upon metal reduction lead to large inner sphere reorganization energies for type 2 Cu sites, whereas retention in the coordination number is observed for type zero Cu sites. These variations in the coordination number are modulated by the outer-sphere coordinating residues Asn47 and Phe114, which are involved in hydrogen bonding with the Asp112 side chain.

Shortening a loop can increase protein native state entropy

Protein loops are essential structural elements that influence not only function but also protein stability and folding rates. It was recently reported that shortening a loop in the AcP protein may increase its native state conformational entropy. This effect on the entropy of the folded state can be much larger than the lower entropic penalty of ordering a shorter loop upon folding, and can therefore result in a more pronounced stabilization than predicted by polymer model for loop closure entropy. In this study, which aims at generalizing the effect of loop length shortening on native state dynamics, we use all-atom molecular dynamics simulations to study how gradual shortening a very long or solvent-exposed loop region in four different proteins can affect their stability. For two proteins, AcP and Ubc7, we show an increase in native state entropy in addition to the known effect of the loop length on the unfolded state entropy. However, for two permutants of SH3 domain, shortening a loop results only with the expected change in the entropy of the unfolded state, which nicely reproduces the observed experimental stabilization. Here, we show that an increase in the native state entropy following loop shortening is not unique to the AcP protein, yet nor is it a general rule that applies to all proteins following the truncation of any loop. This modification of the loop length on the folded state and on the unfolded state may result with a greater effect on protein stability.

Spectroscopic characterization of a green copper site in a single-domain cupredoxin

Cupredoxins are widespread copper-binding proteins, mainly involved in electron transfer pathways. They display a typical rigid greek key motif consisting of an eight stranded β-sandwich. A fascinating feature of cupredoxins is the natural diversity of their copper center geometry. These geometry variations give rise to drastic changes in their color, such as blue, green, red or purple. Based on several spectroscopic and structural analyses, a connection between the geometry of their copper-binding site and their color has been proposed. However, little is known about the relationship between such diversity of copper center geometry in cupredoxins and possible implications for function. This has been difficult to assess, as only a few naturally occurring green and red copper sites have been described so far. We report herein the spectrocopic characterization of a novel kind of single domain cupredoxin of green color, involved in a respiratory pathway of the acidophilic organism Acidithiobacillus ferrooxidans. Biochemical and spectroscopic characterization coupled to bioinformatics analysis reveal the existence of some unusual features for this novel member of the green cupredoxin sub-family. This protein has the highest redox potential reported to date for a green-type cupredoxin. It has a constrained green copper site insensitive to pH or temperature variations. It is a green-type cupredoxin found for the first time in a respiratory pathway. These unique properties might be explained by a region of unknown function never found in other cupredoxins, and by an unusual length of the loop between the second and the fourth copper ligands. These discoveries will impact our knowledge on non-engineered green copper sites, whose involvement in respiratory chains seems more widespread than initially thought.

Molecular dynamics simulations of apocupredoxins: insights into the formation and stabilization of copper sites under entatic control

Cupredoxins perform copper-mediated long-range electron transfer (ET) in biological systems. Their copper-binding sites have evolved to force copper ions into ET-competent systems with decreased reorganization energy, increased reduction potential, and a distinct electronic structure compared with those of non-ET-competent copper complexes. The entatic or rack-induced state hypothesis explains these special properties in terms of the strain that the protein matrix exerts on the metal ions. This idea is supported by X-ray structures of apocupredoxins displaying "closed" arrangements of the copper ligands like those observed in the holoproteins; however, it implies completely buried copper-binding atoms, conflicting with the notion that they must be exposed for copper loading. On the other hand, a recent work based on NMR showed that the copper-binding regions of apocupredoxins are flexible in solution. We have explored five cupredoxins in their "closed" apo forms through molecular dynamics simulations. We observed that prearranged ligand conformations are not stable as the X-ray data suggest, although they do form part of the dynamic landscape of the apoproteins. This translates into variable flexibility of the copper-binding regions within a rigid fold, accompanied by fluctuations of the hydrogen bonds around the copper ligands. Major conformations with solvent-exposed copper-binding atoms could allow initial binding of the copper ions. An eventual subsequent incursion to the closed state would result in binding of the remaining ligands, trapping the closed conformation thanks to the additional binding energy and the fastening of noncovalent interactions that make up the rack.

How the dynamics of the metal-binding loop region controls the acid transition in cupredoxins

Many reduced cupredoxins undergo a pH-dependent structural rearrangement, triggered by protonation of the His ligand belonging to the C-terminal hydrophobic loop, usually termed the acid transition. At variance with several members of the cupredoxin family, the acid transition is not observed for azurin (AZ). We have addressed this issue by performing molecular dynamics simulations of AZ and four mutants, in which the C-terminal loop has been replaced with those of other cupredoxins or with polyalanine loops. All of the loop mutants undergo the acid transition in the pH range of 4.4-5.5. The main differences between AZ and its loop mutants are the average value of the active site solvent accessible surface area and the extent of its fluctuations with time, together with an altered structure of the water layer around the copper center. Using functional mode analysis, we found that these variations arise from changes in nonbonding interactions in the second coordination sphere of the copper center, resulting from the loop mutation. Our results strengthen the view that the dynamics at the site relevant for function and its surroundings are crucial for protein activity and for metal-containing electron transferases.

Enhancement of catalysis and functional expression of a bacterial laccase by single amino acid replacement

Structure-function relationships underlying laccases properties are very limited that makes these enzymes interesting for protein engineering approaches. Therefore in the current study, a thermostable laccase that was isolated from Bacillus sp. HR03 with the ability of bilirubin oxidation besides its laccase and tyrosinase activity is used. The extensive application of this enzyme is limited by its low expression level in Escherichia coli. Based on sequence alignments and structural studies, three single amino acid substitutions, D500G, D500E, D500S and a glycine insertion, are introduced using site-directed mutagenesis to evaluate the role of Asp(500) located in the C-terminal segment close to the T1 copper center. Substitution of aspartic acid with less sterically hindered, conserved residue such as glycine increase kcat (2.3 fold) and total activity (7.3 fold) which is accompanied by a significant increase in the expression level up to 3 fold. Biochemical characterization and structural studies using far-UV CD and fluorescence spectroscopy reveal the importance of C-terminal copper-binding loop in the laccase functional expression and catalytic efficiency. Kinetic characterization of the purified mutants toward 2,2'-azino-bis (3-ethylbenzothiazoline-6-sulfonic acid) (ABTS), syringaldazine (SGZ) and bilirubin, shows that substrate specificity is left unchanged.

The Active Site Loop Modulates the Reorganization Energy of Blue Copper Proteins by Controlling the Dynamic Interplay with Solvent

Understanding the factors governing the rate of electron transfer processes in proteins is crucial not only to a deeper understanding of redox processes in living organisms but also for the design of efficient devices featuring biological molecules. Here, molecular dynamics simulations performed on native azurin and four chimeric cupredoxins allow for the calculation of the reorganization energy and of structure-related quantities that were used to clarify the molecular determinants to the dynamics/function relationship in blue copper proteins. We find that the dynamics of the small, metal-binding loop region controls the outer-sphere reorganization energy not only by determining the exposure of the active site to solvent but also through the modulation of the redox-dependent rearrangement of the whole protein scaffold and of the surrounding water molecules.

The influence of protein folding on the copper affinities of trafficking and target sites

The relative influence of protein unfolding on the Cu(I) affinity of trafficking and target sites for copper has been determined. For the copper metallochaperone Atx1 from Synechocystis PCC 6803 (a cyanobacterium), Saccharomyces cerevisiae and humans unfolding in urea results in a decrease in the Cu(I) affinity from (4-5) × 10(17) M(-1) to (1-3) × 10(16) M(-1) at pH 7. The affinities of the unfolded Atx1s are similar to those for CXXC-containing peptides. Partial unfolding, due to the loop 5 His61Lys mutation in Synechocystis Atx1, gives rise to a more limited decrease in Cu(I) affinity. For the copper target protein plastocyanin from Synechocystis, chemical unfolding results in the Cu(I) affinity decreasing by 5-orders of magnitude. This differential influence of protein unfolding on Cu(I) affinity is due to a more complex copper site structure in the target protein, including numerous interactions of non-coordinating residues with ligating amino acids. This second-coordination sphere is much simpler in the Atx1s with the main interaction provided by the loop 5 residue that tunes the Cu(I) affinity by altering the pK(a) of the C-terminal Cys ligand of the CXXC motif. This interaction and others are absent in the unfolded Atx1s and the two Cys ligands have pK(a) values reminiscent of free thiols (>8) resulting in lowered Cu(I) affinities at pH 7. Residues close to the active site of the thiol-disulfide oxidoreductase thioredoxin appear to lower the Cu(I) affinity of its CXXC motif to 3.1 × 10(15) M(-1) at pH 7, presumably to prevent copper binding in vivo. The structure of a copper site, including the number and relative position of ligands in the primary structure and the complexity of the second-coordination sphere, results in dramatically different effects of unfolding on Cu(I) affinity that has important implications for copper homeostasis.

Functional versatility of a single protein surface in two protein:protein interactions

The ability of the Escherichia coli protein BirA to function as both a metabolic enzyme and a transcription repressor relies on the use of a single surface for two distinct protein:protein interactions. BirA forms a heterodimer with the biotin acceptor protein of acetyl-coenzyme A carboxylase and catalyzes posttranslational biotinylation. Alternatively, it forms a homodimer that binds sequence-specifically to DNA to repress transcription initiation at the biotin biosynthetic operon. Several surface loops on BirA, two of which exhibit sequence conservation in all biotin protein ligases and the remainder of which are highly variable, are located at the two interfaces. The function of these loops in both homodimerization and biotin transfer was investigated by characterizing alanine-substituted variants at 18 positions of one constant and three variable loops. Sedimentation equilibrium measurements reveal that 11 of the substitutions, which are distributed throughout conserved and variable loops, significantly alter homodimerization energetics. By contrast, steady-state and single-turnover kinetic measurements indicate that biotin transfer to biotin carboxyl carrier protein is impacted by seven substitutions, the majority of which are in the constant loop. Furthermore, constant loop residues that function in biotin transfer also support homodimerization. The results reveal clues about the evolution of a single protein surface for use in two distinct functions.

The influence of active site loop mutations on the thermal stability of azurin from Pseudomonas aeruginosa

The copper site and overall structures of azurin (AZ) variants in which the amicyanin (AMI) and plastocyanin (PC) metal binding loops have been introduced, AZAMI and AZPC, respectively, are similar to that of AZ, whereas the loop conformations resemble those in the native proteins. To assess the influence of these loop mutations on stability, the thermal unfolding of AZAMI and AZPC has been investigated by differential scanning calorimetry, absorption and fluorescence spectroscopy. The calorimetric profiles of both variants exhibit a complex shape consisting of two endothermic peaks and an exothermic peak. The temperature of the maximum heat of absorption for the single endothermic peak is 82.7°C for AZ, whereas for AZAMI and AZPC the most intense endothermic peaks are at 74.9 and 68.1°C comparable to values for AMI and PC, respectively. Denaturation investigated using the temperature dependence of the absorbance at ∼600nm and Trp emission, also demonstrates decreased stability for both loop mutants. The thermal transition between the native and the denaturated states is irreversible, scan rate dependent and consistent with the two-state irreversible model. The structure of the active-site loop has a dramatic effect on the kinetic stability and the unfolding pathway of cupredoxins.