Vincenzo Venditti

Three-dimensional computation of atom depth in complex molecular structures

MOTIVATION For a complex molecular system the delineation of atom-atom contacts, exposed surface and binding sites represents a fundamental step to predict its interaction with solvent, ligands and other molecules. Recently, atom depth has been also considered as an additional structural descriptor to correlate protein structure with folding and functional properties. The distance between an atom and the nearest water molecule or the closest surface dot has been proposed as a measure of the atom depth, but, in both cases, the 3D character of depth is largely lost. In the present study, a new approach is proposed to calculate atom depths in a way that the molecular shape can be taken into account. RESULTS An algorithm has been developed to calculate intersections between the molecular volume and spheres centered on the atoms whose depth has to be quantified. Many proteins with different size and shape have been chosen to compare the results obtained from distance-based and volume-based depth calculations. From the wealth of experimental data available for hen egg white lysozyme, H/D exchange rates and TEMPOL induced paramagnetic perturbations have been analyzed both in terms of depth indexes and of atom distances to the solvent accessible surface. The algorithm here proposed yields better correlations between experimental data and atom depth, particularly for those atoms which are located near to the protein surface. AVAILABILITY Instructions to obtain source code and the executable program are available either from http://sienabiografix.com or http://sadic.sourceforge.net CONTACT niccolai@unisi.it SUPPLEMENTARY INFORMATION http://www.Sienabiogzefix.com/publication.

Guanidinoneomycin B Recognition of an HIV-1 RNA Helix

Aminoglycoside antibiotics are small‐molecule drugs that bind RNA. The affinity and specificity of aminoglycoside binding to RNA can be increased through chemical modification, such as guanidinylation. Here, we report the binding of guanidinoneomycin B (GNB) to an RNA helix from the HIV‐1 frameshift site. The binding of GNB increases the melting temperature (Tm) of the frameshift‐site RNA by at least 10 °C, to a point at which a melting transition is not even observed in 2 M urea. A structure of the complex was obtained by using multidimensional heteronuclear NMR spectroscopic methods. We also used a novel paramagnetic‐probe assay to identify the site of GNB binding to the surface of the RNA. GNB makes major‐groove contacts to two sets of Watson–Crick bases and is in van der Waals contact with a highly structured ACAA tetraloop. Rings I and II of GNB fit into the major groove and form the binding interface with the RNA, whereas rings III and IV are exposed to the solvent and disordered. The binding of GNB causes a broadening of the major groove across the binding site.

Minimum-energy path for a u6 RNA conformational change involving protonation, base-pair rearrangement and base flipping.

The U6 RNA internal stem-loop (U6 ISL) is a highly conserved domain of the spliceosome that is important for pre-mRNA splicing. The U6 ISL contains an internal loop that is in equilibrium between two conformations controlled by the protonation state of an adenine (pK(a)=6.5). Lower pH favors formation of a protonated C-A(+) wobble pair and base flipping of the adjacent uracil. Higher pH favors stacking of the uracil and allows an essential metal ion to bind at this position. Here, we define the minimal-energy path for this conformational transition. To do this, we solved the U6 ISL structure at higher pH (8.0) in order to eliminate interference from the low-pH conformer. This structure reveals disruption of the protonated C-A(+) pair and formation of a new C-U pair, which explains the preference for a stacked uracil at higher pH. Next, we used nudged elastic band molecular dynamics simulations to calculate the minimum-energy path between the two conformations. Our results indicate that the C-U pair is dynamic, which allows formation of the more stable C-A(+) pair upon adenine protonation. After formation of the C-A(+) pair, the unpaired uracil follows a minor-groove base-flipping pathway. Molecular dynamics simulations suggest that the extrahelical uracil is stabilized by contacts with the adjacent helix.

Protein-thiol substitution or protein dethiolation by thiol/disulfide exchange reactions: the albumin model.

