Optimized particle-mesh Ewald/multiple-time step integration for molecular dynamics simulations

Comparison of two peroxidases with high potential for biotechnology applications - HRP vs. APEX2

Peroxidases are essential elements in many biotechnological applications. An especially interesting concept involves split enzymes, where the enzyme is separated into two smaller and inactive proteins that can dimerize into a fully active enzyme. Such split forms were developed for the horseradish peroxidase (HRP) and ascorbate peroxidase (APX) already. Both peroxidases have a high potential for biotechnology applications. In the present study, we performed biophysical comparisons of these two peroxidases and their split analogues. The active site availability is similar for all four structures. The split enzymes are comparable in stability with their native analogues, meaning that they can be used for further biotechnology applications. Also, the tertiary structures of the two peroxidases are similar. However, differences that might help in choosing one system over another for biotechnology applications were noticed. The main difference between the two systems is glycosylation which is not present in the case of APX/sAPEX2, while it has a high impact on the HRP/sHRP stability. Further differences are calcium ions and cysteine bridges that are present only in the case of HRP/sHRP. Finally, computational results identified sAPEX2 as the systems with the smallest structural variations during molecular dynamics simulations showing its dominant stability comparing to other simulated proteins. Taken all together, the sAPEX2 system has a high potential for biotechnological applications due to the lack of glycans and cysteines, as well as due to high stability.

Repurposing FDA-approved drugs for combating tigecycline resistance in Acinetobacter baumannii: in silico screening against BaeR protein

Acinetobacter baumannii is becoming a gravely threatening nosocomial infection with a higher mortality rate. The present study targets the BaeR protein that mediates resistance to tigecycline antibiotics. The BaeR protein, along with the aid of BaeS, senses the incoming antibiotics and stimulates the expression of resistance proteins. These resistance proteins efflux the antibiotics and protect the cells from its effect. The main goal of the current study is to determine potential inhibitors from already existing FDA-approved drugs that could mitigate the BaeR protein. A range of in silico approaches, including molecular dynamics, virtual screening, SIFT analysis, ADMET, DFT, MM/GBSA, MMPBSA and per residue interaction analysis, were performed to identify inhibitors against this protein. The screening of FDA-approved compounds against the BaeR protein yielded 620 compounds. These compounds were clustered by SIFT to distinguish related compounds, it resulted in 20 different clusters. The top five clusters that can accommodate the binding site with better interaction and score by fulfilling all criteria were selected. The DFT analysis showed a smaller energy gap among all the compounds, indicating the ability of the compound to form firm interactions. All the compounds showed less binding free energy in both MM/GBSA and MM/PBSA analyses. The compounds were observed to be stable throughout the simulation. The per-residue interaction analysis confirmed that interactions with binding site residues were stable throughout the simulation. As a result of the study, four compounds, namely ZINC000003801919, DB01203, DB11217 and ZINC0000000056652, were identified as efficient candidates to deal with antimicrobial resistance in A. baumannii.

Direct mapping of tyrosine sulfation states in native peptides by nanopore

Sulfation is considered the most prevalent post-translational modification (PTM) on tyrosine; however, its importance is frequently undervalued due to difficulties in direct and unambiguous determination from phosphorylation. Here we present a sequence-independent strategy to directly map and quantify the tyrosine sulfation states in universal native peptides using an engineered protein nanopore. Molecular dynamics simulations and nanopore mutations reveal specific interactions between tyrosine sulfation and the engineered nanopore, dominating identification across diverse peptide sequences. We show a nanopore framework to discover tyrosine sulfation in unknown peptide fragments digested from a native protein and determine the sequence of the sulfated fragment based on current blockade enhancement induced by sulfation. Moreover, our method allows direct observation of peptide sulfation in ultra-low abundance, down to 1%, and distinguishes it from isobaric phosphorylation. This sequence-independent strategy suggests the potential of nanopore to explore specific PTMs in real-life samples and at the omics level.

Cation-induced intramolecular coil-to-globule transition in poly(ADP-ribose)

Poly(ADP-ribose) (PAR), a non-canonical nucleic acid, is essential for DNA/RNA metabolism and protein condensation, and its dysregulation is linked to cancer and neurodegeneration. However, key structural insights into PAR's functions remain largely uncharacterized, hindered by the challenges in synthesizing and characterizing PAR, which are attributed to its length heterogeneity. A central issue is how PAR, comprised solely of ADP-ribose units, attains specificity in its binding and condensing proteins based on chain length. Here, we integrate molecular dynamics simulations with small-angle X-ray scattering to analyze PAR structures. We identify diverse structural ensembles of PAR that fall into distinct subclasses and reveal distinct compaction of two different lengths of PAR upon the addition of small amounts of Mg2+ ions. Unlike PAR15, PAR22 forms ADP-ribose bundles via local intramolecular coil-to-globule transitions. Understanding these length-dependent structural changes could be central to deciphering the specific biological functions of PAR.

Molecular interactions of the Omicron, Kappa, and Delta SARS-CoV-2 spike proteins with quantum dots of graphene oxide

CONTEXT

The Omicron, Kappa, and Delta variants are different strains of the SARS-CoV-2 virus. Graphene oxide quantum dots (GOQDs) represent a burgeoning class of oxygen-enriched, zero-dimensional materials characterized by their sub-20-nm dimensions. Exhibiting pronounced quantum confinement and edge effects, GOQDs manifest exceptional physical-chemical attributes. This study delves into the potential of graphene oxide quantum dots, elucidating their inherent properties pertinent to the surface structures of SARS-CoV-2, employing an integrated computational approach for the repositioning of inhibitory agents.

METHODS

Following rigorous adjustment tests, a spectrum of divergent bonding conformations emerged, with particular emphasis placed on identifying the conformation exhibiting optimal adjustment scores and interactions. The investigation employed molecular docking simulations integrating affinity energy evaluations, electrostatic potential clouds, molecular dynamics encompassing average square root calculations, and the computation of Gibbs-free energy. These values quantify the strength of interaction between GOQDs and SARS-CoV-2 spike protein variants. The receptor structures were optimized using the CHARM-GUI server employing force field AMBERFF14SB. The algorithm embedded in CHARMM offers an efficient interpolation scheme and automatic step size selection, enhancing the efficiency of the optimization process. The 3D structures of the ligands are constructed and optimized with density functional theory (DFT) method based on the most stable conformer of each binder. Autodock Vina Software (ADV) was utilized, where essential parameters were specified. Electrostatic potential maps (MEPs) provide a visual depiction of molecules' charge distributions and related properties. After this, molecular dynamics simulations employing the CHARM36 force field in Gromacs 2022.2 were conducted to investigate GOs' interactions with surface macromolecules of SARS-CoV-2 in an explicit aqueous environment. Furthermore, our investigation suggests that lower values indicate stronger binding. Notably, GO-E consistently showed the most negative values across interactions with different variants, suggesting a higher affinity compared to other GOQDs (GO-A to GO-D).

