Conformational transitions of the serotonin 5-HT3 receptor.

The serotonin 5-HT3 receptor is a pentameric ligand-gated ion channel (pLGIC). It belongs to a large family of receptors that function as allosteric signal transducers across the plasma membrane1,2; upon binding of neurotransmitter molecules to extracellular sites, the receptors undergo complex conformational transitions that result in transient opening of a pore permeable to ions. 5-HT3 receptors are therapeutic targets for emesis and nausea, irritable bowel syndrome and depression3. In spite of several reported pLGIC structures4-8, no clear unifying view has emerged on the conformational transitions involved in channel gating. Here we report four cryo-electron microscopy structures of the full-length mouse 5-HT3 receptor in complex with the anti-emetic drug tropisetron, with serotonin, and with serotonin and a positive allosteric modulator, at resolutions ranging from 3.2 Å to 4.5 Å. The tropisetron-bound structure resembles those obtained with an inhibitory nanobody5 or without ligand9. The other structures include an 'open' state and two ligand-bound states. We present computational insights into the dynamics of the structures, their pore hydration and free-energy profiles, and characterize movements at the gate level and cation accessibility in the pore. Together, these data deepen our understanding of the gating mechanism of pLGICs and capture ligand binding in unprecedented detail.

Structural and functional characterization of DdrC, a novel DNA damage-induced nucleoid associated protein involved in DNA compaction.

Deinococcus radiodurans is a spherical bacterium well-known for its outstanding resistance to DNA-damaging agents. Exposure to such agents leads to drastic changes in the transcriptome of D. radiodurans. In particular, four Deinococcus-specific genes, known as DNA Damage Response genes, are strongly up-regulated and have been shown to contribute to the resistance phenotype of D. radiodurans. One of these, DdrC, is expressed shortly after exposure to γ-radiation and is rapidly recruited to the nucleoid. In vitro, DdrC has been shown to compact circular DNA, circularize linear DNA, anneal complementary DNA strands and protect DNA from nucleases. To shed light on the possible functions of DdrC in D. radiodurans, we determined the crystal structure of the domain-swapped DdrC dimer at a resolution of 2.5 Šand further characterized its DNA binding and compaction properties. Notably, we show that DdrC bears two asymmetric DNA binding sites located on either side of the dimer and can modulate the topology and level of compaction of circular DNA. These findings suggest that DdrC may be a DNA damage-induced nucleoid-associated protein that enhances nucleoid compaction to limit the dispersion of the fragmented genome and facilitate DNA repair after exposure to severe DNA damaging conditions.

Structure, substrate binding, and symmetry of the mitochondrial ADP/ATP carrier in its matrix-open state.

The mitochondrial ADP/ATP carrier (AAC) performs the first and last step in oxidative phosphorylation by exchanging ADP and ATP across the mitochondrial inner membrane. Its optimal function has been shown to be dependent on cardiolipins (CLs), unique phospholipids located almost exclusively in the mitochondrial membrane. In addition, AAC exhibits an enthralling threefold pseudosymmetry, a unique feature of members of the SLC25 family. Recently, its conformation poised for binding of ATP was solved by x-ray crystallography referred to as the matrix state. Binding of the substrate leads to conformational changes that export of ATP to the mitochondrial intermembrane space. In this contribution, we investigate the influence of CLs on the structure, substrate-binding properties, and structural symmetry of the matrix state, employing microsecond-scale molecular dynamics simulations. Our findings demonstrate that CLs play a minor stabilizing role on the AAC structure. The interdomain salt bridges and hydrogen bonds forming the cytoplasmic network and tyrosine braces, which ensure the integrity of the global AAC scaffold, highly benefit from the presence of CLs. Under these conditions, the carrier is found to be organized in a more compact structure in its interior, as revealed by analyses of the electrostatic potential, measure of the AAC cavity aperture, and the substrate-binding assays. Introducing a convenient structure-based symmetry metric, we quantified the structural threefold pseudosymmetry of AAC, not only for the crystallographic structure, but also for conformational states of the carrier explored in the molecular dynamics simulations. Our results suggest that CLs moderately contribute to preserve the pseudosymmetric structure of AAC.

Charge Transfer and Chemo-Mechanical Coupling in Respiratory Complex I.