Dethiolation experiments of thiolated albumin with thionitrobenzoic acid and thiols (glutathione, cysteine, homocysteine) were carried out to understand the role of albumin in plasma distribution of thiols and disulfide species by thiol/disulfide (SH/SS) exchange reactions. During these experiments we observed that thiolated albumin underwent thiol substitution (Alb‐SS‐X+RSH↔Alb‐SS‐R+XSH) or dethiolation (Alb‐SS‐X+XSH↔Alb‐SH+XSSX), depending on the different pKa values of thiols involved in protein–thiol mixed disulfides (Alb‐SS‐X). It appeared in these reactions that the compound with lower pKa in mixed disulfide was a good leaving group and that the pKa differences dictated the kind of reaction (substitution or dethiolation). Thionitrobenzoic acid, bound to albumin by mixed disulfide (Alb‐TNB), underwent rapid substitution after thiol addition, forming the corresponding Alb‐SS‐X (peaks at 0.25–1 min). In turn, Alb‐SS‐X were dethiolated by the excess nonprotein SH groups because of the lower pKa value in mixed disulfide with respect to that of other thiols. Dethiolation of Alb‐SS‐X was accompanied by formation of XSSX and Alb‐SH up to equilibrium levels at 35 min, which were different for each thiol. Structures by molecular simulation of thiolated albumin, carried out for understanding the role of sulfur exposure in mixed disulfides in dethiolation process, evidenced that the sulfur exposure is important for the rate but not for determining the kind of reaction (substitution or dethiolation). Our data underline the contribution of SH/SS exchanges to determine levels of various thiols as reduced and oxidized species in human plasma.

Structure, dynamics and biophysics of the cytoplasmic protein-protein complexes of the bacterial phosphoenolpyruvate:sugar phosphotransferase system

The bacterial phosphotransferase system (PTS) couples phosphoryl transfer, via a series of bimolecular protein-protein interactions, to sugar transport across the membrane. The multitude of complexes in the PTS provides a paradigm for studying protein interactions, and for understanding how the same binding surface can specifically recognize a diverse array of targets. Fifteen years of work aimed at solving the solution structures of all soluble protein-protein complexes of the PTS has served as a test bed for developing NMR and integrated hybrid approaches to study larger complexes in solution and to probe transient, spectroscopically invisible states, including encounter complexes. We review these approaches, highlighting the problems that can be tackled with these methods, and summarize the current findings on protein interactions.

Conformational Selection and Substrate Binding Regulate the Monomer/Dimer Equilibrium of the C-terminal domain of Escherichia coli Enzyme I

Background: Conformational changes in the EIC domain of enzyme I upon ligand binding are thought to regulate the phosphotransfer system by modulating the monomer/dimer equilibrium. Results: Binding of phosphoenolpyruvate shifts a preexisting conformational equilibrium in EIC. Conclusion: Conformational selection provides a direct structural link between ligand binding and dimer affinity. Significance: Isolated EIC is an optimal system for investigating dynamic processes regulating EI. The bacterial phosphotransferase system (PTS) is a signal transduction pathway that couples phosphoryl transfer to active sugar transport across the cell membrane. The PTS is initiated by the binding of phosphoenolpyruvate (PEP) to the C-terminal domain (EIC) of enzyme I (EI), a highly conserved protein that is common to all sugar branches of the PTS. EIC exists in a dynamic monomer/dimer equilibrium that is modulated by ligand binding and is thought to regulate the overall PTS. Isolation of EIC has proven challenging, and conformational dynamics within the EIC domain during the catalytic cycle are still largely unknown. Here, we present a robust protocol for expression and purification of recombinant EIC from Escherichia coli and show that isolated EIC is capable of hydrolyzing PEP. NMR analysis and residual dipolar coupling measurements indicate that the isolated EIC domain in solution adopts a stable tertiary fold and quaternary structure that is consistent with previously reported crystallographic data. NMR relaxation dispersion measurements indicate that residues around the PEP binding site and in the β3α3 turn (residues 333–366), which is located at the dimer interface, undergo a rapid transition on the sub-millisecond time scale (with an exchange rate constant of ∼1500 s−1) between major open (∼97%) and minor closed (∼3%) conformations. Upon PEP binding, the β3α3 turn is effectively locked in the closed state by the formation of salt bridges between the phosphate group of PEP and the side chains of Lys340 and Arg358, thereby stabilizing the dimer.

Dynamic equilibrium between closed and partially closed states of the bacterial Enzyme I unveiled by solution NMR and X-ray scattering