An N-terminal acidic β-sheet domain is responsible for the metal-accumulation properties of amyloid-β protofibrils: a molecular dynamics study

The influence of metal ions on the structure of amyloid- β (Aβ) protofibril models was studied through molecular dynamics to explore the molecular mechanisms underlying metal-induced Aβ aggregation relevant in Alzheimer's disease (AD). The models included 36-, 48-, and 188-mers of the Aβ42 sequence and two disease-modifying variants. Primary structural effects were observed at the N-terminal domain, as it became susceptible to the presence of cations. Specially when β-sheets predominate, this motif orients N-terminal acidic residues toward one single face of the β-sheet, resulting in the formation of an acidic region that attracts cations from the media and promotes the folding of the N-terminal region, with implications in amyloid aggregation. The molecular phenotype of the protofibril models based on Aβ variants shows that the AD-causative D7N mutation promotes the formation of N-terminal β-sheets and accumulates more Zn2+, in contrast to the non-amyloidogenic rodent sequence that hinders the β-sheets and is more selective for Na+ over Zn2+ cations. It is proposed that forming an acidic β-sheet domain and accumulating cations is a plausible molecular mechanism connecting the elevated affinity and concentration of metals in Aβ fibrils to their high content of β-sheet structure at the N-terminal sequence.

The structure and physical properties of a packaged bacteriophage particle

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

Length-dependent Intramolecular Coil-to-Globule Transition in Poly(ADP-ribose) Induced by Cations

Poly(ADP-ribose) (PAR), as part of a post-translational modification, serves as a flexible scaffold for noncovalent protein binding. Such binding is influenced by PAR chain length through a mechanism yet to be elucidated. Structural insights have been elusive, partly due to the difficulties associated with synthesizing PAR chains of defined lengths. Here, we employ an integrated approach combining molecular dynamics (MD) simulations with small-angle X-ray scattering (SAXS) experiments, enabling us to identify highly heterogeneous ensembles of PAR conformers at two different, physiologically relevant lengths: PAR 15 and PAR 22 . Our findings reveal that numerous factors including backbone conformation, base stacking, and chain length contribute to determining the structural ensembles. We also observe length-dependent compaction of PAR upon the addition of small amounts of Mg 2+ ions, with the 22-mer exhibiting ADP-ribose bundles formed through local intramolecular coil-to-globule transitions. This study illuminates how such bundling could be instrumental in deciphering the length-dependent action of PAR.

GRAPHICAL ABSTRACT

Solvent-Induced Lag Phase during the Formation of Lysozyme Amyloid Fibrils Triggered by Sodium Dodecyl Sulfate: Biophysical Experimental and In Silico Study of Solvent Effects

Amyloid aggregates arise from either the partial or complete loss of the native protein structure or the inability of proteins to attain their native conformation. These aggregates have been linked to several diseases, including Alzheimer's, Parkinson's, and lysozyme amyloidosis. A comprehensive dataset was recently reported, demonstrating the critical role of the protein's surrounding environment in amyloid formation. In this study, we investigated the formation of lysozyme amyloid fibrils induced by sodium dodecyl sulfate (SDS) and the effect of solvents in the medium. Experimental data obtained through fluorescence spectroscopy revealed a notable lag phase in amyloid formation when acetone solution was present. This finding suggested that the presence of acetone in the reaction medium created an unfavorable microenvironment for amyloid fibril formation and impeded the organization of the denatured protein into the fibril form. The in silico data provided insights into the molecular mechanism of the interaction between acetone molecules and the lysozyme protofibril, once acetone presented the best experimental results. It was observed that the lysozyme protofibril became highly unstable in the presence of acetone, leading to the complete loss of its β-sheet conformation and resulting in an open structure. Furthermore, the solvation layer of the protofibril in acetone solution was significantly reduced compared to that in other solvents, resulting in fewer hydrogen bonds. Consequently, the presence of acetone facilitated the exposure of the hydrophobic portion of the protofibril, precluding the amyloid fibril formation. In summary, our study underscores the pivotal role the surrounding environment plays in influencing amyloid formation.

FA Sliding as the Mechanism for the ANT1-Mediated Fatty Acid Anion Transport in Lipid Bilayers

Mitochondrial adenine nucleotide translocase (ANT) exchanges ADP for ATP to maintain energy production in the cell. Its protonophoric function in the presence of long-chain fatty acids (FA) is also recognized. Our previous results imply that proton/FA transport can be best described with the FA cycling model, in which protonated FA transports the proton to the mitochondrial matrix. The mechanism by which ANT1 transports FA anions back to the intermembrane space remains unclear. Using a combined approach involving measurements of the current through the planar lipid bilayers reconstituted with ANT1, site-directed mutagenesis and molecular dynamics simulations, we show that the FA anion is first attracted by positively charged arginines or lysines on the matrix side of ANT1 before moving along the positively charged protein-lipid interface and binding to R79, where it is protonated. We show that R79 is also critical for the competitive binding of ANT1 substrates (ADP and ATP) and inhibitors (carboxyatractyloside and bongkrekic acid). The binding sites are well conserved in mitochondrial SLC25 members, suggesting a general mechanism for transporting FA anions across the inner mitochondrial membrane.

Protein nanopore reveals the renin-angiotensin system crosstalk with single-amino-acid resolution

The discovery of crosstalk effects on the renin-angiotensin system (RAS) is limited by the lack of approaches to quantitatively monitor, in real time, multiple components with subtle differences and short half-lives. Here we report a nanopore framework to quantitatively determine the effect of the hidden crosstalk between angiotensin-converting enzyme (ACE) and angiotensin-converting enzyme 2 (ACE2) on RAS. By developing an engineered aerolysin nanopore capable of single-amino-acid resolution, we show that the ACE can be selectively inhibited by ACE2 to prevent cleavage of angiotensin I, even when the concentration of ACE is more than 30-fold higher than that of ACE2. We also show that the activity of ACE2 for cleaving angiotensin peptides is clearly suppressed by the spike protein of SARS-CoV-2. This leads to the relaxation of ACE and the increased probability of accumulation of the principal effector angiotensin II. The spike protein of the SARS-CoV-2 Delta variant is demonstrated to have a much greater impact on the crosstalk than the wild type.

Molecular dynamics simulations depict structural motions of the whole human aryl hydrocarbon receptor influencing its binding of ligands and HSP90

