Correction to "A 3D Model of Human Lysyl Oxidase, A Cross-Linking Enzyme".

[This corrects the article DOI: 10.1021/acsomega.9b00317.].

Transmembrane Peptides as a New Strategy to Inhibit Neuraminidase-1 Activation.

Sialidases, or neuraminidases, are involved in several human disorders such as neurodegenerative, infectious and cardiovascular diseases, and cancers. Accumulative data have shown that inhibition of neuraminidases, such as NEU1 sialidase, may be a promising pharmacological target, and selective inhibitors of NEU1 are therefore needed to better understand the biological functions of this sialidase. In the present study, we designed interfering peptides (IntPep) that target a transmembrane dimerization interface previously identified in human NEU1 that controls its membrane dimerization and sialidase activity. Two complementary strategies were used to deliver the IntPep into cells, either flanked to a TAT sequence or non-tagged for solubilization in detergent micelles. Combined with molecular dynamics simulations and heteronuclear nuclear magnetic resonance (NMR) studies in membrane-mimicking environments, our results show that these IntPep are able to interact with the dimerization interface of human NEU1, to disrupt membrane NEU1 dimerization and to strongly decrease its sialidase activity at the plasma membrane. In conclusion, we report here new selective inhibitors of human NEU1 of strong interest to elucidate the biological functions of this sialidase.

A Three-Dimensional Model of Human Lysyl Oxidase, a Cross-Linking Enzyme.

Lysyl oxidase (LOX) is a cross-linking enzyme identified 50 years ago, but its 3D structure is still unknown. We have thus built a 3D model of human LOX by homology modeling using the X-ray structure of human lysyl oxidase-like 2 as a template. This model is the first one to recapitulate all known biochemical features of LOX, namely, the coordination of the copper ion and the formation of the lysine tyrosylquinone cofactor and the disulfide bridges. Furthermore, this model is stable during a 1 µs molecular dynamics simulation. The catalytic site is located in a groove surrounded by two loops. The distance between these loops fluctuated during the simulations, which suggests that the groove forms a hinge with a variable opening, which is able to accommodate the various sizes of LOX substrates. This 3D model is a pre-requisite to perform docking experiments with LOX substrates and other partners to identify binding sites and to design new LOX inhibitors specific for therapeutic purpose.

Identification of a crucial amino acid implicated in the hydroxylation/desaturation ratio of CpFAH12 bifunctional hydroxylase.

Claviceps purpurea bifunctional Δ12-hydroxylase/desaturase, CpFAH12, and monofunctional desaturase CpFAD2, share 86% of sequence identity. To identify the underlying determinants of the hydroxylation/desaturation specificity, chimeras of these two enzymes were tested for their fatty acid production in an engineered Yarrowia lipolytica strain. It reveals that transmembrane helices are not involved in the hydroxylation/desaturation specificity whereas all cytosolic domains have an impact on it. Especially, replacing the CpFAH12 cytosolic part near the second histidine-box by the corresponding CpFAD2 part annihilates all hydroxylation activity. Further mutagenesis experiments within this domain identified isoleucine 198 as the crucial element for the hydroxylation activity of CpFAH12. Monofunctional variants performing the only desaturation were obtained when this position was exchanged by the threonine of CpFAD2. Saturation mutagenesis at this position showed modulation in the hydroxylation/desaturation specificity in the different variants. The WT enzyme was demonstrated as the most efficient for ricinoleic acid production and some variants showed a better desaturation activity. A model based on the recently discovered membrane desaturase structures indicate that these changes in specificity are more likely due to modifications in the di-iron center geometry rather than changes in the substrate binding mode.

Production of Medium Chain Fatty Acids by Yarrowia lipolytica: Combining Molecular Design and TALEN to Engineer the Fatty Acid Synthase.

Yarrowia lipolytica is a promising organism for the production of lipids of biotechnological interest and particularly for biofuel. In this study, we engineered the key enzyme involved in lipid biosynthesis, the giant multifunctional fatty acid synthase (FAS), to shorten chain length of the synthesized fatty acids. Taking as starting point that the ketoacyl synthase (KS) domain of Yarrowia lipolytica FAS is directly involved in chain length specificity, we used molecular modeling to investigate molecular recognition of palmitic acid (C16 fatty acid) by the KS. This enabled to point out the key role of an isoleucine residue, I1220, from the fatty acid binding site, which could be targeted by mutagenesis. To address this challenge, TALEN (transcription activator-like effector nucleases)-based genome editing technology was applied for the first time to Yarrowia lipolytica and proved to be very efficient for inducing targeted genome modifications. Among the generated FAS mutants, those having a bulky aromatic amino acid residue in place of the native isoleucine at position 1220 led to a significant increase of myristic acid (C14) production compared to parental wild-type KS. Particularly, the best performing mutant, I1220W, accumulates C14 at a level of 11.6% total fatty acids. Overall, this work illustrates how a combination of molecular modeling and genome-editing technology can offer novel opportunities to rationally engineer complex systems for synthetic biology.

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics.

