Jerome Baudry

Biasing reaction pathways with mechanical force

During the course of chemical reactions, reactant molecules need to surmount an energy barrier to allow their transformation into products. The energy needed for this process is usually provided by heat, light, pressure or electrical potential, which act either by changing the distribution of the reactants on their ground-state potential energy surface or by moving them onto an excited-state potential energy surface and thereby facilitate movement over the energy barrier. A fundamentally different way of initiating or accelerating a reaction is the use of force to deform reacting molecules along a specific direction of the reaction coordinate. Mechanical force has indeed been shown to activate covalent bonds in polymers, but the usual result is chain scission. Here we show that mechanically sensitive chemical groups make it possible to harness the mechanical forces generated when exposing polymer solutions to ultrasound, and that this allows us to accelerate rearrangement reactions and bias reaction pathways to yield products not obtainable from purely thermal or light-induced reactions. We find that when placed within long polymer strands, the trans and cis isomers of a 1,2-disubstituted benzocyclobutene undergo an ultrasound-induced electrocyclic ring opening in a formally conrotatory and formally disrotatory process, respectively, that yield identical products. This contrasts with reaction initiation by light or heat alone, in which case the isomers follow mutually exclusive pathways to different products. Mechanical forces associated with ultrasound can thus clearly alter the shape of potential energy surfaces so that otherwise forbidden or slow processes proceed under mild conditions, with the directionally specific nature of mechanical forces providing a reaction control that is fundamentally different from that achieved by adjusting chemical or physical parameters. Because rearrangement in our system occurs before chain scission, the effect we describe might allow the development of materials that are activated by mechanical stress fields.

Steered molecular dynamics investigations of protein function

Molecular recognition and mechanical properties of proteins govern molecular processes in the cell that can cause disease and can be targeted for drug design. Single molecule measurement techniques have greatly advanced knowledge but cannot resolve enough detail to be interpreted in terms of protein structure. We seek to complement the observations through so-called Steered Molecular Dynamics (SMD) simulations that link directly to experiments and provide atomic-level descriptions of the underlying events. Such a research program has been initiated in our group and has involved, for example, studies of elastic properties of immunoglobulin and fibronectin domains as well as the binding of biotin and avidin. In this article we explain the SMD method and suggest how it can be applied to the function of three systems that are the focus of modern molecular biology research: force transduction by the muscle protein titin and extracellular matrix protein fibronectin, recognition of antibody-antigene pairs, and ion selective conductivity of the K+ channel.

Human TLRs 10 and 1 Share Common Mechanisms of Innate Immune Sensing but Not Signaling

TLRs are central receptors of the innate immune system that drive host inflammation and adaptive immune responses in response to invading microbes. Among human TLRs, TLR10 is the only family member without a defined agonist or function. Phylogenetic analysis reveals that TLR10 is most related to TLR1 and TLR6, both of which mediate immune responses to a variety of microbial and fungal components in cooperation with TLR2. The generation and analysis of chimeric receptors containing the extracellular recognition domain of TLR10 and the intracellular signaling domain of TLR1, revealed that TLR10 senses triacylated lipopeptides and a wide variety of other microbial-derived agonists shared by TLR1, but not TLR6. TLR10 requires TLR2 for innate immune recognition, and these receptors colocalize in the phagosome and physically interact in an agonist-dependent fashion. Computational modeling and mutational analysis of TLR10 showed preservation of the essential TLR2 dimer interface and lipopeptide-binding channel found in TLR1. Coimmunoprecipitation experiments indicate that, similar to TLR2/1, TLR2/10 complexes recruit the proximal adaptor MyD88 to the activated receptor complex. However, TLR10, alone or in cooperation with TLR2, fails to activate typical TLR-induced signaling, including NF-κB–, IL-8–, or IFN-β–driven reporters. We conclude that human TLR10 cooperates with TLR2 in the sensing of microbes and fungi but possesses a signaling function distinct from that of other TLR2 subfamily members.

Classical Force Field Parameters for the Heme Prosthetic Group of Cytochrome c

Accurate force fields are essential for describing biological systems in a molecular dynamics simulation. To analyze the docking of the small redox protein cytochrome c (cyt c) requires simulation parameters for the heme in both the reduced and oxidized states. This work presents parameters for the partial charges and geometries for the heme in both redox states with ligands appropriate to cyt c. The parameters are based on both protein X‐ray structures and ab initio density functional theory (DFT) geometry optimizations at the B3LYP/6‐31G* level. The simulations with the new parameter set reproduce the geometries of the X‐ray structures and the interaction energies between water and heme prosthetic group obtained from B3LYP/6‐31G* calculations. The parameter set developed here will provide new insights into docking processes of heme containing redox proteins.