The mitochondrial respiratory chain, formed by five protein complexes, utilizes energy from catabolic processes to synthesize ATP. Complex I, the first and the largest protein complex of the chain, harvests electrons from NADH to reduce quinone, while pumping protons across the mitochondrial membrane. Detailed knowledge of the working principle of such coupled charge-transfer processes remains, however, fragmentary due to bottlenecks in understanding redox-driven conformational transitions and their interplay with the hydrated proton pathways. Complex I from Thermus thermophilus encases 16 subunits with nine iron-sulfur clusters, reduced by electrons from NADH. Here, employing the latest crystal structure of T. thermophilus complex I, we have used microsecond-scale molecular dynamics simulations to study the chemo-mechanical coupling between redox changes of the iron-sulfur clusters and conformational transitions across complex I. First, we identify the redox switches within complex I, which allosterically couple the dynamics of the quinone binding pocket to the site of NADH reduction. Second, our free-energy calculations reveal that the affinity of the quinone, specifically menaquinone, for the binding-site is higher than that of its reduced, menaquinol form-a design essential for menaquinol release. Remarkably, the barriers to diffusive menaquinone dynamics are lesser than that of the more ubiquitous ubiquinone, and the naphthoquinone headgroup of the former furnishes stronger binding interactions with the pocket, favoring menaquinone for charge transport in T. thermophilus. Our computations are consistent with experimentally validated mutations and hierarchize the key residues into three functional classes, identifying new mutation targets. Third, long-range hydrogen-bond networks connecting the quinone-binding site to the transmembrane subunits are found to be responsible for proton pumping. Put together, the simulations reveal the molecular design principles linking redox reactions to quinone turnover to proton translocation in complex I.

Modulation of membrane permeability by carbon dioxide.

Promoting drug delivery across the biological membrane is a common strategy to improve bioavailability. Inspired by the observation that carbonated alcoholic beverages can increase the absorption rate of ethanol, we speculate that carbon dioxide (CO2 ) molecules could also enhance membrane permeability to drugs. In the present work, we have investigated the effect of CO2 on the permeability of a model membrane formed by 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine lipids to three drug-like molecules, namely, ethanol, 2',3'-dideoxyadenosine, and trimethoprim. The free-energy and fractional-diffusivity profiles underlying membrane translocation were obtained from µs-timescale simulations and combined in the framework of the fractional solubility-diffusion model. We find that addition of CO2 in the lipid environment results in an increase of the membrane permeability to the three substrates. Further analysis of the permeation events reveals that CO2 expands and loosens the membrane, which, in turn, facilitates permeation of the drug-like molecules. © 2019 Wiley Periodicals, Inc.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies.

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

Binding properties of the quaternary assembly protein SPAG1.

In cells, many constituents are able to assemble resulting in large macromolecular machineries possessing very specific biological and physiological functions, e.g. ribosome, spliceosome and proteasome. Assembly of such entities is commonly mediated by transient protein factors. SPAG1 is a multidomain protein, known to participate in the assembly of both the inner and outer dynein arms. These arms are required for the function of sensitive and motile cells. Together with RUVBL1, RUVBL2 and PIH1D2, SPAG1 is a key element of R2SP, a protein complex assisting the quaternary assembly of specific protein clients in a tissue-specific manner and associating with heat shock proteins (HSPs) and regulators. In this study, we have investigated the role of TPR domains of SPAG1 in the recruitment of HSP chaperones by combining biochemical assays, ITC, NMR spectroscopy and molecular dynamics (MD) simulations. First, we propose that only two, out of the three TPR domains, are able to recruit the protein chaperones HSP70 and HSP90. We then focused on one of these TPR domains and elucidated its 3D structure using NMR spectroscopy. Relying on an NMR-driven docking approach and MD simulations, we deciphered its binding interface with the C-terminal tails of both HSP70 and HSP90. Finally, we addressed the biological function of SPAG1 and specifically demonstrated that a SPAG1 sub-fragment, containing a putative P-loop motif, cannot efficiently bind and hydrolyze GTP in vitro Our data challenge the interpretation of SPAG1 possessing GTPase activity. We propose instead that SPAG1 regulates nucleotide hydrolysis activity of the HSP and RUVBL1/2 partners.

Changes in Microenvironment Modulate the B- to A-DNA Transition.

