Mutational analysis in Corynebacterium stationis MFS transporters for improving nucleotide bioproduction.

Product secretion from an engineered cell can be advantageous for microbial cell factories. Extensive work on nucleotide manufacturing, one of the most successful microbial fermentation processes, has enabled Corynebacterium stationis to transport nucleotides outside the cell by random mutagenesis; however, the underlying mechanism has not been elucidated, hindering its applications in transporter engineering. Herein, we report the nucleotide-exporting major facilitator superfamily (MFS) transporter from the C. stationis genome and its hyperactive mutation at the G64 residue. Structural estimation and molecular dynamics simulations suggested that the activity of this transporter improved via two mechanisms: (1) enhancing interactions between transmembrane helices through the conserved "RxxQG" motif along with substrate binding and (2) trapping substrate-interacting residue for easier release from the cavity. Our results provide novel insights into how MFS transporters change their conformation from inward- to outward-facing states upon substrate binding to facilitate efflux and can contribute to the development of rational design approaches for efflux improvements in microbial cell factories. KEYPOINTS: • An MFS transporter from C. stationis genome and its mutation at residue G64 were assessed • It enhanced the transporter activity by strengthening transmembrane helix interactions and trapped substrate-interacting residues • Our results contribute to rational design approach development for efflux improvement.

Determining Factor of the Quantum Yield of the Cyclization Reaction via Triplet States for Dye-Attached Diarylethene.

The difference in the quantum yields of the cyclization reaction of two isomers of dye-attached diarylethene via the triplet state observed in the experiment [J. Phys. Chem. C 2009, 113, 11623-11627] was theoretically examined by quantum chemical calculations. By evaluating the spin-orbit couplings, we found that the ratio of the rate constants from the S1 state to the T2 state between two isomers agreed with that of the experimental quantum yield of the cyclization reaction. The difference in the spin-orbit couplings is due to the difference in the delocalization of the orbitals between diarylethene and dye. We further found that after the intersystem crossing took place the cyclization reaction via the triplet state occurred in the experiment due to the low energy barrier (∼10 kcal/mol) for the reaction.

Cupid and Psyche system for the diagnosis and treatment of advanced cancer.

In advanced cancer patients, malignant cells invade and disseminate within normal cells and develop resistance to therapy with additional genetic mutations, which makes radical cure very difficult. Precision medicine against advanced cancer is hampered by the lack of systems aimed at multiple target molecules within multiple loci. Here, we report the development of a versatile diagnostic and therapeutic system for advanced cancer, named the Cupid and Psyche system. Based on the strong non-covalent interaction of streptavidin and biotin, a low immunogenic mutated streptavidin, Cupid, and a modified artificial biotin, Psyche, have been designed. Cupid can be fused with various single-chain variable fragment antibodies and forms tetramer to recognize cancer cells precisely. Psyche can be conjugated to a wide range of diagnostic and therapeutic agents against malignant cells. The Cupid and Psyche system can be used in pre-targeting therapy as well as photo-immunotherapy effectively in animal models supporting the concept of a system for precision medicine for multiple targets within multiple loci.

Epiregulin Recognition Mechanisms by Anti-epiregulin Antibody 9E5: STRUCTURAL, FUNCTIONAL, AND MOLECULAR DYNAMICS SIMULATION ANALYSES.

Epiregulin (EPR) is a ligand of the epidermal growth factor (EGF) family that upon binding to its epidermal growth factor receptor (EGFR) stimulates proliferative signaling, especially in colon cancer cells. Here, we describe the three-dimensional structure of the EPR antibody (the 9E5(Fab) fragment) in the presence and absence of EPR. Among the six complementarity-determining regions (CDRs), CDR1-3 in the light chain and CDR2 in the heavy chain predominantly recognize EPR. In particular, CDR3 in the heavy chain dramatically moves with cis-trans isomerization of Pro(103). A molecular dynamics simulation and mutational analyses revealed that Arg(40) in EPR is a key residue for the specific binding of 9E5 IgG. From isothermal titration calorimetry analysis, the dissociation constant was determined to be 6.5 nm. Surface plasmon resonance analysis revealed that the dissociation rate of 9E5 IgG is extremely slow. The superimposed structure of 9E5(Fab)·EPR on the known complex structure of EGF·EGFR showed that the 9E5(Fab) paratope overlaps with Domains I and III on the EGFR, which reveals that the 9E5(Fab)·EPR complex could not bind to the EGFR. The 9E5 antibody will also be useful in medicine as a neutralizing antibody specific for colon cancer.