Molecular dynamics study of the nature and origin of retinal's twisted structure in bacteriorhodopsin.

The planarity of the polyene chain of the retinal chromophore in bacteriorhodopsin is studied using molecular dynamics simulation techniques and applying different force-field parameters and starting crystal structures. The largest deviations from a planar structure are observed for the C(13)==C(14) and C(15)==N(16) double bonds in the retinal Schiff base structure. The other dihedral angles along the polyene chain of the chromophore, although having lower torsional barriers in some cases, do not significantly deviate from the planar structure. The results of the simulations of different mutants of the pigment show that, among the studied amino acids of the binding pocket, the side chain of Trp-86 has the largest impact on the planarity of retinal, and the mutation of this amino acid to alanine leads to chromophore planarity. Deletion of the methyl C(20), removal of a water molecule hydrogen-bonded to H(15), or mutation of other amino acids to alanine did not show any significant influence on the distortion of the chromophore. The results from the present study suggest the importance of the bulky residue of Trp-86 in the isomerization process, in both ground and excited states of the chromophore, and in fine-tuning of the pK(a) of the retinal protonated Schiff base in bacteriorhodopsin. The dark adaptation of the pigment and the last step of the bacteriorhodopsin photocycle imply low barriers against the rotation of the double bonds in the Schiff base region. The twisted double bonds found in the present study are consistent with the proposed mechanism of these ground state isomerization events.

Oligomerization and fibril assembly of the amyloid-β protein

Abstract In this chapter, we attempt to analyze the evolution of the amyloid-β (Aβ) molecular structure from its inception as part of the Aβ precursor protein to its release by the secretases and its extrusion from membrane into an aqueous environment. Biophysical studies suggest that the Aβ peptide sustains a series of transitions from a molecule rich in α-helix to a molecule in which β-strands prevail. It is proposed that initially the extended C-termini of two opposing Aβ dimers form an antiparallel β-sheet and that the subsequent addition of dimers generates a helical Aβ protofilament. Two or more protofilaments create a strand in which the hydrophobic core of the β-sheets is shielded from the aqueous environment by the N-terminal polar domains of the Aβ dimers. Once the nucleation has occurred, the Aβ filament grows in length by the addition of dimers or tetramers.

A survey of aspartate-phenylalanine and glutamate-phenylalanine interactions in the protein data bank: searching for anion-π pairs.

Protein structures are stabilized using noncovalent interactions. In addition to the traditional noncovalent interactions, newer types of interactions are thought to be present in proteins. One such interaction, an anion-π pair, in which the positively charged edge of an aromatic ring interacts with an anion, forming a favorable anion-quadrupole interaction, has been previously proposed [Jackson, M. R., et al. (2007) J. Phys. Chem. B111, 8242-8249]. To study the role of anion-π interactions in stabilizing protein structure, we analyzed pairwise interactions between phenylalanine (Phe) and the anionic amino acids, aspartate (Asp) and glutamate (Glu). Particular emphasis was focused on identification of Phe-Asp or -Glu pairs separated by less than 7 Å in the high-resolution, nonredundant Protein Data Bank. Simplifying Phe to benzene and Asp or Glu to formate molecules facilitated in silico analysis of the pairs. Kitaura-Morokuma energy calculations were performed on roughly 19000 benzene-formate pairs and the resulting energies analyzed as a function of distance and angle. Edgewise interactions typically produced strongly stabilizing interaction energies (-2 to -7.3 kcal/mol), while interactions involving the ring face resulted in weakly stabilizing to repulsive interaction energies. The strongest, most stabilizing interactions were identified as preferentially occurring in buried residues. Anion-π pairs are found throughout protein structures, in helices as well as β strands. Numerous pairs also had nearby cation-π interactions as well as potential π-π stacking. While more than 1000 structures did not contain an anion-π pair, the 3134 remaining structures contained approximately 2.6 anion-π pairs per protein, suggesting it is a reasonably common motif that could contribute to the overall structural stability of a protein.

Key Role of Active-Site Water Molecules in Bacteriorhodopsin Proton-Transfer Reactions