B- to A-DNA transition is known to be sensitive to the macroscopic properties of the solution, such as salt and ethanol concentrations. Microenvironmental effects on DNA conformational transition have been broadly studied. Providing an intuitive picture of how DNA responds to environmental changes is, however, still needed. Analyzing the chemical equilibrium of B-to-A DNA transition at critical concentrations, employing explicit-solvent simulations, is envisioned to help understand such microenvironmental effects. In the present study, free-energy calculations characterizing the B- to A-DNA transition and the distribution of cations were carried out in solvents with different ethanol concentrations. With the addition of ethanol, the most stable structure of DNA changes from the B- to A-form, in agreement with previous experimental observation. In 60% ethanol, a chemical equilibrium is found, showing reversible transition between B- and A-DNA. Analysis of the microenvironment around DNA suggests that with the increase of ethanol concentration, the cations exhibit a significant tendency to move toward the backbone, and mobility of water molecules around the major groove and backbone decreases gradually, leading eventually to a B-to-A transition. The present results provide a free-energy view of DNA microenvironment and of the role of cation motion in the conformational transition.

C-glyco"RGD" as αIIbß3 and αvß integrin ligands for imaging applications: Synthesis, in vitro evaluation and molecular modeling.

The design of conjugates displaying simultaneously high selectivity and high affinity for different subtypes of integrins is a current challenge. The arginine-glycine-aspartic acid amino acid sequence (RGD) is one of the most efficient short peptides targeting these receptors. We report herein the development of linear and cyclic fluoro-C-glycoside"RGD" conjugates, taking advantage of the robustness and hydrophilicity of C-glycosides. As attested by in vitro evaluation, the design of these C-glyco"RGD" with a flexible three-carbon triazolyl linker allows distinct profiles towards αIIbß3 and αvß3 integrins. Molecular-dynamics simulations confirm the suitability of cyclic C-glyco-c(RGDfC) to target αvß3 integrin. These C-glyco"RGD" could become promising biological tools in particular for Positron Emission Tomography imaging.

Conformational polymorphism or structural invariance in DNA photoinduced lesions: implications for repair rates.

DNA photolesions constitute a particularly deleterious class of molecular defects responsible for the insurgence of a vast majority of skin malignant tumors. Dimerization of two adjacent thymines or cytosines mostly gives rise to cyclobutane pyrimidine dimers (CPD) and pyrimidine(6-4)pyrimidone 64-PP as the most common defects. We perform all-atom classical simulations, up to 2 µs, of CPD and 64-PP embedded in a 16-bp duplex, which reveal the constrasted behavior of the two lesions. In particular we evidence a very limited structural deformation induced by CPD while 64-PP is characterized by a complex structural polymorphism. Our simulations also allow to unify the contrasting experimental structural results obtained by nuclear magnetic resonance or Förster Resonant Energy Transfer method, showing that both low and high bent structures are indeed accessible. These contrasting behaviors can also explain repair resistance or the different replication obstruction, and hence the genotoxicity of these two photolesions.

Effects of hydration on the protonation state of a lysine analog crossing a phospholipid bilayer - insights from molecular dynamics and free-energy calculations.

The low bioavailability of most therapeutic compounds is often counterbalanced by association with molecular vectors capable of crossing cell membranes. Previous studies demonstrated that for vectors bearing titratable chemical groups, the translocation process might be accompanied by a change in the protonation state. For simple compounds e.g. a lysine analog, free energy calculations, using a single collective variable, namely the insertion depth, suggest that such a transition could only take place if the amino acid diffuses deep enough into the hydrophobic core of the membrane, a situation thermodynamically unfavorable. Here, we determined the 2D potential of mean force associated with the translocation of lysine across a model membrane using as reaction coordinates not only its location in the bilayer but also its hydration. Our results cogently demonstrate that the change in protonation can result from a small fluctuation in the latter, even at low insertion depth.

Conformational changes of DNA induced by a trans-azobenzene derivative via non-covalent interactions.

In biological environments and in aqueous solution, DNA generally adopts the canonical B conformation. Recently, an azobenzene photoswitch containing a polyamine chain with three positive charges was shown to induce a reversible conformational transition between the A and B forms of DNA, the transition being triggered by trans-cis isomerization of the photoswitch upon non-covalent intercalation. It was proposed that, in its trans conformation, azobenzene stabilizes the A form of DNA. The structural details and the mechanism upon which trans-azobenzene induces the B-to-A DNA transition remain, however, unclear. In the present work, two possible intercalating modes of trans-azobenzene, from the minor groove and from the major groove, were investigated with all-atom molecular-dynamics simulations. Intercalation from the major groove was found to be the most probable binding mode due to favorable electrostatic and π-π stacking interactions. The free-energy profile associated with the B-to-A conformational transition reveals that intercalation from the major groove leads to a conformational change of DNA, showing a slight tendency to interconvert from B- to A-DNA, in agreement with the CD spectrum obtained from the experiment. However, the presence of only one interacting azobenzene is not sufficient to lead to a global conformational change to A-DNA. The present results are expected to serve in the design of DNA switches, which can induce reversible DNA conformational changes.