Significance The bacterial phosphotransferase system couples phosphoryl transfer to sugar transport across the cell membrane. The first protein in the pathway, Enzyme I (EI), undergoes two large rigid body domain reorientations between an autophosphorylation-competent closed state and an open state that allows subsequent phosphoryl transfer to its downstream protein partner. Simultaneous use of solution X-ray scattering and NMR dipolar coupling data to guide simulated annealing refinement reveals the existence of a dynamic equilibrium between closed and partially closed conformations in a complex of a mutant of EI with phosphoenolpyruvate. The partially closed conformation represents an intermediate in the open-to-closed transition. Enzyme I (EI) is the first component in the bacterial phosphotransferase system, a signal transduction pathway in which phosphoryl transfer through a series of bimolecular protein–protein interactions is coupled to sugar transport across the membrane. EI is a multidomain, 128-kDa homodimer that has been shown to exist in two conformational states related to one another by two large (50–90°) rigid body domain reorientations. The open conformation of apo EI allows phosphoryl transfer from His189 located in the N-terminal domain α/β (EINα/β) subdomain to the downstream protein partner bound to the EINα subdomain. The closed conformation, observed in a trapped phosphoryl transfer intermediate, brings the EINα/β subdomain into close proximity to the C-terminal dimerization domain (EIC), thereby permitting in-line phosphoryl transfer from phosphoenolpyruvate (PEP) bound to EIC to His189. Here, we investigate the solution conformation of a complex of an active site mutant of EI (H189A) with PEP. Simulated annealing refinement driven simultaneously by solution small angle X-ray scattering and NMR residual dipolar coupling data demonstrates unambiguously that the EI(H189A)–PEP complex exists in a dynamic equilibrium between two approximately equally populated conformational states, one corresponding to the closed structure and the other to a partially closed species. The latter likely represents an intermediate in the open-to-closed transition.

Tertiary structure prediction of SARS coronavirus helicase

Abstract SARS coronavirus, SCV, has been recently responsible of a sudden and widespread infection which caused almost 800 victims. The limited amount of SCV protein structural information is partially responsible of the lack of specific drugs against the virus. Coronavirus helicases are very conserved and peculiar proteins which have been proposed as suitable targets for antiviral drugs, such as bananins, which have been recently shown to inhibit the SCV helicase in vitro. Here, the quaternary structure of SCV helicase has been predicted, which will provide a solid foundation for the rational design of other antiviral helicase inhibitors.

NMR studies on the surface accessibility of the archaeal protein Sso7d by using TEMPOL and Gd(III)(DTPA-BMA) as paramagnetic probes

Understanding how proteins are approached by surrounding molecules is fundamental to increase our knowledge of life at atomic resolution. Here, the surface accessibility of a multifunctional small protein, the archaeal protein Sso7d from Sulfolobus solfataricus, has been investigated by using TEMPOL and Gd(III)(DTPA-BMA) as paramagnetic probes. The DNA binding domain of Sso7d appears very accessible both to TEMPOL and Gd(III)(DTPA-BMA). Differences in paramagnetic attenuation profiles of (1)H-(15)N HSQC protein backbone amide correlations, observed in the presence of the latter paramagnetic probes, are consistent with the hydrogen bond acceptor capability of the N-oxyl moiety of TEMPOL to surface exposed Sso7d amide groups. By using the gadolinium complex as a paramagnetic probe a better agreement between Sso7d structural features and attenuation profile is achieved. It is interesting to note that the protein P-loop region, in spite of the high surface exposure predicted by the available protein structures, is not approached by TEMPOL and only partially by Gd(III)(DTPA-BMA).

Large interdomain rearrangement triggered by suppression of micro- to millisecond dynamics in bacterial Enzyme I

Enzyme I (EI), the first component of the bacterial phosphotransfer signal transduction system, undergoes one of the largest substrate-induced interdomain rearrangements documented to date. Here, we characterize the perturbations generated by two small molecules, the natural substrate phosphoenolpyruvate (PEP) and the inhibitor α-ketoglutarate (αKG), on the structure and dynamics of EI using NMR, small-angle X-ray scattering (SAXS) and biochemical techniques. The results indicate unambiguously that the open-to-closed conformational switch of EI is triggered by complete suppression of micro- to millisecond dynamics within the C-terminal domain of EI. Indeed, we show that a ligand-induced transition from a dynamic to a more rigid conformational state of the C-terminal domain stabilizes the interface between the N- and C-terminal domains observed in the structure of the closed state, thereby promoting the resulting conformational switch and autophosphorylation of EI. The mechanisms described here may be common to several other multidomain proteins and allosteric systems.