The aryl hydrocarbon receptor (AhR) has broad biological functions when its ligands activate it; the non-binding interactions with AhR have not been fully elucidated due to the absence of a complete tridimensional (3D) structure. Therefore, utilization of the whole 3D structure from Homo sapiens AhR by in silico studies will allow us to better study and analyze the binding mode of its full and partial agonists, and antagonists, as well as its interaction with the HSP90 chaperone. The 3D AhR structure was obtained from I-TASSER and subjected to molecular dynamics (MD) simulations to obtain different structural conformations and determine the most populated AhR conformer by clustering analyses. The AhR-3D structures selected from MD simulations and those from clustering analyses were used to achieve docking studies with some of its ligands and protein-protein docking with HSP90. Once the AhR-3D structure was built, its Ramachandran maps and energy showed a well-qualified 3D model. MD simulations showed that the per-Arnt-Sim homology (PAS) PAS A, PAS B, and Q domains underwent conformational changes, identifying the conformation when agonists were binding also, and HSP90 was binding near the PAS A, PAS B, and Q domains. However, when antagonists are binding, HSP90 does not bind near the PAS A, PAS B, and Q domains. These studies show that the complex agonist-AhR-HSP90 can be formed, but this complex is not formed when an antagonist is binding. Knowing the conformations when the ligands bind to AHR and the behavior of HSP90 allows for an understanding of its activity.Communicated by Ramaswamy H. Sarma.

All-Atom Molecular Dynamics Simulations of Membrane-Spanning DNA Origami Nanopores

Building on the recent technological advances, all-atom molecular dynamics (MD) simulations have become an indispensable tool to study the molecular behavior at nanoscale. Molecular simulations have been used to characterize the structure, dynamics, and mechanical and electrical properties of DNA origami objects. In this chapter we describe a method to build all-atom model of lipid-spanning DNA origami nanopores and perform molecular dynamics simulations in explicit electrolyte solutions.

Quantitative characterization of the path of glucose diffusion facilitated by human glucose transporter 1

Glucose transporter GLUT1 is ubiquitously expressed in the human body from the red cells to the blood-brain barrier to the skeletal muscles. It is physiologically relevant to understand how GLUT1 facilitates diffusion of glucose across the cell membrane. It is also pathologically relevant because GLUT1 deficiency causes neurological disorders and anemia and because GLUT1 overexpression fuels the abnormal growth of cancer cells. This article presents a quantitative investigation of GLUT1 based on all-atom molecular-dynamics (MD) simulations of the transporter embedded in lipid bilayers of asymmetric inner-and-outer-leaflet lipid compositions, subject to asymmetric intra-and-extra-cellular environments. This is in contrast with the current literature of MD studies that have not considered both of the aforementioned asymmetries of the cell membrane. The equilibrium (unbiased) dynamics of GLUT1 shows that it can facilitate glucose diffusion across the cell membrane without undergoing large-scale conformational motions. The Gibbs free-energy profile, which is still lacking in the current literature of GLUT1, quantitatively characterizes the diffusion path of glucose from the periplasm, through an extracellular gate of GLUT1, on to the binding site, and off to the cytoplasm. This transport mechanism is validated by the experimental data that GLUT1 has low water-permeability, uptake-efflux symmetry, and 10 kcal/mol Arrhenius activation barrier around 37 °C.

A Computational–Experimental Investigation of the Molecular Mechanism of Interleukin-6-Piperine Interaction

Herein, we elucidate the biophysical aspects of the interaction of an important protein, Interleukin-6 (IL6), which is involved in cytokine storm syndrome, with a natural product with anti-inflammatory activity, piperine. Despite the role of piperine in the inhibition of the transcriptional protein NF-κB pathway responsible for activation of IL6 gene expression, there are no studies to the best of our knowledge regarding the characterisation of the molecular interaction of the IL6-piperine complex. In this context, the characterisation was performed with spectroscopic experiments aided by molecular modelling. Fluorescence spectroscopy alongside van’t Hoff analyses showed that the complexation event is a spontaneous process driven by non-specific interactions. Circular dichroism aided by molecular dynamics revealed that piperine caused local α-helix reduction. Molecular docking and molecular dynamics disclosed the microenvironment of interaction as non-polar amino acid residues. Although piperine has three available hydrogen bond acceptors, only one hydrogen-bond was formed during our simulation experiments, reinforcing the major role of non-specific interactions that we observed experimentally. Root mean square deviation (RMSD) and hydrodynamic radii revealed that the IL6-piperine complex was stable during 800 ns of simulation. Taken together, these results can support ongoing IL6 drug discovery efforts.

The Ability of Chlorophyll to Trap Carcinogen Aflatoxin B1: A Theoretical Approach

The coordination of one and two aflatoxin B1 (AFB1, a potent carcinogen) molecules with chlorophyll a (chl a) was studied at a theoretical level. Calculations were performed using the M06-2X method in conjunction with the 6-311G(d,p) basis set, in both gas and water phases. The molecular electrostatic potential map shows the chemical activity of various sites of the AFB1 and chl a molecules. The energy difference between molecular orbitals of AFB1 and chl a allowed for the establishment of an intermolecular interaction. A charge transfer from AFB1 to the central cation of chl a was shown. The energies of the optimized structures for chl a show two configurations, unfolded and folded, with a difference of 15.41 kcal/mol. Chl a appeared axially coordinated to the plane (α-down or β-up) of the porphyrin moiety, either with the oxygen atom of the ketonic group, or with the oxygen atom of the lactone moiety of AFB1. The complexes of maximum stability were chl a 1-α-E-AFB1 and chl a 2-β-E-AFB1, at -36.4 and -39.2 kcal/mol, respectively. Additionally, with two AFB1 molecules were chl a 1-D-2AFB1 and chl a 2-E-2AFB1, at -60.0 and -64.8 kcal/mol, respectively. Finally, biosorbents containing chlorophyll could improve AFB1 adsorption.

Single-molecule biophysics experiments in silico: Toward a physical model of a replisome

The interpretation of single-molecule experiments is frequently aided by computational modeling of biomolecular dynamics. The growth of computing power and ongoing validation of computational models suggest that it soon may be possible to replace some experiments outright with computational mimics. Here, we offer a blueprint for performing single-molecule studies in silico using a DNA-binding protein as a test bed. We demonstrate how atomistic simulations, typically limited to sub-millisecond durations and zeptoliter volumes, can guide development of a coarse-grained model for use in simulations that mimic single-molecule experiments. We apply the model to recapitulate, in silico, force-extension characterization of protein binding to single-stranded DNA and protein and DNA replacement assays, providing a detailed portrait of the underlying mechanics. Finally, we use the model to simulate the trombone loop of a replication fork, a large complex of proteins and DNA.

•

Coarse-grained model derived from all-atom simulation recapitulates experiments

•

Model reproduces the elastic response to force and exchange dynamics

•

Model reveals structure of intermediate states usually inaccessible to experiment

•

Model applied to viral replisome with trombone loop containing tens of SSB proteins

Discovery and Characterization of a Cryptic Secondary Binding Site in the Molecular Chaperone HSP70