Correlation of bistranded clustered abasic DNA lesion processing with structural and dynamic DNA helix distortion.

Clustered apurinic/apyrimidinic (AP; abasic) DNA lesions produced by ionizing radiation are by far more cytotoxic than isolated AP lesion entities. The structure and dynamics of a series of seven 23-bp oligonucleotides featuring simple bistranded clustered damage sites, comprising of two AP sites, zero, one, three or five bases 3' or 5' apart from each other, were investigated through 400 ns explicit solvent molecular dynamics simulations. They provide representative structures of synthetically engineered multiply damage sites-containing oligonucleotides whose repair was investigated experimentally (Nucl. Acids Res. 2004, 32:5609-5620; Nucl. Acids Res. 2002, 30: 2800-2808). The inspection of extrahelical positioning of the AP sites, bulge and non Watson-Crick hydrogen bonding corroborates the experimental measurements of repair efficiencies by bacterial or human AP endonucleases Nfo and APE1, respectively. This study provides unprecedented knowledge into the structure and dynamics of clustered abasic DNA lesions, notably rationalizing the non-symmetry with respect to 3' to 5' position. In addition, it provides strong mechanistic insights and basis for future studies on the effects of clustered DNA damage on the recognition and processing of these lesions by bacterial or human DNA repair enzymes specialized in the processing of such lesions.

Accurate Estimation of the Standard Binding Free Energy of Netropsin with DNA.

DNA is the target of chemical compounds (drugs, pollutants, photosensitizers, etc.), which bind through non-covalent interactions. Depending on their structure and their chemical properties, DNA binders can associate to the minor or to the major groove of double-stranded DNA. They can also intercalate between two adjacent base pairs, or even replace one or two base pairs within the DNA double helix. The subsequent biological effects are strongly dependent on the architecture of the binding motif. Discriminating between the different binding patterns is of paramount importance to predict and rationalize the effect of a given compound on DNA. The structural characterization of DNA complexes remains, however, cumbersome at the experimental level. In this contribution, we employed all-atom molecular dynamics simulations to determine the standard binding free energy of DNA with netropsin, a well-characterized antiviral and antimicrobial drug, which associates to the minor groove of double-stranded DNA. To overcome the sampling limitations of classical molecular dynamics simulations, which cannot capture the large change in configurational entropy that accompanies binding, we resort to a series of potentials of mean force calculations involving a set of geometrical restraints acting on collective variables.

Mitochondrial ADP/ATP Carrier in Dodecylphosphocholine Binds Cardiolipins with Non-native Affinity.

Biophysical investigation of membrane proteins generally requires their extraction from native sources using detergents, a step that can lead, possibly irreversibly, to protein denaturation. The propensity of dodecylphosphocholine (DPC), a detergent widely utilized in NMR studies of membrane proteins, to distort their structure has been the subject of much controversy. It has been recently proposed that the binding specificity of the yeast mitochondrial ADP/ATP carrier (yAAC3) toward cardiolipins is preserved in DPC, thereby suggesting that DPC is a suitable environment in which to study membrane proteins. In this communication, we used all-atom molecular dynamics simulations to investigate the specific binding of cardiolipins to yAAC3. Our data demonstrate that the interaction interface observed in a native-like environment differs markedly from that inferred from an NMR investigation in DPC, implying that in this detergent, the protein structure is distorted. We further investigated yAAC3 solubilized in DPC and in the milder dodecylmaltoside with thermal-shift assays. The loss of thermal transition observed in DPC confirms that the protein is no longer properly folded in this environment.

Effect of Meso vs Macro Size of Hierarchical Porous Silica on the Adsorption and Activity of Immobilized ß-Galactosidase.