The functional mechanism of the light-driven proton pump protein bacteriorhodopsin depends on the location of water molecules in the active site at various stages of the photocycle and on their roles in the proton-transfer steps. Here, free energy computations indicate that electrostatic interactions favor the presence of a cytoplasmic-side water molecule hydrogen bonding to the retinal Schiff base in the state preceding proton transfer from the retinal Schiff base to Asp85. However, the nonequilibrium nature of the pumping process means that the probability of occupancy of a water molecule in a given site depends both on the free energies of insertion of the water molecule in this and other sites during the preceding photocycle steps and on the kinetic accessibility of these sites on the time scale of the reaction steps. The presence of the cytoplasmic-side water molecule has a dramatic effect on the mechanism of proton transfer: the proton is channeled on the Thr89 side of the retinal, whereas the transfer on the Asp212 side is hindered. Reaction-path simulations and molecular dynamics simulations indicate that the presence of the cytoplasmic-side water molecule permits a low-energy bacteriorhodopsin conformer in which the water molecule bridges the twisted retinal Schiff base and the proton acceptor Asp85. From this low-energy conformer, proton transfer occurs via a concerted mechanism in which the water molecule participates as an intermediate proton carrier.

Arabidopsis thaliana NIP7;1: an anther-specific boric acid transporter of the aquaporin superfamily regulated by an unusual tyrosine in helix 2 of the transport pore.

Plant nodulin-26 intrinsic proteins (NIPs) are members of the aquaporin superfamily that serve as multifunctional transporters of uncharged metabolites. In Arabidopsis thaliana, a specific NIP pore subclass, known as the NIP II proteins, is represented by AtNIP5;1 and AtNIP6;1, which encode channel proteins expressed in roots and leaf nodes, respectively, that participate in the transport of the critical cell wall nutrient boric acid. Modeling of the protein encoded by the AtNIP7;1 gene shows that it is a third member of the NIP II pore subclass in Arabidopsis. However, unlike AtNIP5;1 and AtNIP6;1 proteins, which form constitutive boric acid channels, AtNIP7;1 forms a channel with an extremely low intrinsic boric acid transport activity. Molecular modeling and molecular dynamics simulations of AtNIP7;1 suggest that a conserved tyrosine residue (Tyr81) located in transmembrane helix 2 adjacent to the aromatic arginine (ar/R) pore selectivity region stabilizes a closed pore conformation through interaction with the canonical Arg220 in ar/R region. Substitution of Tyr81 with a Cys residue, characteristic of established NIP boric acid channels, results in opening of the AtNIP7;1 pore that acquires a robust, transport activity for boric acid as well as other NIP II test solutes (glycerol and urea). Substitution of a Phe for Tyr81 also opens the channel, supporting the prediction from MD simulations that hydrogen bond interaction between the Tyr81 phenol group and the ar/R Arg may contribute to the stabilization of a closed pore state. Expression analyses show that AtNIP7;1 is selectively expressed in developing anther tissues of young floral buds of A. thaliana, principally in developing pollen grains of stage 9-11 anthers. Because boric acid is both an essential nutrient as well as a toxic compound at high concentrations, it is proposed that Tyr81 modulates transport and may provide an additional level of regulation for this transporter in male gametophyte development.

Task-parallel message passing interface implementation of Autodock4 for docking of very large databases of compounds using high-performance super-computers

A message passing interface (MPI)‐based implementation (Autodock4.lga.MPI) of the grid‐based docking program Autodock4 has been developed to allow simultaneous and independent docking of multiple compounds on up to thousands of central processing units (CPUs) using the Lamarkian genetic algorithm. The MPI version reads a single binary file containing precalculated grids that represent the protein–ligand interactions, i.e., van der Waals, electrostatic, and desolvation potentials, and needs only two input parameter files for the entire docking run. In comparison, the serial version of Autodock4 reads ASCII grid files and requires one parameter file per compound. The modifications performed result in significantly reduced input/output activity compared with the serial version. Autodock4.lga.MPI scales up to 8192 CPUs with a maximal overhead of 16.3%, of which two thirds is due to input/output operations and one third originates from MPI operations. The optimal docking strategy, which minimizes docking CPU time without lowering the quality of the database enrichments, comprises the docking of ligands preordered from the most to the least flexible and the assignment of the number of energy evaluations as a function of the number of rotatable bounds. In 24 h, on 8192 high‐performance computing CPUs, the present MPI version would allow docking to a rigid protein of about 300K small flexible compounds or 11 million rigid compounds.