Two distinct amyloid beta-protein (Abeta) assembly pathways leading to oligomers and fibrils identified by combined fluorescence correlation spectroscopy, morphology, and toxicity analyses.

Nonfibrillar assemblies of amyloid ß-protein (Aß) are considered to play primary roles in Alzheimer disease (AD). Elucidating the assembly pathways of these specific aggregates is essential for understanding disease pathogenesis and developing knowledge-based therapies. However, these assemblies cannot be monitored in vivo, and there has been no reliable in vitro monitoring method at low protein concentration. We have developed a highly sensitive in vitro monitoring method using fluorescence correlation spectroscopy (FCS) combined with transmission electron microscopy (TEM) and toxicity assays. Using Aß labeled at the N terminus or Lys(16), we uncovered two distinct assembly pathways. One leads to highly toxic 10-15-nm spherical Aß assemblies, termed amylospheroids (ASPDs). The other leads to fibrils. The first step in ASPD formation is trimerization. ASPDs of ∼330 kDa in mass form from these trimers after 5 h of slow rotation. Up to at least 24 h, ASPDs remain the dominant structures in assembly reactions. Neurotoxicity studies reveal that the most toxic ASPDs are ∼128 kDa (∼32-mers). In contrast, fibrillogenesis begins with dimer formation and then proceeds to formation of 15-40-nm spherical intermediates, from which fibrils originate after 15 h. Unlike ASPD formation, the Lys(16)-labeled peptide disturbed fibril formation because the Aß(16-20) region is critical for this final step. These differences in the assembly pathways clearly indicated that ASPDs are not fibril precursors. The method we have developed should facilitate identifying Aß assembly steps at which inhibition may be beneficial.

Engineered MATE multidrug transporters reveal two functionally distinct ion-coupling pathways in NorM from Vibrio cholerae.

Multidrug and toxic compound extrusion (MATE) transport proteins confer multidrug resistance on pathogenic microorganisms and affect pharmacokinetics in mammals. Our understanding of how MATE transporters work, has mostly relied on protein structures and MD simulations. However, the energetics of drug transport has not been studied in detail. Many MATE transporters utilise the electrochemical H+ or Na+ gradient to drive substrate efflux, but NorM-VC from Vibrio cholerae can utilise both forms of metabolic energy. To dissect the localisation and organisation of H+ and Na+ translocation pathways in NorM-VC we engineered chimaeric proteins in which the N-lobe of H+-coupled NorM-PS from Pseudomonas stutzeri is fused to the C-lobe of NorM-VC, and vice versa. Our findings in drug binding and transport experiments with chimaeric, mutant and wildtype transporters highlight the versatile nature of energy coupling in NorM-VC, which enables adaptation to fluctuating salinity levels in the natural habitat of V. cholerae.

A New Lipid Force Field (FUJI).