Heat Shock Protein 70s (HSP70s) are key molecular chaperones that are overexpressed in many cancers and often associated with metastasis and poor prognosis. It has proven difficult to develop ATP-competitive, drug-like small molecule inhibitors of HSP70s due to the flexible and hydrophilic nature of the HSP70 ATP-binding site and its high affinity for endogenous nucleotides. The aim of this study was to explore the potential for the inhibition of HSP70 through alternative binding sites using fragment-based approaches. A surface plasmon resonance (SPR) fragment screen designed to detect secondary binding sites in HSP70 led to the identification by X-ray crystallography of a cryptic binding site in the nucleotide-binding domain (NBD) of HSP70 adjacent to the ATP-binding site. Fragment binding was confirmed and characterized as ATP-competitive using SPR and ligand-observed NMR methods. Molecular dynamics simulations were applied to understand the interactions with the protein upon ligand binding, and local secondary structure changes consistent with interconversion between the observed crystal structures with and without the cryptic pocket were detected. A virtual high-throughput screen (vHTS) against the cryptic pocket was conducted, and five compounds with diverse chemical scaffolds were confirmed to bind to HSP70 with micromolar affinity by SPR. These results identified and characterized a new targetable site on HSP70. While targeting HSP70 remains challenging, the new site may provide opportunities to develop allosteric ATP-competitive inhibitors with differentiated physicochemical properties from current series.

Unravelling the Interaction of Piperlongumine with the Nucleotide-Binding Domain of HSP70: A Spectroscopic and In Silico Study

Piperlongumine (PPL) is an alkaloid extracted from several pepper species that exhibits anti-inflammatory and anti-carcinogenic properties. Nevertheless, the molecular mode of action of PPL that confers such powerful pharmacological properties remains unknown. From this perspective, spectroscopic methods aided by computational modeling were employed to characterize the interaction between PPL and nucleotide-binding domain of heat shock protein 70 (NBD/HSP70), which is involved in the pathogenesis of several diseases. Steady-state fluorescence spectroscopy along with time-resolved fluorescence revealed the complex formation based on a static quenching mechanism. Van't Hoff analyses showed that the binding of PPL toward NBD is driven by equivalent contributions of entropic and enthalpic factors. Furthermore, IDF and Scatchard methods applied to fluorescence intensities determined two cooperative binding sites with K_b of (6.3 ± 0.2) × 10⁴ M^-1. Circular dichroism determined the thermal stability of the NBD domain and showed that PPL caused minor changes in the protein secondary structure. Computational simulations elucidated the microenvironment of these interactions, showing that the binding sites are composed mainly of polar amino acids and the predominant interaction of PPL with NBD is Van der Waals in nature.

Revisiting the Origin of Nanopore Current Blockage for Volume Difference Sensing at the Atomic Level

Changes in the nanopore ionic current during entry of a target molecule underlie the sensing capability and dominate the intensity and extent of applications of the nanopore approach. The volume exclusion model has been proposed and corrected to describe the nanopore current blockage. However, increasing evidence shows nonconformity with this model, suggesting that the ionic current within a nanopore should be entirely reconsidered. Here, we revisit the origin of nanopore current blockage from a theoretical perspective and propose that the noncovalent interactions between a nanopore and a target molecule affect the conductance of the solution inside the nanopore, leading to enhanced current blockage. Moreover, by considering the example of an aerolysin nanopore discriminating the cytosine DNA and methylcytosine DNA that differ by a single methyl group, we completely demonstrate, by nanopore experiments and molecular dynamics simulations, the essential nature of this noncovalent interaction for discrimination. Our conductance model suggests multiplicative effects of both volume exclusion and noncovalent interaction on the current blockage and provides a new strategy to achieve volume difference sensing at the atomic level with highly specific current events, which would promote the nanopore protein sequencing and its applications in real-life systems.

Optimized Hydrogen Mass Repartitioning Scheme Combined with Accurate Temperature/Pressure Evaluations for Thermodynamic and Kinetic Properties of Biological Systems

In recent years, molecular dynamics (MD) simulations with larger time steps have been performed with the hydrogen-mass-repartitioning (HMR) scheme, where the mass of each hydrogen atom is increased to reduce high-frequency motion while the mass of a non-hydrogen atom bonded to a hydrogen atom is decreased to keep the total molecular mass unchanged. Here, we optimize the scaling factors in HMR and combine them with previously developed accurate temperature/pressure evaluations. The heterogeneous HMR scaling factors are useful to avoid the structural instability of amino acid residues having a five- or six-membered ring in MD simulations with larger time steps. It also reproduces kinetic properties, namely translational and rotational diffusions, if the HMR scaling factors are applied to only solute molecules. The integration scheme is tested for biological systems that include soluble/membrane proteins and lipid bilayers for about 200 μs MD simulations in total and give consistent results in MD simulations with both a small time step of 2.0 fs and a large, multiple time step integration with time steps of 3.5 fs (for fast motions) and 7.0 fs (for slower motions). We also confirm that the multiple time step integration scheme used in this study provides more accurate energy conservations than the RESPA/C1 and is comparable to the RESPA/C2 in NAMD. In summary, the current integration scheme combining the optimized HMR with accurate temperature/pressure evaluations can provide stable and reliable MD trajectories with a larger time step, which are computationally more than 2-fold efficient compared to the conventional methods.

Protein Sumoylation Is Crucial for Phagocytosis in Entamoeba histolytica Trophozoites

Posttranslational modifications provide Entamoeba histolytica proteins the timing and signaling to intervene during different processes, such as phagocytosis. However, SUMOylation has not been studied in E. histolytica yet. Here, we characterized the E. histolytica SUMO gene, its product (EhSUMO), and the relevance of SUMOylation in phagocytosis. Our results indicated that EhSUMO has an extended N-terminus that differentiates SUMO from ubiquitin. It also presents the GG residues at the C-terminus and the ΨKXE/D binding motif, both involved in target protein contact. Additionally, the E. histolytica genome possesses the enzymes belonging to the SUMOylation-deSUMOylation machinery. Confocal microscopy assays disclosed a remarkable EhSUMO membrane activity with convoluted and changing structures in trophozoites during erythrophagocytosis. SUMOylated proteins appeared in pseudopodia, phagocytic channels, and around the adhered and ingested erythrocytes. Docking analysis predicted interaction of EhSUMO with EhADH (an ALIX family protein), and immunoprecipitation and immunofluorescence assays revealed that the association increased during phagocytosis; whereas the EhVps32 (a protein of the ESCRT-III complex)-EhSUMO interaction appeared stronger since basal conditions. In EhSUMO knocked-down trophozoites, the bizarre membranous structures disappeared, and EhSUMO interaction with EhADH and EhVps32 diminished. Our results evidenced the presence of a SUMO gene in E. histolytica and the SUMOylation relevance during phagocytosis. This is supported by bioinformatics screening of many other proteins of E. histolytica involved in phagocytosis, which present putative SUMOylation sites and the ΨKXE/D binding motif.

