Molecular dynamics study of the solution of semiflexible telechelic polymer chains with strongly associating endgroups J.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 it possible to reconcile some apparently contradictory results of previous studies. The present constant-pH MD method is shown to require a moderate increase of computational cost when compared to the usual MD method.

The Bacillus subtilis endospore coat protein CotA shows laccase activity. By using comparative modeling techniques, we were able to derive a model for CotA based on the known x-ray structures of zucchini ascorbate oxidase and Cuprinus cereneus laccase. This model of CotA contains all the structural features of a laccase, including the reactive surface-exposed copper center (T1) and two buried copper centers (T2 and T3). Single amino acid substitutions in the CotA T1 copper center (H497A, or M502L) did not prevent assembly of the mutant proteins into the coat and did not alter the pattern of extractable coat polypeptides. However, in contrast to a wild type strain, both mutants produced unpigmented colonies and spores unable to oxidize syringaldazine (SGZ) and 22-azino-bis-(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS). The CotA protein was purified to homogeneity from an overproducing Escherichia coli strain. The purified CotA shows an absorbance and a EPR spectra typical of blue multicopper oxidases. Optimal enzymatic activity was found at

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.

We present in this work the first molecular simulation study of an enzyme, the serine protease cutinase from Fusarium solani pisi, in two ionic liquids (ILs): 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF6]) and 1-butyl-3-methylimidazolium nitrate ([BMIM][NO(3)]). We tested different water contents in these ILs at room temperature (298 K) and high temperature (343 K), and we observe that the enzyme structure is highly dependent on the amount of water present in the IL media. We show that the enzyme is preferentially stabilized in [BMIM][PF6] at 5-10% (w/w) (weight of water over protein) water content at room temperature. [BMIM][PF6] renders a more nativelike enzyme structure at the same water content of 5-10% (w/w) as previously found for hexane, and the system displays a similar bell-shape-like dependence with the water content in the IL media. [BMIM][PF6] is shown to increase significantly the protein thermostability at high temperatures, especially at low hydration. Our analysis indicates that the enzyme is less stabilized in [BMIM][NO(3)] relative to [BMIM][PF6] at both temperatures, most likely due to the strong affinity of the [NO(3)]- anion toward the protein main chain. These findings are in accordance with the experimental knowledge for these two ionic liquids. We also show that these ILs "strip off" most of the water from the enzyme surface in a degree similar to that found for polar organic solvents such as acetonitrile, and that the remaining waters at the enzyme surface are organized in many small clusters.

The primary and three-dimensional structures of a [NiFe] hydrogenase isolated from D. desulfitricans ATCC 27774 were determined, by nucleotide analysis and single-crystal X-ray crystallography. The three-dimensional structural model was refined to R=0.167 and Rfree=0.223 using data to 1.8 A resolution. Two unique structural features are observed: the [4Fe-4S] cluster nearest the [NiFe] centre has been modified [4Fe-3S-3O] by loss of one sulfur atom and inclusion of three oxygen atoms; a three-fold disorder was observed for Cys536 which binds to the nickel atom in the [NiFe] centre. Also, the bridging sulfur atom that caps the active site was found to have partial occupancy, thus corresponding to a partly activated enzyme. These structural features may have biological relevance. In particular, the two less-populated rotamers of Cys536 may be involved in the activation process of the enzyme, as well as in the catalytic cycle. Molecular modelling studies were carried out on the interaction between this [NiFe] hydrogenase and its physiological partner, the tetrahaem cytochrome c3 from the same organism. The lowest energy docking solutions were found to correspond to an interaction between the haem IV region in tetrahaem cytochrome c3 with the distal [4Fe-4S] cluster in [NiFe] hydrogenase. This interaction should correspond to efficient electron transfer and be physiologically relevant, given the proximity of the two redox centres and the fact that electron transfer decay coupling calculations show high coupling values and a short electron transfer pathway. On the other hand, other docking solutions have been found that, despite showing low electron transfer efficiency, may give clues on possible proton transfer mechanisms between the two molecules.

A united-atom model of 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF(6)]) and 1-butyl-3-methylimidazolium nitrate ([BMIM][NO(3)]) is developed in the framework of the GROMOS96 43A1(1) force field. These two ionic liquids are parametrized, and their equilibrium properties in the 298-363 K temperature range are subjected to validation against known experimental properties, namely, density, self-diffusion, shear viscosity, and isothermal compressibility. The ionic radial/spatial distributions, pi interaction, gauche/trans populations of the butyl tail, and enthalpies of vaporization are also reported. The properties obtained from the molecular dynamics simulations agree with experimental data and have the same temperature dependence. The strengths and weakness of our model are discussed.

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 wtand 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. endoplasmic reticulum quality control ͉ folding ͉ membrane traffic ͉ arginine-framed motifs ͉ trafficking signals

Cutinase from Fusarium solani pisi was genetically modified near the active site, by site-directed mutagenesis, to enhance its activity towards polyethylene terephthalate (PET) and polyamide 6,6 (PA 6,6) fibers. The mutations L81A, N84A, L182A, V184A and L189A were done to enlarge the active site in order to better fit a larger polymer chain. Modeling studies have shown enhanced free energy stabilization of model substrate tetrahedral intermediate (TI) bound at the enzyme active site for all mutants, for both model polymers. L81A and L182A showed an activity increase of four- and five-fold, respectively, when compared with the wild type, for PET fibers. L182A showed the one- and two-fold higher ability to biodegrade aliphatic polyamide substrates. Further studies in aliphatic polyesters seem to indicate that cutinase has higher ability to recognize aliphatic substrates.