To explore inhomogeneous and anisotropic systems such as lipid bilayers, the Lennard-Jones particle mesh Ewald (LJ-PME) method has been applied without a conventional isotropic dispersion correction. As the popular AMBER and CHARMM lipid force fields were developed using a cutoff scheme, their lipid bilayers unacceptably shrink when using the LJ-PME method. In this study, a new all-atom lipid force field (FUJI) was developed on the basis of the AMBER force-field scheme including the Lipid14 van der Waals parameters. Point charges were calculated using the restrained electrostatic potentials of many lipid conformers. Further, torsion energy profiles were calculated using high-level ab initio molecular orbitals (LCCSD(T)/Aug-cc-pVTZ//LMP2/Aug-cc-pVTZ), following which the molecular mechanical dihedral parameters were derived through a fast Fourier transform. By incorporation of these parameters into a new lipid force field without fitting experimental data, the desired lipid characteristics such as the area per lipid and lateral diffusion coefficients were obtained through GROMACS molecular dynamics simulations using the LJ-PME method and virtual hydrogen sites. The calculated area per lipid and lateral diffusion coefficients showed satisfactory agreement with experimental data. Furthermore, the electron-density profiles along the membrane normal were calculated for pure lipid bilayers, and the resulting membrane thicknesses agreed well with the experimental values. As the new lipid force field is compatible with FUJI for protein and small molecules, the new FUJI force field will offer accurate modeling for complex systems consisting of various membrane proteins and lipids.

Efficient free energy calculation of water across lipid membranes.

An efficient free energy (FE) calculation of a water molecule to go across lipid membranes is presented. Both overlapping distribution and cavity insertion Widom methods are complementarily used. The former is useful for a dense region where water molecules are found, i.e., from the interfacial to bulk water region, while the latter works well in the low density region, i.e., the hydrocarbon region. Since both methods evaluate the excess chemical potential of water, the obtained FE profile is free from the fitting problem usually arisen when a FE difference method is used. A diphytanyl phosphatidylcholine bilayer is used for our test calculations. An excellent and fast convergence of the chemical potential is obtained when each method is applied for the appropriate region. The estimated FE barrier using the Ewald method for the electrostatic interaction is approximately 7.2 kcal/mol, which is higher than that using the interaction cutoff of 20 A by about 0.9 kcal/mol.

Initiation of prolyl cis-trans isomerisation in the CDR-H3 loop of an antibody in response to antigen binding.

Proline cis-trans isomerisation is a regulatory mechanism used in a range of biological processes, and is related to various diseases such as Alzheimers disease and cancer. However, the details of the exact molecular mechanism by which it occurs are not known. Using X-ray crystallography, proline isomerisation has been shown to occur following formation of an antigen-antibody complex between the target epiregulin (EPR) and the antibody 9E5, at proline (Pro103), located in the third complementarity-determining region (CDR) of the heavy chain of 9E5. To obtain an accurate description of the pathway involved in cis-trans isomerisation in this system, we performed ten independent long molecular dynamics (MD) simulations starting at a stable transient bound structure obtained from many short binding MD simulations. As a result, we were able to describe the process by which cis-trans isomerisation is initiated, and suggest a catalysis mechanism for cis-trans isomerization in this antigen-antibody system. We found that Asp102, which is immediately adjacent to Pro103, rotates while changing its interacting partner residues in the light chain of 9E5, and at the same time EPR polar residues help to stabilise the intermediate states in the isomerisation process by interacting strongly with Asp102.

Beta1-adrenoceptor-mediated relaxation with isoprenaline and the role of MaxiK channels in guinea-pig esophageal smooth muscle.

The possible functional coupling between beta1-adrenoceptor and MaxiK channels which results in smooth muscle relaxation was examined in the guinea-pig esophageal muscularis mucosae. Isoprenaline-elicited relaxation of esophageal smooth muscle was confirmed to be mediated through beta1-adrenoceptors as the response was competitively antagonized by a beta1-selective antagonist atenolol with a pA2 value of 7.01. Iberiotoxin (IbTx, 10(-7) M), a selective MaxiK channel inhibitor, substantially diminished the relaxant response to isoprenaline. The extent of the MaxiK channel contribution to the relaxant response was 15-40% of the control response when estimated as the E50%-Emax responses to isoprenaline. The relaxation to isoprenaline was also attenuated by high-KCl (80 mM) to the same degree as the relaxant response generated in the presence of IbTx, and thus the estimated extent of the K+ channel contribution was 10-40%. These findings indicate that beta1-adrenoceptors are substantially coupled with MaxiK channels to produce relaxation of esophageal smooth muscle in the guinea-pig. Although MaxiK channels account for the contribution of K+ channels to the beta1-adrenoceptor-mediated relaxation in this smooth muscle preparation, their contribution seems to be less when compared to the beta2-adrenoceptor-mediated relaxation of tracheal smooth muscle.

π-Conjugation of two nitronyl nitroxides-attached diarylethenes.

Unusual blue shift of the absorption maxima of two nitronyl nitroxide attached diarylethene through phenyl units (DAE-phe-NN) with increasing number of phenyl units is examined by time dependent density functional theory (TDDFT). The extended π-conjugation between nitronyl nitroxide and diarylethene is rather suppressed by the bridge phenyl units. In comparison, the red shift found in two nitronyl nitroxide attached diarylethenes through thiophene units (DAE-thio-NN) with increasing number of thiophene units is due to the longer π-conjugation induced by smaller dihedral angles between diarylethene and bridge and between bridges.

Molecular dynamics simulation of an archaeal lipid bilayer with sodium chloride.

We have performed molecular dynamics simulations of a bilayer formed by the synthetic archaeal lipid, diphytanyl phosphatidylcholine, in NaCl electrolyte solution at four different concentrations (0-4 M) to investigate how structural and dynamic properties of the model archaeal membrane are changed due to the ionic strength of the solution. The archaeal lipid bilayer shows minor changes in their physical properties, indicating the unusual high stability of the membrane against salt, though small reductions of molecular area and lateral diffusion of the lipid are detected at the highest electrolyte concentration of 4 M. Sodium ions penetrate to the ether-rich region, where the ions are likely bound to the ether oxygen in the sn-1 chain rather than to that in the sn-2 chain. The observed salt bridges among two or three neighboring lipids account for the small reduction in the molecular area. The bound ions together with the counter (chloride) ions give rise to a diffusive electric double layer; as a result, the membrane dipole potential is slightly increased with increasing NaCl concentration.

Molecular dynamics study of bipolar tetraether lipid membranes.

Membranes composed of bipolar tetraether lipids have been studied by a series of 25-ns molecular dynamics simulations to understand the microscopic structure and dynamics as well as membrane area elasticity. By comparing macrocyclic and acyclic tetraether and diether archaeal lipids, the effect of tail linkage of the two phytanyl-chained lipids on the membrane properties is elucidated. Tetraether lipids show smaller molecular area and lateral mobility. For the latter, calculated diffusion coefficients are indeed one order-of-magnitude smaller than that of the diether lipid. These two tetraether membranes are alike in many physical properties except for membrane area elasticity. The macrocyclic tetraether membrane shows a higher elastic area expansion modulus than its acyclic counterpart by a factor of three. Free energy profiles of a water molecule crossing the membranes show no major difference in barrier height; however, a significant difference is observed near the membrane center due to the lack of the slip-plane in tetraether membranes.

Comparative molecular dynamics study of ether- and ester-linked phospholipid bilayers.

The lipid membranes found in archaea have high bilayer stability and low permeability. The molecular structure of their constituent lipids is characterized by ether-linked, branched hydrophobic chains, whereas the conventional lipids obtained from eukaryotic or eubacterial sources have ester linked straight chains. In order to elucidate the influence of the ether linkage, instead of an ester one, on the physical properties of the lipid bilayers, we have carried out comparative 10 ns molecular dynamics simulations of diphytanyl phosphatidylcholine (ether-DPhPC) and diphytanoyl phosphatidylcholine (ester-DPhPC) bilayers in water, respectively. We analyze bilayer structures, hydration of the lipids, membrane dipole potentials, and free energy profiles of water and oxygen across the bilayers. We observe that the membrane dipole potential for the ether-DPhPC bilayer, which arises mainly from the ether linkage, is about half of that of the ester-DPhPC. The calculated free energy barrier for a water molecule in the ether-DPhPC bilayer system is slightly higher than that in the ester-DPhPC counterpart, which is in accord with experimental data.