Analysis of the Relationship Between Metalloprotease-9 and Tau Protein in Alzheimer's Disease

BACKGROUND

Neurofibrillary tangles (NFTs) and amyloid plaques are the neuropathological hallmarks in brains with Alzheimer's disease (AD). Post-translational modifications of tau, such as phosphorylation and truncation, have been proposed as initiators in the assembly of the abnormal paired helical filaments that constitute the NFTs. Neurons and NFTs are sites of matrix metalloproteinases (MMPs).

OBJECTIVE

The aim of this study was to analyze the relationship of MMP-9 and tau protein in brain samples with AD.

METHODS

This study was performed on brain tissue samples from patients with early, moderate, and late AD. MMPs and tau levels were analyzed by western blot and gelatin-substrate zymography. Immunofluorescence techniques and confocal microscopy were used to analyze the presence of both proteins in NFTs. Further, molecular dynamics simulations (MDS) and protein-protein docking were conducted to predict interaction between MMP-9 and tau protein.

RESULTS

MMP-9 expression was greatest in moderate and late AD, whereas MMP-2 expression was only increased in late-stage AD. Interestingly, confocal microscopy revealed NFTs in which there was co-localization of MMP-9 and tau protein. MDS and protein-protein docking predictions indicate that a high-affinity complex can be formed between MMP-9 and full-length tau protein.

CONCLUSION

These observations provide preliminary evidence of an interaction between these two proteins. Post-translational modifications of tau protein, such as C-terminal truncation or phosphorylation of amino acid residues in the MMP-9 recognition site and conformational changes in the protein, such as folding of the N-terminal sequence over the three-repeat domain, could preclude the interaction between MMP-9 and tau protein during stages of NFT development.

Biomolecular modeling thrives in the age of technology

The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.

Determining the In-Plane Orientation and Binding Mode of Single Fluorescent Dyes in DNA Origami Structures

We present a technique to determine the orientation of single fluorophores attached to DNA origami structures based on two measurements. First, the orientation of the absorption transition dipole of the molecule is determined through a polarization-resolved excitation measurement. Second, the orientation of the DNA origami structure is obtained from a DNA-PAINT nanoscopy measurement. Both measurements are performed consecutively on a fluorescence wide-field microscope. We employed this approach to study the orientation of single ATTO 647N, ATTO 643, and Cy5 fluorophores covalently attached to a 2D rectangular DNA origami structure with different nanoenvironments, achieved by changing both the fluorophores' binding position and immediate vicinity. Our results show that when fluorophores are incorporated with additional space, for example, by omitting nucleotides in an elsewise double-stranded environment, they tend to stick to the DNA and to adopt a preferred orientation that depends more on the specific molecular environment than on the fluorophore type. With the aid of all-atom molecular dynamics simulations, we rationalized our observations and provide insight into the fluorophores' probable binding modes. We believe this work constitutes an important step toward manipulating the orientation of single fluorophores in DNA origami structures, which is vital for the development of more efficient and reproducible self-assembled nanophotonic devices.

Carbon nanotubes for rapid capturing of SARS-COV-2 virus: revealing a mechanistic aspect of binding based on computational studies

We investigate the binding interactions of synthesized multi-walled carbon nanotubes (MWCNTs) with SARS-CoV-2 virus. Two essential components of the SARS-CoV-2 structure i.e.6LU7 (main protease of SARS-CoV-2) and 6LZG (spike receptor-binding domain complexed with its receptor ACE2) were used for computational studies. MWCNTs of different morphologies (zigzag, armchair and chiral) were synthesized through a thermal chemical vapour deposition process as a function of pyrolysis temperature. A direct correlation between radius to volume ratio of the synthesized MWCNTs and the binding energies for all three (zigzag, armchair and chiral) conformations were observed in our computational studies. Our result suggests that MWCNTs interact with the active sites of the main protease along with the host angiotensin-converting enzyme2 (ACE2) receptors. Furthermore, it is also observed that MWCNTs have significant binding affinities towards SARS-CoV-2. However, the highest free binding energy of -87.09 kcal mol-1 with 6LZG were shown by the armchair MWCNTs with SARS-CoV-2 through the simulated molecular dynamic trajectories, which could alter the SARS-CoV-2 structure with higher accuracy. The radial distribution function also confirms the density variation as a function of distance from a reference particle of MWCNTs for the study of interparticle interactions of the MWCNT and SARS-CoV-2. Due to these interesting attributes, such MWCNTs could find potential application in personal protective equipment (PPE) and diagnostic kits.

Detailed Characterization of the Cooperative Binding of Piperine with Heat Shock Protein 70 by Molecular Biophysical Approaches

In this work, for the first time, details of the complex formed by heat shock protein 70 (HSP70) independent nucleotide binding domain (NBD) and piperine were characterized through experimental and computational molecular biophysical methods. Fluorescence spectroscopy results revealed positive cooperativity between the two binding sites. Circular dichroism identified secondary conformational changes. Molecular dynamics along with molecular mechanics Poisson Boltzmann surface area (MM/PBSA) reinforced the positive cooperativity, showing that the affinity of piperine for NBD increased when piperine occupied both binding sites instead of one. The spontaneity of the complexation was demonstrated through the Gibbs free energy (∆G < 0 kJ/mol) for different temperatures obtained experimentally by van’t Hoff analysis and computationally by umbrella sampling with the potential of mean force profile. Furthermore, the mean forces which drove the complexation were disclosed by van’t Hoff and MM/PBSA as being the non-specific interactions. In conclusion, the work revealed characteristics of NBD and piperine interaction, which may support further drug discover studies.

High-Fidelity Capture, Threading, and Infinite-Depth Sequencing of Single DNA Molecules with a Double-Nanopore System

Nanopore sequencing of nucleic acids has an illustrious history of innovations that eventually made commercial nanopore sequencing possible. Nevertheless, the present nanopore sequencing technology leaves much room for improvement, especially with respect to accuracy of raw reads and detection of nucleotide modifications. Double-nanopore sequencing-an approach where a DNA molecule is pulled back and forth by a tug-of-war of two nanopores-could potentially improve single-molecule read accuracy and modification detection by offering multiple reads of the same DNA fragment. One principle difficulty in realizing such a technology is threading single-stranded DNA through both nanopores. Here, we describe and demonstrate through simulations a nanofluidic system for loading and threading DNA strands through a double-nanopore setup with nearly 100% fidelity. The high-efficiency loading is realized by using hourglass-shaped side channels that not only deliver the molecules to the nanopore but also retain molecules that missed the nanopore at the first passage to attempt the nanopore capture again. The second nanopore capture is facilitated by an orthogonal microfluidic flow that unravels the molecule captured by the first nanopore and delivers it to the capture volume of the second nanopore. We demonstrate the potential utility of our double-nanopore system for DNA sequencing by simulating repeat back-and-forth motion-flossing-of a DNA strand through the double-nanopore system. We show that repeat exposure of the same DNA fragments to the nanopore sensing volume considerably increases accuracy of the nucleotide sequence determination and that correlated displacement of ssDNA through the two nanopores may facilitate recognition of homopolymer fragments.

