Cláudio M. Soares

Constant-pH molecular dynamics using stochastic titration

A new method is proposed for performing constant-pH molecular dynamics (MD) simulations, that is, MD simulations where pH is one of the external thermodynamic parameters, like the temperature or the pressure. The protonation state of each titrable site in the solute is allowed to change during a molecular mechanics (MM) MD simulation, the new states being obtained from a combination of continuum electrostatics (CE) calculations and Monte Carlo (MC) simulation of protonation equilibrium. The coupling between the MM/MD and CE/MC algorithms is done in a way that ensures a proper Markov chain, sampling from the intended semigrand canonical distribution. This stochastic titration method is applied to succinic acid, aimed at illustrating the method and examining the choice of its adjustable parameters. The complete titration of succinic acid, using constant-pH MD simulations at different pH values, gives a clear picture of the coupling between the trans/gauche isomerization and the protonation process, making i...

Revertant mutants G550E and 4RK rescue cystic fibrosis mutants in the first nucleotide-binding domain of CFTR by different mechanisms

The revertant mutations G550E and 4RK [the simultaneous mutation of four arginine-framed tripeptides (AFTs): R29K, R516K, R555K, and R766K] rescue the cell surface expression and function of F508del-cystic fibrosis (CF) transmembrane conductance regulator (-CFTR), the most common CF mutation. Here, we investigate their mechanism of action by using biochemical and functional assays to examine their effects on F508del and three CF mutations (R560T, A561E, and V562I) located within a conserved region of the first nucleotide-binding domain (NBD1) of CFTR. Like F508del, R560T and A561E disrupt CFTR trafficking. G550E rescued the trafficking defect of A561E but not that of R560T. Of note, the processing and function of V562I were equivalent to that of wild-type (wt)-CFTR, suggesting that V562I is not a disease-causing mutation. Biochemical studies revealed that 4RK generates higher steady-state levels of mature CFTR (band C) for wt- and V562I-CFTR than does G550E. Moreover, functional studies showed that the revertants rescue the gating defect of F508del-CFTR with different efficacies. 4RK modestly increased F508del-CFTR activity by prolonging channel openings, whereas G550E restored F508del-CFTR activity to wt levels by altering the duration of channel openings and closings. Thus, our data suggest that the revertants G550E and 4RK might rescue F508del-CFTR by distinct mechanisms. G550E likely alters the conformation of NBD1, whereas 4RK allows F508del-CFTR to escape endoplasmic reticulum retention/retrieval mediated by AFTs. We propose that AFTs might constitute a checkpoint for endoplasmic reticulum quality control.

Crystal Structure of Cardosin A, a Glycosylated and Arg-Gly-Asp-containing Aspartic Proteinase from the Flowers of Cynara cardunculus L.*

Aspartic proteinases (AP) have been widely studied within the living world, but so far no plant AP have been structurally characterized. The refined cardosin A crystallographic structure includes two molecules, built up by two glycosylated peptide chains (31 and 15 kDa each). The fold of cardosin A is typical within the AP family. The glycosyl content is described by 19 sugar rings attached to Asn-67 and Asn-257. They are localized on the molecular surface away from the conserved active site and show a new glycan of the plant complex type. A hydrogen bond between Gln-126 and Manβ4 renders the monosaccharide oxygen O-2 sterically inaccessible to accept a xylosyl residue, therefore explaining the new type of the identified plant glycan. The Arg-Gly-Asp sequence, which has been shown to be involved in recognition of a putative cardosin A receptor, was found in a loop between two β-strands on the molecular surface opposite the active site cleft. Based on the crystal structure, a possible mechanism whereby cardosin A might be orientated at the cell surface of the style to interact with its putative receptor from pollen is proposed. The biological implications of these findings are also discussed.

Simulation of Electron-Proton Coupling with a Monte Carlo Method: Application to Cytochrome c3 Using Continuum Electrostatics

A new method is presented for simulating the simultaneous binding equilibrium of electrons and protons on protein molecules, which makes it possible to study the full equilibrium thermodynamics of redox and protonation processes, including electron-proton coupling. The simulations using this method reflect directly the pH and electrostatic potential of the environment, thus providing a much closer and realistic connection with experimental parameters than do usual methods. By ignoring the full binding equilibrium, calculations usually overlook the twofold effect that binding fluctuations have on the behavior of redox proteins: first, they affect the energy of the system by creating partially occupied sites; second, they affect its entropy by introducing an additional empty/occupied site disorder (here named occupational entropy). The proposed method is applied to cytochrome c3 of Desulfovibrio vulgaris Hildenborough to study its redox properties and electron-proton coupling (redox-Bohr effect), using a continuum electrostatic method based on the linear Poisson-Boltzmann equation. Unlike previous studies using other methods, the full reduction order of the four hemes at physiological pH is successfully predicted. The sites more strongly involved in the redox-Bohr effect are identified by analysis of their titration curves/surfaces and the shifts of their midpoint redox potentials and pKa values. Site-site couplings are analyzed using statistical correlations, a method much more realistic than the usual analysis based on direct interactions. The site found to be more strongly involved in the redox-Bohr effect is propionate D of heme I, in agreement with previous studies; other likely candidates are His67, the N-terminus, and propionate D of heme IV. Even though the present study is limited to equilibrium conditions, the possible role of binding fluctuations in the concerted transfer of protons and electrons under nonequilibrium conditions is also discussed. The occupational entropy contributions to midpoint redox potentials and pKa values are computed and shown to be significant.

