Guilherme M. Arantes

Some aspects of acylation of cellulose under homogeneous solution conditions

Commercially available cellulose (Avicell PH101) was successfully acylated under homogeneous solution conditions by the following procedure: 2.0 g of cellulose were stirred with 75 mL of N,N-dimethylacetamide for 1 h at 150°C, 3.5 g of LiCl were added, the temperature was raised to 170°C, ca. 18.5 mL of the solvent were distilled and the suspension was cooled to room temperature and stirred overnight. The temperature of the clear cellulose solution was raised to 110°C, kept at that temperature for 1 h, an acid anhydride was added and the solution stirred at 110°C for additional 4 h. Acetates, propionates, butyrates, and acetate/propionate mixed ester were prepared with excellent control of the degree of substitution, DS, 1 to 3 for acetates, 2 and 3 for propionates and butyrates, and 3 for acetate/propionate. The degree of polymerization of cellulose is negligibly affected under these reaction conditions. The distribution of the acetyl moiety among the three OH groups of the anhydroglucose unit shows a preference for the C6 position.

A Component of the Xanthomonadaceae Type IV Secretion System Combines a VirB7 Motif with a N0 Domain Found in Outer Membrane Transport Proteins

Type IV secretion systems (T4SS) are used by Gram-negative bacteria to translocate protein and DNA substrates across the cell envelope and into target cells. Translocation across the outer membrane is achieved via a ringed tetradecameric outer membrane complex made up of a small VirB7 lipoprotein (normally 30 to 45 residues in the mature form) and the C-terminal domains of the VirB9 and VirB10 subunits. Several species from the genera of Xanthomonas phytopathogens possess an uncharacterized type IV secretion system with some distinguishing features, one of which is an unusually large VirB7 subunit (118 residues in the mature form). Here, we report the NMR and 1.0 Å X-ray structures of the VirB7 subunit from Xanthomonas citri subsp. citri (VirB7XAC2622) and its interaction with VirB9. NMR solution studies show that residues 27–41 of the disordered flexible N-terminal region of VirB7XAC2622 interact specifically with the VirB9 C-terminal domain, resulting in a significant reduction in the conformational freedom of both regions. VirB7XAC2622 has a unique C-terminal domain whose topology is strikingly similar to that of N0 domains found in proteins from different systems involved in transport across the bacterial outer membrane. We show that VirB7XAC2622 oligomerizes through interactions involving conserved residues in the N0 domain and residues 42–49 within the flexible N-terminal region and that these homotropic interactions can persist in the presence of heterotropic interactions with VirB9. Finally, we propose that VirB7XAC2622 oligomerization is compatible with the core complex structure in a manner such that the N0 domains form an extra layer on the perimeter of the tetradecameric ring.

Specific parametrisation of a hybrid potential to simulate reactions in phosphatases

Phosphatases are key biomolecules because they regulate many cellular processes. These enzymes have been studied for many years, but there are still doubts about the catalytic mechanism. Computer simulations can be used to shed light on these questions. Here we develop a new and specific parametrisation, and present extensive tests of a hybrid potential that can be used to reliably simulate reactions catalysed by phosphatases. High level ab initio data for phosphate ester thiolysis and alcoholysis is used in the training set. The parametrised quantum mechanical Hamiltonian reproduces ab initio energies with a root mean-squared deviation of 3 kcal mol(-1) for species along the pathway of various phosphate ester reactions. Preliminary results for simulation with the calibrated hybrid potential of catalysis by the phosphatase VHR indicate the calculated reaction barriers are in very good agreement with experiment.

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties.

Flexibility and inhibitor binding in cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti‐cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo‐ or sulfate‐bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer‐generated conformational ensembles shown here indicate that the last 30–40 residues in the C‐terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C‐terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B‐inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins. Proteins 2010.

Protein thermal denaturation is modulated by central residues in the protein structure network.