The ENUF method-Ewald summation based on nonuniform fast Fourier transform: Implementation, parallelization, and application

Computer simulations of model systems are widely used to explore striking phenomena in promising applications spanning from physics, chemistry, biology, to materials science and engineering. The long range electrostatic interactions between charged particles constitute a prominent factor in determining structures and states of model systems. How to efficiently calculate electrostatic interactions in simulation systems subjected to partial or full periodic boundary conditions has been a grand challenging task. In the past decades, a large variety of computational schemes has been proposed, among which the Ewald summation method is the most reliable route to accurately deal with electrostatic interactions between charged particles in simulation systems. In addition, extensive efforts have been done to improve computational efficiencies of the Ewald summation based methods. Representative examples are approaches based on cutoffs, reaction fields, multi-poles, multi-grids, and particle-mesh schemes. We sketched an ENUF method, an abbreviation for the Ewald summation method based on the nonuniform fast Fourier transform technique, and have implemented this method in particle-based simulation packages to calculate electrostatic energies and forces at micro- and mesoscopic levels. Extensive computational studies of conformational properties of polyelectrolytes, dendrimer-membrane complexes, and ionic fluids demonstrated that the ENUF method and its derivatives conserve both energy and momentum to floating point accuracy, and exhibit a computational complexity of O N log N with optimal physical parameters. These ENUF based methods are attractive alternatives in molecular simulations where high accuracy and efficiency of simulation methods are needed to accelerate calculations of electrostatic interactions at extended spatiotemporal scales.

The Cytokine IL-1β and Piperine Complex Surveyed by Experimental and Computational Molecular Biophysics

The bioactive piperine, a compound found in some pepper species, has been widely studied because of its therapeutic properties that include the inhibition of an important inflammation pathway triggered by interleukin-1 beta (IL-1β). However, investigation into the molecular interactions between IL-1β and piperine is not reported in the literature. Here, we present for the first time the characterisation of the complex formed by IL-1β and piperine through experimental and computational molecular biophysical analyses. Fluorescence spectroscopy unveiled the presence of one binding site for piperine with an affinity constant of 14.3 × 10⁴ M⁻¹ at 298 K. The thermodynamic analysis indicated that the interaction with IL-1β was spontaneous (∆G = −25 kJ/mol) and, when split into enthalpic and entropic contributions, the latter was more significant. Circular dichroism spectroscopy showed that piperine did not affect IL-1β secondary structure (~2%) and therefore its stability. The set of experimental data parameterized the computational biophysical approach. Through molecular docking, the binding site micro-environment was revealed to be composed mostly by non-polar amino acids. Furthermore, molecular dynamics, along with umbrella sampling, are in agreement with the thermodynamic parameters obtained by fluorescence assays and showed that large protein movements are not present in IL-1β, corroborating the circular dichroism data.

Molecular Transport across the Ionic Liquid-Aqueous Electrolyte Interface in a MoS2 Nanopore

Nanopore sequencing of DNA has been enabled by the use of a biological enzyme to thread DNA through an engineered biological nanopore while recording the ionic current flowing through the nanopore. Efforts to realize a similar concept using a solid-state nanopore have been met with several technical challenges, one of which is the high speed of DNA translocation and the other the low ionic current contrast among individual nucleotides. A promising avenue to addressing both problems is using an ionic liquid to slow DNA translocation and a tiny nanopore in the MoS₂ membrane to distinguish individual nucleotides. The physical mechanisms enabling these technical advances have remained elusive. Here, we characterize the ion and DNA transport through the ionic liquid/aqueous electrolyte interface, with and without a MoS₂ nanopore, using the all-atom molecular dynamics method. We find that the partial miscibility of the ionic liquid and the aqueous electrolyte considerably alters the physics of the nanopore translocation process. Thus, the interface of the two phases generates a contact potential of 600 mV, the ionic current is dominated by the motion of ionic liquid molecules through the aqueous solution phase, and the DNA nucleotides exhibit preferential partitioning into the aqueous electrolyte, which leads to spontaneous transport of DNA polymers from the ionic liquid to the aqueous solution compartment in the absence of external voltage bias. The complex physics of the two-phase nanopore system offers a multitude of opportunities for extending the functionality of nanopore-sensing platforms.

Mechanistic Insights into the Effect of Ligands on Structural Stability and Selectivity of Sulfotransferase 2A1 (SULT2A1)