ß-Galactosidase (ß-Gal) is one of the most important enzymes used in milk processing for improving their nutritional quality and digestibility. Herein, ß-Gal has been entrapped into a meso-macroporous material (average pore size 9 and 200 nm, respectively) prepared by a sol-gel method from a silica precursor and a dispersion of solid lipid nanoparticles in a micelle phase. The physisorption of the enzyme depends on the concentration of the feed solution and on the pore size of the support. The enzyme is preferentially adsorbed either in mesopores or in macropores, depending on its initial concentration. Moreover, this selective adsorption, arising from the oligomeric complexation of the enzyme (monomer/dimer/tetramer), has an effect on the catalytic activity of the material. Indeed, the enzyme encapsulated in macropores is more active than the enzyme immobilized in mesopores. Designed materials containing ß-Gal are of particular interest for food applications and potentially extended to bioconversion, bioremediation, or biosensing when coupling the designed support with other enzymes.

Deciphering the photosensitization mechanisms of hypericin towards biological membranes.

Resorting to state-of-the art molecular modeling and simulation techniques we provide full characterization of the photophysical properties of the naturally occurring hypericin chromophore, currently used in photodynamic therapy. In particular, we reveal the different photophysical pathways leading to intersystem-crossing and hence, triplet manifold population that is necessary for the subsequent production of singlet oxygen. In particular we identify an extended region of quasi-degeneracy between the first singlet excited state and three triplet state surfaces. This energetic factor allows the occurrence of intersystem-crossing even in the presence of a relatively small spin-orbit coupling. Furthermore, thanks to extended all-atom molecular dynamics simulations we provide insight into the interaction of hypericin with lipid bilayers. We demonstrate the formation of stable interactions with the membrane and, in particular, the penetration of hypericin into its hydrophobic core. This organization allows a spatial overlap between hypericin and the lipid oxidizable double bond pointing towards the production of singlet oxygen in close spatial proximity to its reactant, hence favoring photosensitization.

Structural and energetic study of cation-π-cation interactions in proteins.

Cation-π interactions of aromatic rings and positively charged groups are among the most important interactions in structural biology. The role and energetic characteristics of these interactions are well established. However, the occurrence of cation-π-cation interactions is an unexpected motif, which raises intriguing questions about its functional role in proteins. We present a statistical analysis of the occurrence, composition and geometrical preferences of cation-π-cation interactions identified in a set of non-redundant protein structures taken from the Protein Data Bank. Our results demonstrate that this structural motif is observed at a small, albeit non-negligible frequency in proteins, and suggest a preference to establish cation-π-cation motifs with Trp, followed by Tyr and Phe. Furthermore, we have found that cation-π-cation interactions tend to be highly conserved, which supports their structural or functional role. Finally, we have performed an energetic analysis of a representative subset of cation-π-cation complexes combining quantum-chemical and continuum solvation calculations. Our results point out that the protein environment can strongly screen the cation-cation repulsion, leading to an attractive interaction in 64% of the complexes analyzed. Together with the high degree of conservation observed, these results suggest a potential stabilizing role in the protein fold, as demonstrated recently for a miniature protein (Craven et al., J. Am. Chem. Soc. 2016, 138, 1543). From a computational point of view, the significant contribution of non-additive three-body terms challenges the suitability of standard additive force fields for describing cation-π-cation motifs in molecular simulations.

AlaScan: A Graphical User Interface for Alanine Scanning Free-Energy Calculations.

Computation of the free-energy changes that underlie molecular recognition and association has gained significant importance due to its considerable potential in drug discovery. The massive increase of computational power in recent years substantiates the application of more accurate theoretical methods for the calculation of binding free energies. The impact of such advances is the application of parent approaches, like computational alanine scanning, to investigate in silico the effect of amino-acid replacement in protein-ligand and protein-protein complexes, or probe the thermostability of individual proteins. Because human effort represents a significant cost that precludes the routine use of this form of free-energy calculations, minimizing manual intervention constitutes a stringent prerequisite for any such systematic computation. With this objective in mind, we propose a new plug-in, referred to as AlaScan, developed within the popular visualization program VMD to automate the major steps in alanine-scanning calculations, employing free-energy perturbation as implemented in the widely used molecular dynamics code NAMD. The AlaScan plug-in can be utilized upstream, to prepare input files for selected alanine mutations. It can also be utilized downstream to perform the analysis of different alanine-scanning calculations and to report the free-energy estimates in a user-friendly graphical user interface, allowing favorable mutations to be identified at a glance. The plug-in also assists the end-user in assessing the reliability of the calculation through rapid visual inspection.