Network structural analysis, known as residue interaction networks or graphs (RIN or RIG, respectively) or protein structural networks or graphs (PSN or PSG, respectively), comprises a useful tool for detecting important residues for protein function, stability, folding and allostery. In RIN, the tertiary structure is represented by a network in which residues (nodes) are connected by interactions (edges). Such structural networks have consistently presented a few central residues that are important for shortening the pathways linking any two residues in a protein structure. To experimentally demonstrate that central residues effectively participate in protein properties, mutations were directed to seven central residues of the β‐glucosidase Sfβgly (β‐d‐glucoside glucohydrolase; EC 3.2.1.21). These mutations reduced the thermal stability of the enzyme, as evaluated by changes in transition temperature (Tm) and the denaturation rate at 45 °C. Moreover, mutations directed to the vicinity of a central residue also caused significant decreases in the Tm of Sfβgly and clearly increased the unfolding rate constant at 45 °C. However, mutations at noncentral residues or at surrounding residues did not affect the thermal stability of Sfβgly. Therefore, the data reported in the present study suggest that the perturbation of the central residues reduced the stability of the native structure of Sfβgly. These results are in agreement with previous findings showing that networks are robust, whereas attacks on central nodes cause network failure. Finally, the present study demonstrates that central residues underlie the functional properties of proteins.

Partition, orientation and mobility of ubiquinones in a lipid bilayer

Ubiquinone is the universal mobile charge carrier involved in biological electron transfer processes. Its redox properties and biological function depend on the molecular partition and lateral diffusion over biological membranes. However, ubiquinone localization and dynamics within lipid bilayers are long debated and still uncertain. Here we present molecular dynamics simulations of several ubiquinone homologs with variable isoprenoid tail lengths complexed to phosphatidylcholine bilayers. Initially, a new force-field parametrization for ubiquinone is derived from and compared to high level quantum chemical data. Free energy profiles for ubiquinone insertion in the lipid bilayer are obtained with the new force-field. The profiles allow for the determination of the equilibrium location of ubiquinone in the membrane as well as for the validation of the simulation model by direct comparison with experimental partition coefficients. A detailed analysis of structural properties and interactions shows that the ubiquinone polar head group is localized at the water-bilayer interface at the same depth of the lipid glycerol groups and oriented normal to the membrane plane. Both the localization and orientation of ubiquinone head groups do not change significantly when increasing the number of isoprenoid units. The isoprenoid tail is extended and packed with the lipid acyl chains. For ubiquinones with long tails, the terminal isoprenoid units have high flexibility. Calculated ubiquinone diffusion coefficients are similar to that found for the phosphatidylcholine lipid. These results may have further implications for the mechanisms of ubiquinone transport and binding to respiratory and photosynthetic protein complexes.

Force-induced chemical reactions on the metal centre in a single metalloprotein molecule

Metalloproteins play indispensable roles in biology owing to the versatile chemical reactivity of metal centres. However, studying their reactivity in many metalloproteins is challenging, as protein three-dimensional structure encloses labile metal centres, thus limiting their access to reactants and impeding direct measurements. Here we demonstrate the use of single-molecule atomic force microscopy to induce partial unfolding to expose metal centres in metalloproteins to aqueous solution, thus allowing for studying their chemical reactivity in aqueous solution for the first time. As a proof-of-principle, we demonstrate two chemical reactions for the FeS4 centre in rubredoxin: electrophilic protonation and nucleophilic ligand substitution. Our results show that protonation and ligand substitution result in mechanical destabilization of the FeS4 centre. Quantum chemical calculations corroborated experimental results and revealed detailed reaction mechanisms. We anticipate that this novel approach will provide insights into chemical reactivity of metal centres in metalloproteins under biologically more relevant conditions.

The Catalytic Acid in the Dephosphorylation of the Cdk2-pTpY/CycA Protein Complex by Cdc25B Phosphatase

The development of anticancer therapeutics that target Cdc25 phosphatases is now an active area of research. A complete understanding of the Cdc25 catalytic mechanism would certainly allow a more rational inhibitor design. However, the identity of the catalytic acid used by Cdc25 has been debated and not established unambiguously. Results of molecular dynamics simulations with a calibrated hybrid potential for the first reaction step catalyzed by Cdc25B in complex with its natural substrate, the Cdk2-pTpY/CycA protein complex, are presented here. The calculated reaction free-energy profiles are in very good agreement with experimental measurements and are used to discern between different proposals for the general acid. In addition, the simulations give useful insight on interactions that can be explored for the design of inhibitors specific to Cdc25.

Theoretical modeling of low-energy electronic absorption bands in reduced cobaloximes

The reduced Co(I) states of cobaloximes are powerful nucleophiles that play an important role in the hydrogen-evolving catalytic activity of these species. In this work we analyze the low-energy electronic absorption bands of two cobaloxime systems experimentally and use a variety of density functional theory and molecular orbital ab initio quantum chemical approaches. Overall we find a reasonable qualitative understanding of the electronic excitation spectra of these compounds but show that obtaining quantitative results remains a challenging task.