Improved modeling of side-chains in proteins with rotamer-based methods: A flexible rotamer model

Side‐chain modeling has a widespread application in many current methods for protein tertiary structure determination, prediction, and design. Of the existing side‐chain modeling methods, rotamer‐based methods are the fastest and most efficient. Classically, a rotamer is conceived as a single, rigid conformation of an amino acid side‐chain. Here, we present a flexible rotamer model in which a rotamer is a continuous ensemble of conformations that cluster around the classic rigid rotamer. We have developed a thermodynamically based method for calculating effective energies for the flexible rotamer. These energies have a one‐to‐one correspondence with the potential energies of the rigid rotamer. Therefore, the flexible rotamer model is completely general and may be used with any rotamer‐based method in substitution of the rigid rotamer model. We have compared the performance of the flexible and rigid rotamer models with one side‐chain modeling method in particular (the self‐consistent mean field theory method) on a set of 20 high quality crystallographic protein structures. For the flexible rotamer model, we obtained average predictions of 85.8% for χ 1, 76.5% for χ 1+2 and 1.34 Å for root‐mean‐square deviation (RMSD); the corresponding values for core residues were 93.0%, 87.7% and 0.70 Å, respectively. These values represent improvements of 7.3% for χ 1, 8.1% for χ 1+2 and 0.23 Å for RMSD over the predictions obtained with the rigid rotamer model under otherwise identical conditions; the corresponding improvements for core residues were 6.9%, 10.5% and 0.43 Å, respectively. We found that the predictions obtained with the flexible rotamer model were also significantly better than those obtained for the same set of proteins with another state‐of‐the‐art side‐chain placement method in the literature, especially for core residues. The flexible rotamer model represents a considerable improvement over the classic rigid rotamer model. It can, therefore, be used with considerable advantage in all rotamer‐based methods commonly applied to protein tertiary structure determination, prediction, and design and also in predictions of free energies in mutational studies. Proteins 1999;37:530–543. ©1999 Wiley‐Liss, Inc.

Modeling hydration mechanisms of enzymes in nonpolar and polar organic solvents.

A comprehensive study of the hydration mechanism of an enzyme in nonaqueous media was done using molecular dynamics simulations in five organic solvents with different polarities, namely, hexane, 3‐pentanone, diisopropyl ether, ethanol, and acetonitrile. In these solvents, the serine protease cutinase from Fusarium solani pisi was increasingly hydrated with 12 different hydration levels ranging from 5% to 100% (w/w) (weight of water/weight of protein). The ability of organic solvents to ‘strip off’ water from the enzyme surface was clearly dependent on the nature of the organic solvent. The rmsd of the enzyme from the crystal structure was shown to be lower at specific hydration levels, depending on the organic solvent used. It was also shown that organic solvents determine the structure and dynamics of water at the enzyme surface. Nonpolar solvents enhance the formation of large clusters of water that are tightly bound to the enzyme, whereas water in polar organic solvents is fragmented in small clusters loosely bound to the enzyme surface. Ions seem to play an important role in the stabilization of exposed charged residues, mainly at low hydration levels. A common feature is found for the preferential localization of water molecules at particular regions of the enzyme surface in all organic solvents: water seems to be localized at equivalent regions of the enzyme surface independently of the organic solvent employed.

Hydrogenases in Desulfovibrio vulgaris Hildenborough: structural and physiologic characterisation of the membrane-bound [NiFeSe] hydrogenase

The genome of Desulfovibrio vulgaris Hildenborough (DvH) encodes for six hydrogenases (Hases), making it an interesting organism to study the role of these proteins in sulphate respiration. In this work we address the role of the [NiFeSe] Hase, found to be the major Hase associated with the cytoplasmic membrane. The purified enzyme displays interesting catalytic properties, such as a very high H2 production activity, which is dependent on the presence of phospholipids or detergent, and resistance to oxygen inactivation since it is isolated aerobically in a Ni(II) oxidation state. Evidence was obtained that the [NiFeSe] Hase is post-translationally modified to include a hydrophobic group bound to the N-terminal, which is responsible for its membrane association. Cleavage of this group originates a soluble, less active form of the enzyme. Sequence analysis shows that [NiFeSe] Hases from Desulfovibrionacae form a separate family from the [NiFe] enzymes of these organisms, and are more closely related to [NiFe] Hases from more distant bacterial species that have a medial [4Fe4S]2+/1+ cluster, but not a selenocysteine. The interaction of the [NiFeSe] Hase with periplasmic cytochromes was investigated and is similar to the [NiFe]1 Hase, with the Type I cytochrome c3 as the preferred electron acceptor. A model of the DvH [NiFeSe] Hase was generated based on the structure of the Desulfomicrobium baculatum enzyme. The structures of the two [NiFeSe] Hases are compared with the structures of [NiFe] Hases, to evaluate the consensual structural differences between the two families. Several conserved residues close to the redox centres were identified, which may be relevant to the higher activity displayed by [NiFeSe] Hases.