The accurate prediction of the structure and dynamics of DNA remains a major challenge in computational biology due to the dearth of precise experimental information on DNA free in solution and limitations in the DNA force-fields underpinning the simulations. A new generation of force-fields has been developed to better represent the sequence-dependent B-DNA intrinsic mechanics, in particular with respect to the BI ↔ BII backbone equilibrium, which is essential to understand the B-DNA properties. Here, the performance of MD simulations with the newly updated force-fields Parmbsc0εζOLI and CHARMM36 was tested against a large ensemble of recent NMR data collected on four DNA dodecamers involved in nucleosome positioning. We find impressive progress towards a coherent, realistic representation of B-DNA in solution, despite residual shortcomings. This improved representation allows new and deeper interpretation of the experimental observables, including regarding the behavior of facing phosphate groups in complementary dinucleotides, and their modulation by the sequence. It also provides the opportunity to extensively revisit and refine the coupling between backbone states and inter base pair parameters, which emerges as a common theme across all the complementary dinucleotides. In sum, the global agreement between simulations and experiment reveals new aspects of intrinsic DNA mechanics, a key component of DNA-protein recognition.

Energetic and molecular water permeation mechanisms of the human red blood cell urea transporter B.

Urea transporter B (UT-B) is a passive membrane channel that facilitates highly efficient permeation of urea. In red blood cells (RBC), while the major function of UT-B is to transport urea, it is assumed that this protein is able to conduct water. Here, we have revisited this last issue by studying RBCs and ghosts from human variants with defects of aquaporin 1 (AQP1) or UT-B. We found that UT-B's osmotic water unit permeability (pfunit) is similar to that of AQP1. The determination of diffusional permeability coefficient (Pd) allowed the calculation of the Pf/Pd ratio, which is consistent with a single-file water transport. Molecular dynamic simulations of water conduction through human UT-B confirmed the experimental finding. From these results, we propose an atomistic description of water-protein interactions involved in this permeation. Inside the UT-B pore, five water molecules were found to form a single-file and move rapidly along a channel by hydrogen bond exchange involving two critical threonines. We further show that the energy barrier for water located in the central region coincides with a water dipole reorientation, which can be related to the proton exclusion observed experimentally. In conclusion, our results indicate that UT-B should be considered as a new member of the water channel family.

Simulating DNA by molecular dynamics: aims, methods, and validation.

The structure and dynamics of the B-DNA double helix involves subtle sequence-dependent effects which are decisive for its function, but difficult to characterize. These structural and dynamic effects can be addressed by simulations of DNA sequences in explicit solvent. Here, we present and discuss the state-of-art of B-DNA molecular dynamics simulations with the major force fields in use today. We explain why a critical analysis of the MD trajectories is required to assess their reliability, and estimate the value and limitations of these models. Overall, simulations of DNA bear great promise towards deciphering the structural and physical subtleties of this biopolymer, where much remains to be understood.

Mg2+ in the major groove modulates B-DNA structure and dynamics.

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.

How cations can assist DNase I in DNA binding and hydrolysis.

DNase I requires Ca²+ and Mg²+ for hydrolyzing double-stranded DNA. However, the number and the location of DNase I ion-binding sites remain unclear, as well as the role of these counter-ions. Using molecular dynamics simulations, we show that bovine pancreatic (bp) DNase I contains four ion-binding pockets. Two of them strongly bind Ca²+ while the other two sites coordinate Mg²+. These theoretical results are strongly supported by revisiting crystallographic structures that contain bpDNase I. One Ca²+ stabilizes the functional DNase I structure. The presence of Mg²+ in close vicinity to the catalytic pocket of bpDNase I reinforces the idea of a cation-assisted hydrolytic mechanism. Importantly, Poisson-Boltzmann-type electrostatic potential calculations demonstrate that the divalent cations collectively control the electrostatic fit between bpDNase I and DNA. These results improve our understanding of the essential role of cations in the biological function of bpDNase I. The high degree of conservation of the amino acids involved in the identified cation-binding sites across DNase I and DNase I-like proteins from various species suggests that our findings generally apply to all DNase I-DNA interactions.

Multi-resolution approach for interactively locating functionally linked ion binding sites by steering small molecules into electrostatic potential maps using a haptic device.

Metal ions drive important parts of biology, yet it remains experimentally challenging to locate their binding sites. Here we present an innovative computational approach. We use interactive steering of charged ions or small molecules in an electrostatic potential map in order to identify potential binding sites. The user interacts with a haptic device and experiences tactile feedback related to the strength of binding at a given site. The potential field is the first level of resolution used in this model. Any type of potential field can be used, implicitly taking into account various conditions such as ionic strength, dielectric constants or the presence of a membrane. At a second level, we represent the accessibility of all binding sites by modelling the shape of the target macromolecule via non-bonded van der Waals interactions between its static atomic or coarse-grained structure and the probe molecule(s). The third independent level concerns the representation of the molecular probe itself. Ion selectivity can be assessed by using multiple interacting ions as probes. This method was successfully applied to the DNase I enzyme, where we recently identified two new cation binding sites by computationally expensive extended molecular dynamics simulations.