Cytosolic sulfotransferases (SULTs) acting as phase II metabolic enzymes can be used in the sulfonation of small molecules by transferring a sulfonate group from the unique co-factor 3'-phosphoadenosine 5'-phosphosulfate (PAPS) to the substrates. In the present study, molecular dynamics (MD) simulations and ensemble docking study were employed to theoretically characterize the mechanism for the effect of co-factor (PAP) and ligands (LCA, raloxifene, α-hydroxytamoxifen, ouabain, and 3'-phosphoadenylyl sulfate) on structural stability and selectivity of SULT2A1 from the perspective of the dynamic behavior of SULT2A1 structures. Structural stability and network analyses indicated that the cooperation between PAP and LCA may enhance the thermal stability and compact communication in enzymes. During the MD simulations, the obviously rigid region and inward displacement were detected in the active-site cap (loop16) of the conformation containing PAP, which may be responsible for the significant changes in substrate accessibility and catalytic activity. The smaller substrates such as LCA could bind stably to the active pocket in the presence of PAP. However, the substrates or inhibitors with a large spatial structure needed to bind to the open conformation (without PAP) prior to PAPS binding.

Automated sequence design of 2D wireframe DNA origami with honeycomb edges

Wireframe DNA origami has emerged as a powerful approach to fabricating nearly arbitrary 2D and 3D geometries at the nanometer-scale. Complex scaffold and staple routing needed to design wireframe DNA origami objects, however, render fully automated, geometry-based sequence design approaches essential for their synthesis. And wireframe DNA origami structural fidelity can be limited by wireframe edges that are composed only of one or two duplexes. Here we introduce a fully automated computational approach that programs 2D wireframe origami assemblies using honeycomb edges composed of six parallel duplexes. These wireframe assemblies show enhanced structural fidelity from electron microscopy-based measurement of programmed angles compared with identical geometries programmed using dual-duplex edges. Molecular dynamics provides additional theoretical support for the enhanced structural fidelity observed. Application of our top-down sequence design procedure to a variety of complex objects demonstrates its broad utility for programmable 2D nanoscale materials.

Wireframe DNA origami is a powerful approach to creating 2D and 3D geometries. Here the authors introduce an automated computational design approach that programs structures with high structural fidelity.

Crystallographic structure versus homology model: a case study of molecular dynamics simulation of human and zebrafish histone deacetylase 10

Histone deacetylase (HDAC) 10 has been implicated in the pathology of various cancers and neurodegenerative disorders, making the discovery of novel inhibitors of the isoform an important endeavor. However, the unavailability of crystallographic structure of human HDAC10 (hHDAC10) hinders structure-based drug design effort. Previously, we reported the homology modeled structure of human HDAC10 built using the crystallographic structure of Danio rerio (zebrafish) HDAC10 (zHDAC10) (Protein Data Bank (PDB) ID; 5TD7, released on 24 May 2017) as a template. Here, in continuation with our study, both hHDAC10 and zHDAC10, and their respective complexes with trichostatin A (TSA), quisinostat, and the native ligand (in 5TD7), 7-[(3-aminopropyl)amino]-1,1,1-trifluoroheptane-2,2-diol (PDB ID; FKS) were submitted to 100 ns-long unrestrained molecular dynamics (MD) simulations. Comparative analyses of the MD trajectories revealed that zHDAC10 and its complexes displayed higher stability than hHDAC10 and its corresponding complexes over time. Nonetheless, docking of active and inactive set molecules revealed that more reliable conformations of hHDAC10 could be obtained at an extended time period. This study may shed more light on the reliability of hHDAC10 modeled structure for use in selective inhibitor design.Communicated by Ramaswamy H. Sarma.

A theoretical and experimental approach to evaluate zein-calcium interaction in nixtamalization process

The possible interactions between α-zein and Ca²⁺ in nixtamalization process were analyzed from a multidisciplinary approach, considering the effect of these interactions on the thermal properties of the nixtamalized flour. SDS-PAGE under reducing and non-reducing conditions did not reveal differences between patterns of zeins from nixtamalized and control samples. However, analysis from affinity capillary electrophoresis indicated an increment in protein volume when calcium is added to zein extracted from nixtamalized flour. In addition, the binding constant for the zein-calcium interaction was calculated indicating a higher affinity for calcium by zein from nixtamalized samples. Molecular dynamics simulations indicated that the interaction α-zein-Ca²⁺ through C-ter was more favorable than Glu48. However, in excess of Ca²⁺ ions, each site could bind one calcium atom at the same time, confirming that aggregation of α-zein through calcium bridges is possible, expanding the technological applications of this protein.

EhRabB mobilises the EhCPADH complex through the actin cytoskeleton during phagocytosis of Entamoeba histolytica

Movement and phagocytosis are clue events in colonisation and invasion of tissues by Entamoeba histolytica, the protozoan causative of human amoebiasis. During phagocytosis, EhRab proteins interact with other functional molecules, conducting them to the precise cellular site. The gene encoding EhrabB is located in the complementary chain of the DNA fragment containing Ehcp112 and Ehadh genes, which encode for the proteins of the EhCPADH complex, involved in phagocytosis. This particular genetic organisation suggests that the three corresponding proteins may be functionally related. Here, we studied the relationship of EhRabB with EhCPADH and actin during phagocytosis. First, we obtained the EhRabB 3D structure to carry out docking analysis to predict the interaction sites involved in the EhRabB protein and the EhCPADH complex contact. By confocal microscopy, transmission electron microscopy, and immunoprecipitation assays, we revealed the interaction among these proteins when they move through different vesicles formed during phagocytosis. The role of the actin cytoskeleton in this event was also confirmed using Latrunculin A to interfere with actin polymerisation. This affected the movement of EhRabB and EhCPADH, as well as the rate of phagocytosis. Mutant trophozoites, silenced in EhrabB gene, evidenced the interaction of this molecule with EhCPADH and strengthened the role of actin during erythrophagocytosis.

Confining a Protein-Containing Water Nanodroplet inside Silica Nanochannels

Incorporation of biological systems in water nanodroplets has recently emerged as a new frontier to investigate structural changes of biomolecules, with perspective applications in ultra-fast drug delivery. We report on the molecular dynamics of the digestive protein Pepsin subjected to a double confinement. The double confinement stemmed from embedding the protein inside a water nanodroplet, which in turn was caged in a nanochannel mimicking the mesoporous silica SBA-15. The nano-bio-droplet, whose size fits with the pore diameter, behaved differently depending on the protonation state of the pore surface silanols. Neutral channel sections allowed for the droplet to flow, while deprotonated sections acted as anchoring piers for the droplet. Inside the droplet, the protein, not directly bonded to the surface, showed a behavior similar to that reported for bulk water solutions, indicating that double confinement should not alter its catalytic activity. Our results suggest that nanobiodroplets, recently fabricated in volatile environments, can be encapsulated and stored in mesoporous silicas.

Translocation intermediates of ubiquitin through an α-hemolysin nanopore: implications for detection of post-translational modifications

Nanopore based sensors constitute a promising approach to single molecule protein characterization being able, in principle, to detect sequences, structural elements and folding states of proteins and polypeptide chains. In narrow nanopores, one of the open issues concerns the coupling between unfolding and translocation. Here, we studied the ubiquitin translocation in an α-hemolysin nanopore, the most widely used pore for nanopore sensing, via all-atom molecular dynamics simulations. We completely characterize the co-translocational unfolding pathway finding that robust translocation intermediates are associated with the rearrangement of secondary structural elements, as also confirmed by coarse grained simulations. An interesting recurrent pattern is the clogging of the α-hemolysin constriction by an N-terminal β-hairpin. This region of ubiquitin is the target of several post-translational modifications. We propose a strategy to detect post-translational modifications at the N-terminal using the α-hemolysin nanopore based on the comparison of the co-translocational unfolding signals associated with modified and unmodified proteins.

Automated Sequence Design of 3D Polyhedral Wireframe DNA Origami with Honeycomb Edges

3D polyhedral wireframe DNA nanoparticles (DNA-NPs) fabricated using scaffolded DNA origami offer complete and independent control over NP size, structure, and asymmetric functionalization on the 10-100 nm scale. However, the complex DNA sequence design needed for the synthesis of these versatile DNA-NPs has limited their widespread use to date. While the automated sequence design algorithms DAEDALUS and vHelix-BSCOR apply to DNA-NPs synthesized using either uniformly dual or hybrid single-dual duplex edges, respectively, these DNA-NPs are relatively compliant mechanically and are therefore of limited utility for some applications. Further, these algorithms are incapable of handling DNA-NP edge designs composed of more than two duplexes, which are needed to enhance DNA-NP mechanical stiffness. As an alternative, here we introduce the scaffolded DNA origami sequence design algorithm TALOS, which is a generalized procedure for the fully automated design of wireframe 3D polyhedra composed of edges of any cross section with an even number of duplexes, and apply it to DNA-NPs composed uniformly of single honeycomb edges. We also introduce a multiway vertex design that enables the fabrication of DNA-NPs with arbitrary edge lengths and vertex angles and apply it to synthesize a highly asymmetric origami object. Sequence designs are demonstrated to fold robustly into target DNA-NP shapes with high folding efficiency and structural fidelity that is verified using single particle cryo-electron microscopy and 3D reconstruction. In order to test its generality, we apply TALOS to design an  in silico library of over 200 DNA-NPs of distinct symmetries and sizes, and for broad impact, we also provide the software as open source for the generation of custom NP designs.

Vorinostat, a possible alternative to metronidazole for the treatment of amebiasis caused by Entamoeba histolytica

AbbreviationsSAHAsuberoylanilide hydroxamic acidEhHDACHistone Deacetylase from Entamoeba histolyticaRgRadius of gyrationRMSDroot-mean-square deviationRMSFroot-mean-square fluctuationMDSmolecular dynamics simulationVMDVisual Molecular DynamicsNAMDNanoscale Molecular DynamicsPBCperiodic boundary conditionsPMEParticle Mesh Ewald3Dthree-dimensionalCαalpha carbonFDAFood and Drug AdministrationnsnanosecondsGPU CUDAGraphics Processing Unit Compute Unified Device ArchitectureCommunicated by Ramaswamy H. Sarma.

Evaluation of the Antimicrobial Activity of Cationic Peptides Loaded in Surface-Modified Nanoliposomes against Foodborne Bacteria

Bacteria are a common group of foodborne pathogens presenting public health issues with a large economic burden for the food industry. Our work focused on a solution to this problem by evaluating antibiotic activity against two bacteria (Listeria monocytogenes and Escherichia coli) of relevance in the field of foodstuffs. We used two approaches: (i) structural modification of the antimicrobial peptides and (ii) nano-vehiculisation of the modified peptides into polymer-coated liposomes. To achieve this, two antimicrobial peptides, herein named ‘peptide +2′ and ‘peptide +5′ were synthesised using the solid phase method. The physicochemical characterisation of the peptides was carried out using measurements of surface tension and dynamic light scattering. Additionally, nanoliposomes were elaborated by the ethanol injection method and coated with a cationic polymer (Eudragit E-100) through the layer-by-layer process. Liposome characterisation, in terms of size, polydispersity and zeta potential, was undertaken using dynamic light scattering. The results show that the degree of hydrophilic modification in the peptide leads to different characteristics of amphipathicity and subsequently to different physicochemical behaviour. On the other hand, antibacterial activity against both bacteria was slightly altered after modifying peptide sequence. Nonetheless, after the encapsulation of the peptides into polymer-coated nano-liposomes, the antibacterial activity increased approximately 2000-fold against that of L. monocytogenes.

Pharmacokinetics and tissue distribution of N-(2-hydroxyphenyl)-2-propylpentanamide in Wistar Rats and its binding properties to human serum albumin

N-(2-hydroxyphenyl)-2-propylpentanamide (HO-AAVPA) is a novel valproic acid derivative that has shown anti-proliferative activity against epitheloid cervix carcinoma (HeLa), rhabdomyosarcoma (A204), and several breast cancer cell lines. The aim of this research was to evaluate the pharmacokinetic profile and tissue distribution of HO-AAVPA in Wistar rats, as well as its human serum albumin binding potential by experimental and in silico methods. A single dose of HO-AAVPA was given to male rats by intravenous, intragastric or intraperitoneal routes at doses of 25, 100, and 100 mg/kg, respectively. Then, blood samples were drawn at predetermined intervals of time, and the HO-AAVPA concentration in the plasma was quantified with a validated HPLC method. The elimination half-life (t_1/2) was approximately 222 min, and the systemic clearance (CL) and apparent volume of distribution (Vd) were 2.20 mL/min/kg and 0.70 L/kg, respectively. The absolute oral bioavailability of HO-AAVPA was 33.8%, and the binding rate of HO-AAVPA with rat plasma proteins was between 66.2% and 83.0%. Additionally, in silico, UV and Raman spectroscopy data showed weak interactions between the test compound and human serum albumin. Thus, the results that were obtained demonstrated that despite its low oral bioavailability, the potential anticancer agent HO-AAVPA exhibits acceptable pharmacokinetic properties that would allow it to reach its site of action and exert its pharmacological effect in Wistar Rats, and it has a convenient profile for future assays to evaluate its human applications.

Cold-induced aldimine bond cleavage by Tris in Bacillus subtilis alanine racemase

The commonly used Tris buffer acts as a surrogate substrate and deformylates pyridoxal phosphate in a bacterial alanine racemase at subzero temperatures.

Accelerated Construction of Kinetic Network Model of Biomolecules Using Steered Molecular Dynamics

A new class of rare event acceleration techniques based on steered molecular dynamics (SMD) simulations is introduced. A stretching force applied on a biomolecule causes it to access large end-to-end distances. Under these conditions the biomolecule undergoes rapid conformational changes that are rare at zero-force conditions. A theory describing kinetics of a biomolecule at various stretching forces is presented. Using the theory, a master-Markov state model (master-MSM) is constructed from rates frequently accessed over a small range of force conditions. The master-MSM is shown to be applicable over a wide range of force conditions. We demonstrate application of the theory to three different biomolecular systems, namely, deca-alanine, TBA (thrombin binding aptamer), and a RNA hairpin. The master-MSM is used to estimate the kinetics at zero-force conditions, i.e., on the unbiased free-energy landscape, resulting inasmuch as 2-6 orders-of-magnitude speed-up over standard molecular dynamics.

Structural Articulation of Biochemical Reactions Using Restrained Geometries and Topology Switching

A strategy named "restrained geometries and topology switching" (RGATS) is presented to obtain detailed trajectories for complex biochemical reactions using molecular mechanics (MM) methods. It enables prediction of realistic dynamical pathways for chemical reactions, especially for accurately characterizing the structural adjustments of highly complex environments to any proximal biochemical reaction. It can be used to generate reactive conformations, model stepwise or concerted reactions in complex environments, and probe the influence of changes in the environment. Its ability to take reactively nonoptimal conformations and generate favorable starting conformations for a biochemical reaction is illustrated for a proton transfer between two model compounds. Its ability to study concerted reactions in explicit solvent is illustrated using proton transfers between an ammonium ion and two conserved histidines in an ammonia transporter channel embedded in a lipid membrane. Its ability to characterize the changes induced by subtle differences in the active site environment is illustrated using nucleotide addition by a DNA polymerase in the presence of two versus three Mg²⁺ ions. RGATS can be employed within any MM program and requires no additional software implementation. This allows the full assortment of computational methods implemented in all available MM programs to be used to tackle virtually any question about biochemical reactions that is answerable without using a quantum mechanical (QM) model. It can also be applied to generate reasonable starting structures for more detailed and expensive QM or QM/MM methods. In particular, this strategy enables rapid prediction of reactant, intermediary, or product state structures in any macromolecular context, with the only requirement being that the structure in any one of these states is either known or can be accurately modeled.