Functional role of N-glycosylation from ADAM10 in processing, localization and activity of the enzyme.

A disintegrin and metalloprotease 10 (ADAM10) is a type I transmembrane glycoprotein with four potential N-glycosylation sites (N267, N278, N439 and N551), that cleaves several plasma membrane proteins. In this work, ADAM10 was found to contain high-mannose and complex-type glycans. Individual N-glycosylation site mutants S269A, T280A, S441A, T553A were constructed, and results indicated that all sites were occupied. T280A was found to accumulate in the endoplasmic reticulum as the non-processed precursor of the enzyme. Furthermore, it exhibited only residual levels of metalloprotease activity in vivo towards the L1 cell adhesion molecule, as well as in vitro, using a ProTNF-alpha peptide as substrate. S441A showed increased ADAM10 susceptibility to proteolysis. Mutation of N267, N439 and N551 did not completely abolish enzyme activity, however, reduced levels were found. ADAM10 is sorted into secretory vesicles, the exosomes. Here, a fraction of ADAM10 from exosomes was found to contain more processed N-linked glycans than the cellular enzyme. In conclusion, N-glycosylation is crucial for ADAM10 processing and resistance to proteolysis, and results suggest that it is required for full-enzyme activity.

Mind the gap: cytochrome interactions reveal electron pathways across the periplasm of Shewanella oneidensis MR-1

Extracellular electron transfer is the key metabolic trait that enables some bacteria to play a significant role in the biogeochemical cycling of metals and in bioelectrochemical devices such as microbial fuel cells. In Shewanella oneidensis MR-1, electrons generated in the cytoplasm by catabolic processes must cross the periplasmic space to reach terminal oxidoreductases found at the cell surface. Lack of knowledge on how these electrons flow across the periplasmic space is one of the unresolved issues related with extracellular electron transfer. Using NMR to probe protein-protein interactions, kinetic measurements of electron transfer and electrostatic calculations, we were able to identify protein partners and their docking sites, and determine the dissociation constants. The results showed that both STC (small tetrahaem cytochrome c) and FccA (flavocytochrome c) interact with their redox partners, CymA and MtrA, through a single haem, avoiding the establishment of stable redox complexes capable of spanning the periplasmic space. Furthermore, we verified that the most abundant periplasmic cytochromes STC, FccA and ScyA (monohaem cytochrome c5) do not interact with each other and this is likely to be the consequence of negative surface charges in these proteins. This reveals the co-existence of two non-mixing redox pathways that lead to extracellular electron transfer in S. oneidensis MR-1 established through transient protein interactions.

Conformational changes induced by ATP‐hydrolysis in an ABC transporter: A molecular dynamics study of the Sav1866 exporter

ATP‐Binding Cassette (ABC) transporters are ubiquitous membrane proteins that use energy from ATP binding or/and hydrolysis to actively transport allocrites across membranes. In this study, we identify ATP‐hydrolysis induced conformational changes in a complete ABC exporter (Sav1866) from Staphylococcus aureaus, using molecular dynamics (MD) simulations. By performing MD simulations on the ATP and ADP+IP bound states, we identify the conformational consequences of hydrolysis, showing that the major rearrangements are not restricted to the NBDs, but extend to the transmembrane domains (TMDs) external regions. For the first time, to our knowledge, we see, within the context of a complete transporter, NBD dimer opening in the ADP+IP state in contrast with all ATP‐bound states. This opening results from the dissociation of the ABC signature motif from the nucleotide. In addition, in both states, we observe the opening of a gate entrance in the intracellular loop region leading to the exposure of the TMDs internal cavity to the cytoplasm. To see if this opening was large enough to allow allocrite transport, the adiabatic energy profile for doxorubicin passage was determined. For both states, this profile, although an approximation, is overall downhill from the cytoplasmatic to the extracellular side, and the local energy barriers along the TMDs are relatively small, evidencing the exporter nature of Sav1866. The major difference between states is an energy barrier located in the cytoplasmic gate region, which becomes reduced upon hydrolysis, suggesting that allocrite passage is facilitated, and evidencing a possible molecular mechanism for the active transport in these proteins. Proteins 2011;