A protocol for preparing explicitly solvated systems for stable molecular dynamics simulations

modXNA: A Modular Approach to Parametrization of Modified Nucleic Acids for Use with Amber Force Fields

Modified nucleic acids have surged as a popular therapeutic route, emphasizing the importance of nucleic acid research in drug discovery and development. Beyond well-known RNA vaccines, antisense oligonucleotides and aptamers can incorporate various modified nucleic acids to target specific biomolecules for various therapeutic activities. Molecular dynamics simulations can accelerate the design and development of these systems with noncanonical nucleic acids by observing intricate dynamic properties and relative stability on the all-atom level. However, modeling these modified systems is challenging due to the time and resources required to parametrize components outside default force field parameters. Here, we present modXNA, a tool to derive and build modified nucleotides for use with Amber force fields. Several nucleic acid systems varying in size and number of modification sites were used to evaluate the accuracy of modXNA parameters, and results indicate the dynamics and structure are preserved throughout the simulations. We detail the protocol for quantum mechanics charge derivation and describe a workflow for implementing modXNA in Amber molecular dynamics simulations, which includes updates and added features to CPPTRAJ.

Fluctuations in Medium Viscosity May Affect the Stability of the CAG Tract in the ATXN2 Gene

Background: Trinucleotide repeats are the cause of many neurodegenerative diseases that are currently incurable. In this regard, the question of the causes of occurrence and methods of prevention or treatment of diseases caused by the expansion of repeats in the CAG tract of the ATXN2 gene remains relevant. Previously, it was shown that the frequency of occurrence of additional OS (open states) zones increases with increasing length of the CAG tract, and the value inverse to the frequency correlates with the age of disease onset. Methods: In this work, the influence of the viscosity of the medium and the external torque on the stability of the CAG tract in the ATXN2 gene was studied using mathematical modeling methods. Results: It has been established that the probability of the appearance of additional OS zones of significant size increases with an increase in the CAG of the tract (k > 40 CAG repeats) for all viscosity values, however, at k ≤ 40, the change in viscosity does not significantly affect the probability of additional OS zones in the tract. Conclusions: It was found that under normal conditions (absence of pathology), viscosity does not have a reliable effect on the stability of the DNA molecule, but when pathology appears, an increase in viscosity contributes to an increase in DNA stability, and, accordingly, a decrease has a negative effect on the stabilization of the DNA molecule. In the zone of close to incomplete penetrance of the disease, viscosity does not have a reliable effect on the stability of the CAG tract.

Anisotropic coarse-grain Monte Carlo simulations of lysozyme, lactoferrin, and NISTmAb by precomputing atomistic models

We develop a multiscale coarse-grain model of the NIST Monoclonal Antibody Reference Material 8671 (NISTmAb) to enable systematic computational investigations of high-concentration physical instabilities such as phase separation, clustering, and aggregation. Our multiscale coarse-graining strategy captures atomic-resolution interactions with a computational approach that is orders of magnitude more efficient than atomistic models, assuming the biomolecule can be decomposed into one or more rigid bodies with known, fixed structures. This method reduces interactions between tens of thousands of atoms to a single anisotropic interaction site. The anisotropic interaction between unique pairs of rigid bodies is precomputed over a discrete set of relative orientations and stored, allowing interactions between arbitrarily oriented rigid bodies to be interpolated from the precomputed table during coarse-grained Monte Carlo simulations. We present this approach for lysozyme and lactoferrin as a single rigid body and for the NISTmAb as three rigid bodies bound by a flexible hinge with an implicit solvent model. This coarse-graining strategy predicts experimentally measured radius of gyration and second osmotic virial coefficient data, enabling routine Monte Carlo simulation of medically relevant concentrations of interacting proteins while retaining atomistic detail. All methodologies used in this work are available in the open-source software Free Energy and Advanced Sampling Simulation Toolkit.

Screening of Potential Inhibitors Targeting the Main Protease Structure of SARS-CoV-2 via Molecular Docking, and Approach with Molecular Dynamics, RMSD, RMSF, H-Bond, SASA and MMGBSA

Severe Acute Respiratory Syndrome caused by a coronavirus is a recent viral infection. There is no scientific evidence or clinical trials to indicate that possible therapies have demonstrated results in suspected or confirmed patients. This work aims to perform a virtual screening of 1430 ligands through molecular docking and to evaluate the possible inhibitory capacity of these drugs about the Mpro protease of Covid-19. The selected drugs were registered with the FDA and available in the virtual drug library, widely used by the population. The simulation was performed using the MolAiCalD algorithm, with a Lamarckian genetic model (GA) combined with energy estimation based on rigid and flexible conformation grids. In addition, molecular dynamics studies were also performed to verify the stability of the receptor-ligand complexes formed through analyses of RMSD, RMSF, H-Bond, SASA, and MMGBSA. Compared to the binding energy of the synthetic redocking coupling (-6.8 kcal/mol/RMSD of 1.34 Å), which was considerably higher, it was then decided to analyze the parameters of only three ligands: ergotamine (-9.9 kcal/mol/RMSD of 2.0 Å), dihydroergotamine (-9.8 kcal/mol/RMSD of 1.46 Å) and olysio (-9.5 kcal/mol/RMSD of 1.5 Å). It can be stated that ergotamine showed the best interactions with the Mpro protease of Covid-19 in the in silico study, showing itself as a promising candidate for treating Covid-19.

Molecular Modeling of Single- and Double-Hydrocarbon-Stapled Coiled-Coil Inhibitors against Bcr-Abl: Toward a Treatment Strategy for CML

The chimeric oncoprotein Bcr-Abl is the causative agent of virtually all chronic myeloid leukemias and a subset of acute lymphoblastic leukemias. As a result of the so-called Philadelphia chromosome translocation t(9;22), Bcr-Abl manifests as a constitutively active tyrosine kinase, which promotes leukemogenesis by activation of cell cycle signaling pathways. Constitutive and oncogenic activation is mediated by an N-terminal coiled-coil oligomerization domain in Bcr (Bcr-CC), presenting a therapeutic target for inhibition of Bcr-Abl activity toward the treatment of Bcr-Abl+ leukemias. Previously, we demonstrated that a rationally designed Bcr-CC mutant, CCmut3, exerts a dominant negative effect upon Bcr-Abl activity by preferential oligomerization with Bcr-CC. Moreover, we have shown that conjugation to a leukemia-specific cell-penetrating peptide (CPP-CCmut3) improves intracellular delivery and activity. However, our full-length CPP-CCmut3 construct (81 aa) is encumbered by an intrinsically high degree of conformational variability and susceptibility to proteolytic degradation relative to traditional small-molecule therapeutics. Here, we iterate a new generation of CCmut3 inhibitors against Bcr-CC-mediated Bcr-Abl assembly designed to address these constraints through incorporation of all-hydrocarbon staples spanning i and i + 7 positions in α-helix 2 (CPP-CCmut3-st). We utilize computational modeling and biomolecular simulation to evaluate single- and double-stapled CCmut3 candidates in silico for dynamics and binding energetics. We further model a truncated system characterized by the deletion of α-helix 1 and the flexible loop linker, which are known to impart high conformational variability. To study the impact of the N-terminal cyclic CPP toward model stability and inhibitor activity, we also model the full-length and truncated systems devoid of the CPP, with a cyclized CPP, and with an open-configuration CPP, for a total of six systems that comprise our library. From this library, we present lead-stapled peptide candidates to be synthesized and evaluated experimentally as our next iteration of inhibitors against Bcr-Abl.

De novo design of proteins housing excitonically coupled chlorophyll special pairs

Natural photosystems couple light harvesting to charge separation using a 'special pair' of chlorophyll molecules that accepts excitation energy from the antenna and initiates an electron-transfer cascade. To investigate the photophysics of special pairs independently of the complexities of native photosynthetic proteins, and as a first step toward creating synthetic photosystems for new energy conversion technologies, we designed C2-symmetric proteins that hold two chlorophyll molecules in closely juxtaposed arrangements. X-ray crystallography confirmed that one designed protein binds two chlorophylls in the same orientation as native special pairs, whereas a second designed protein positions them in a previously unseen geometry. Spectroscopy revealed that the chlorophylls are excitonically coupled, and fluorescence lifetime imaging demonstrated energy transfer. The cryo-electron microscopy structure of a designed 24-chlorophyll octahedral nanocage with a special pair on each edge closely matched the design model. The results suggest that the de novo design of artificial photosynthetic systems is within reach of current computational methods.

A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta

Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-β (Aβ42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the β-content of Aβ42; however, it increases the tendency of the monomeric Aβ42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aβ42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aβ42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.

PBP-A, a cyanobacterial DD-peptidase with high specificity for amidated muropeptides, exhibits pH-dependent promiscuous activity harmful to Escherichia coli

Penicillin binding proteins (PBPs) are involved in biosynthesis, remodeling and recycling of peptidoglycan (PG) in bacteria. PBP-A from Thermosynechococcus elongatus belongs to a cyanobacterial family of enzymes sharing close structural and phylogenetic proximity to class A β-lactamases. With the long-term aim of converting PBP-A into a β-lactamase by directed evolution, we simulated what may happen when an organism like Escherichia coli acquires such a new PBP and observed growth defect associated with the enzyme activity. To further explore the molecular origins of this harmful effect, we decided to characterize deeper the activity of PBP-A both in vitro and in vivo. We found that PBP-A is an enzyme endowed with DD-carboxypeptidase and DD-endopeptidase activities, featuring high specificity towards muropeptides amidated on the D-iso-glutamyl residue. We also show that a low promiscuous activity on non-amidated peptidoglycan deteriorates E. coli's envelope, which is much higher under acidic conditions where substrate discrimination is mitigated. Besides expanding our knowledge of the biochemical activity of PBP-A, this work also highlights that promiscuity may depend on environmental conditions and how it may hinder rather than promote enzyme evolution in nature or in the laboratory.

On the Cholesterol Raising Effect of Coffee Diterpenes Cafestol and 16-O-Methylcafestol: Interaction with Farnesoid X Receptor

The diterpene cafestol represents the most potent cholesterol-elevating compound known in the human diet, being responsible for more than 80% of the effect of coffee on serum lipids, with a mechanism still not fully clarified. In the present study, the interaction of cafestol and 16-O-methylcafestol with the stabilized ligand-binding domain (LBD) of the Farnesoid X Receptor was evaluated by fluorescence and circular dichroism. Fluorescence quenching was observed with both cafestol and 16-O-methylcafestol due to an interaction occurring in the close environment of the tryptophan W454 residue of the protein, as confirmed by docking and molecular dynamics. A conformational change of the protein was also observed by circular dichroism, particularly for cafestol. These results provide evidence at the molecular level of the interactions of FXR with the coffee diterpenes, confirming that cafestol can act as an agonist of FXR, causing an enhancement of the cholesterol level in blood serum.

MPI-parallelization of the grid inhomogeneous solvation theory calculation

The grid inhomogeneous solvation theory (GIST) method requires the often time-consuming calculation of water-water and water-solute energy on a grid. Previous efforts to speed up this calculation include using OpenMP, GPUs, and particle mesh Ewald. This article details how the speed of this calculation can be increased by parallelizing it with MPI, where trajectory frames are divided among multiple processors. This requires very little communication between individual processes during trajectory processing, meaning the calculation scales well to large processor counts. This article also details how the entropy calculation, which must happen after trajectory processing since it requires information from all trajectory frames, is parallelized via MPI. This parallelized GIST method has been implemented in the freely-available CPPTRAJ analysis software.

Enhancing biocatalyst performance through immobilization of lipase (Eversa® Transform 2.0) on hybrid amine-epoxy core-shell magnetic nanoparticles

Magnetic nanoparticles were functionalized with polyethylenimine (PEI) and activated with epoxy. This support was used to immobilize Lipase (Eversa® Transform 2.0) (EVS), optimization using the Taguchi method. XRF, SEM, TEM, XRD, FTIR, TGA, and VSM performed the characterizations. The optimal conditions were immobilization yield (I.Y.) of 95.04 ± 0.79 %, time of 15 h, ionic load of 95 mM, protein load of 5 mg/g, and temperature of 25 °C. The maximum loading capacity was 25 mg/g, and its stability in 60 days of storage showed a negligible loss of only 9.53 % of its activity. The biocatalyst demonstrated better stability at varying temperatures than free EVS, maintaining 28 % of its activity at 70 °C. It was feasible to esterify free fatty acids (FFA) from babassu oil with the best reaction of 97.91 % and ten cycles having an efficiency above 50 %. The esterification of produced biolubricant was confirmed by NMR, and it displayed kinematic viscosity and density of 6.052 mm2/s and 0.832 g/cm3, respectively, at 40 °C. The in-silico study showed a binding affinity of -5.8 kcal/mol between EVS and oleic acid, suggesting a stable substrate-lipase combination suitable for esterification.

An antibacterial packaging film based on amylose starch with quaternary ammonium salt chitosan and its application for meat preservation

A new generation of food packaging films is gradually replacing traditional plastic packaging films because of their biodegradability, safety, and some functional properties such as anti-bacterial and oxidant resistance. In the present work, an antibacterial packing film based on amylose starch and 2-hydroxypropyl-trimethylammonium chloride chitosan (HTCC) was prepared for meat preservation. The interfacial bonding mechanism between amylose, HTCC, and glutaraldehyde (GA) was determined experimentally and through molecular dynamics (MD) simulation. The macromolecular chains of amylose starch and HTCC became entangled via inter-molecular H-bonds and then cross-linked with GA via the Schiff base reaction. The interaction of amylose starch and HTCC improved the mechanical properties of the amylose films. Compared with the amylose films, the tensile strength and elongation at break of the optimal HTCC/amylose films reached to 16.13 MPa (an increase of 206.65 %) and 53.86 % (an increase of 109.49 %). The HTCC/amylose films were found to provide obvious bacteriostatic performance, a relatively low cytotoxicity, the lower transmittance in the UV region, and thus the ability to enhance the preservation of fresh meat. These excellent characteristics therefore suggest that HTCC/amylose films might be promising candidates for application in antibacterial food packaging films.

Synthesis, structural, characterization, antibacterial and antibiotic modifying activity, ADMET study, molecular docking and dynamics of chalcone (E)-1-(4-aminophenyl)-3-(4-nitrophenyl)prop-2-en-1-one in strains of Staphylococcus aureus carrying NorA and MepA efflux pumps

Chalcones have an open chain flavonoid structure that can be obtained from natural sources or by synthesis and are widely distributed in fruits, vegetables, and tea. They have a simple and easy to handle structure due to the α-β-unsaturated bridge responsible for most biological activities. The facility to synthesize chalcones combined with its efficient in combating serious bacterial infections make these compounds important agents in the fight against microorganisms. In this work, the chalcone (E)-1-(4-aminophenyl)-3-(4-nitrophenyl)prop-2-en-1-one (HDZPNB) was characterized by spectroscopy and electronic methods. In addition, microbiological tests were performed to investigate the modulator potential and efflux pump inhibition on S. aureus multi-resistant strains. The modulating effect of HDZPNB chalcone in association with the antibiotic norfloxacin, on the resistance of the S. aureus 1199 strain, resulted in increase the MIC. In addition, when HDZPNB was associated with ethidium bromide (EB), it caused an increase in the MIC value, thus not inhibiting the efflux pump. For the strain of S. aureus 1199B, carrying the NorA pump, the HDZPNB associated with norfloxacin showed no modulatory, and when the chalcone was used in association with EB, it had no inhibitory effect on the efflux pump. For the tested strain of S. aureus K2068, which carries the MepA pump, it can be observed that the chalcone together the antibiotic resulted in an increase the MIC. On the other hand, when chalcone was used in association with EB, it caused a decrease in bromide MIC, equal to the reduction caused by standard inhibitors. Thus, these results indicate that the HDZPNB could also act as an inhibitor of the S. aureus gene overexpressing pump MepA. The molecular docking reveals that chalcone has a good binding energies -7.9 for HDZPNB/MepA complexes, molecular dynamics simulations showed that Chalcone/MetA complexes showed good stability of the structure in an aqueous solution, and ADMET study showed that the chalcone has a good oral bioavailability, high passive permeability, low risk of efflux, low clearance rate and low toxic risk by ingestion. The microbiological tests show that the chalcone can be used as a possible inhibitor of the Mep A efflux pump.Communicated by Ramaswamy H. Sarma.

Anxiolytic effects of N-(4,5-dihydro-5-oxo-1,2-dithiolo-[4,3,b]-pyrrole-6-yl)-N-methylformamide, a pyrroloformamide isolated from a marine Streptomyces sp., in adult zebrafish by the 5-HT system

General anxiety disorders are among the most prevalent mental health problems worldwide. The emergence and development of anxiety disorders can be due to genetic (30-50%) or non-genetic (50-70%) factors. Despite medical progress, available pharmacotherapies are sometimes ineffective or can cause undesirable side effects. Thus, it becomes necessary to discover new safe and effective drugs against anxiety. This study evaluated the anxiolytic effect in adult zebrafish (Danio rerio) of a natural pyrroloformamide (PFD), N-(4,5-dihydro-5-oxo-1,2-dithiolo-[4,3,b]-pyrrole-6-yl)-N-methylformamide, isolated from a Streptomyces sp. bacterium strain recovered from the ascidian Eudistoma vannamei. The complete structure of PFD was determined by a detailed NMR analysis, including 1H-13C and 1H-15N-HBMC data. In addition, conformational and DFT computational studies also were performed. A group of fishes (n = 6) was treated orally with PFD (0.1, 0.5 and 1.0 mg/mL; 20 μL) and subjected to locomotor activity and light/dark tests, as well as, acute toxicity 96 h. The involvement of the GABAergic and serotonergic (5-HT) systems was investigated using flumazenil (a silent modulator of GABA receptor) and 5-HT1, 5-HT2A/2C and 5-HTR3A/3B receptors antagonists, known as pizotifen, granisetron and cyproheptadine, respectively. PFD was nontoxic, reduced locomotor activity and promoted the anxiolytic effect in zebrafish. Flumazenil did not inhibit the anxiolytic effect of the PFD via the GABAergic system. This effect was reduced by a pretreatment with pizotifen and granisetron, and was not reversed after treatment with cyproheptadine. Molecular docking and dynamics studies confirmed the interaction of PFD with the 5-HT receptor.Communicated by Ramaswamy H. Sarma.

Computational Modeling of Stapled Coiled-Coil Inhibitors Against Bcr-Abl: Toward a Treatment Strategy for CML

The chimeric oncoprotein Bcr-Abl is the causative agent of virtually all chronic myeloid leukemias (CML) and a subset of acute lymphoblastic leukemias (ALL). As a result of the so-called Philadelphia Chromosome translocation t(9;22), Bcr-Abl manifests as a constitutively active tyrosine kinase which promotes leukemogenesis by activation of cell cycle signaling pathways. Constitutive and oncogenic activation is mediated by an N-terminal coiled-coil oligomerization domain in Bcr (Bcr-CC), presenting a therapeutic target for inhibition of Bcr-Abl activity toward the treatment of Bcr-Abl+ leukemias. Previously, we demonstrated that a rationally designed Bcr-CC mutant, CCmut3, exerts a dominant negative effect upon Bcr-Abl activity by preferential oligomerization with Bcr-CC. Moreover, we have shown conjugation to a leukemia-specific cell-penetrating peptide (CPP-CCmut3) improves intracellular delivery and activity. However, our full-length CPP-CCmut3 construct (81 aa) is encumbered by an intrinsically high degree of conformational variability and susceptibility to proteolytic degradation, relative to traditional small molecule therapeutics. Here, we iterate a new generation of our inhibitor against Bcr-CC mediated Bcr-Abl assembly that is designed to address these constraints through incorporation of all-hydrocarbon staples spanning i, i + 7 positions in helix α2 (CPP-CCmut3-st). We utilize computational modeling and biomolecular simulation to design and characterize single and double staple candidates in silico, evaluating binding energetics and building upon our seminal work modeling single hydrocarbon staples when applied to a truncated Bcr-CC sequence. This strategy enables us to efficiently build, characterize, and screen lead single/double stapled CPP-CCmut3-st candidates for experimental studies and validation in vitro and in vivo. In addition to full-length CPP-CCmut, we model a truncated system characterized by deletion of helix α1 and the flexible-loop linker, which are known to impart high conformational variability. To study the impact of the N-terminal cyclic CPP toward model stability and inhibitor activity, we also model the full-length and truncated systems without CPP, with cyclized CPP, and with linear CPP, for a total of six systems which comprise our library. From this library, we present lead stapled peptide candidates to be synthesized and evaluated experimentally as our next-generation inhibitors against Bcr-Abl.

Identification of Dual-Target Inhibitors for Epidermal Growth Factor Receptor and AKT: Virtual Screening Based on Structure and Molecular Dynamics Study

Epidermal growth factor EGFR is an important target for non-small cell lung (NSCL) cancer, and inhibitors of the AKT protein have been used in many cancer treatments, including those for NSCL cancer. Therefore, searching small molecular inhibitors which can target both EGFR and AKT may help cancer treatment. In this study, we applied a ligand-based pharmacophore model, molecular docking, and MD simulation methods to search for potential inhibitors of EGFR and then studied dual-target inhibitors of EGFR and AKT by screening the immune-oncology Chinese medicine (TCMIO) database and the human endogenous database (HMDB). It was found that TCMIO89212, TCMIO90156, and TCMIO98874 had large binding free energies with EGFR and AKT, and HMDB0012243 also has the ability to bind to EGFR and AKT. These results may provide valuable information for further experimental study.

Divalent ions as mediators of carbonylation in cardiac myosin binding protein C

The dosing and efficacy of chemotherapeutic drugs can be limited by toxicity caused by off-pathway reactions. One hypothesis for how such toxicity arises is via metal-catalyzed oxidative damage of cardiac myosin binding protein C (cMyBP-C) found in cardiac tissue. Previous research indicates that metal ion mediated reactive oxygen species induce high levels of protein carbonylation, changing the structure and function of this protein. In this work, we use long timescale all-atom molecular dynamics simulations to investigate the ion environment surrounding the C0 and C1 subunits of cMyBP-C responsible for actin binding. We show that divalent cations are co-localized with protein carbonylation-prone amino acid residues and that carbonylation of these residues can lead to site-specific interruption to the actin-cMyBP-C binding.

Effects of glycans and hinge on dynamics in the IgG1 Fc

The crystallizable fragment (Fc) domain of immunoglobulin subclass IgG1 antibodies is engineered for a wide variety of pharmaceutical applications. Two important structural variables in Fc constructs are the hinge region connecting the Fc to the antigen binding fragments (Fab) and the glycans present in various glycoforms. These components affect receptor binding interactions that mediate immune activation. To design new antibody drugs, a robust in silico method for linking stability to structural changes is necessary. In this work, all-atom simulations were used to compare the dynamic behavior of the four structural variants arising from presence or absence of the hinge and glycans. We expressed the simplest of these constructs, the 'minimal Fc' with no hinge and no glycans, in Escherichia coli and report its crystal structure. The 'maximal Fc' that includes full hinge and G0F/G1F glycans is based on a previously reported structure, Protein Data Bank (PDB) ID: 5VGP. These, along with two intermediate structures (with only the glycans or with only the hinge) were used to independently measure the stability effects of the two structural variables using umbrella sampling simulations. Principal component analysis (PCA) was used to determine free energy effects along the Fc's dominant mode of motion. This work provides a comprehensive picture of the effects of hinge and glycans on Fc dynamics and stability.Communicated by Ramaswamy H. Sarma.

Computational completion of the Aurora interaction region of N-Myc in the Aurora a kinase complex

Inhibiting protein-protein interactions of the Myc family is a viable pharmacological strategy for modulation of the levels of Myc oncoproteins in cancer. Aurora A kinase (AurA) and N-Myc interaction is one of the most attractive targets of this strategy because formation of this complex blocks proteasomal degradation of N-Myc in neuroblastoma. Two crystallization studies have captured this complex (PDB IDs: 5g1x, 7ztl), partially resolving the AurA interaction region (AIR) of N-Myc. Prompted by the missing N-Myc fragment in these crystal structures, we modeled the complete structure between AurA and N-Myc, and comprehensively analyzed how the incomplete and complete N-Myc behave in complex by molecular dynamics simulations. Molecular dynamics simulations of the incomplete PDB complex (5g1x) repeatedly showed partial dissociation of the short N-Myc fragment (61-89) from the kinase. The missing N-Myc (19-60) fragment was modeled utilizing the N-terminal lobe of AurA as the protein-protein interaction surface, wherein TPX2, a well-known partner of AurA, also binds. Binding free energy calculations along with flexibility analysis confirmed that the complete AIR of N-Myc stabilizes the complex, accentuating the N-terminal lobe of AurA as a binding site for the missing N-Myc fragment (19-60). We further generated additional models consisting of only the missing N-Myc (19-60), and the fused form of TPX2 (7-43) and N-Myc (61-89). These partners also formed more stable interactions with the N-terminal lobe of AurA than did the incomplete N-Myc fragment (61-89) in the 5g1x complex. Altogether, this study provides structural insights into the involvement of the N-terminus of the AIR of N-Myc and the N-terminal lobe of AurA in formation of a stable complex, reflecting its potential for effective targeting of N-Myc.

Local structural flexibility drives oligomorphism in computationally designed protein assemblies

Many naturally occurring protein assemblies have dynamic structures that allow them to perform specialized functions. For example, clathrin coats adopt a wide variety of architectures to adapt to vesicular cargos of various sizes. Although computational methods for designing novel self-assembling proteins have advanced substantially over the past decade, most existing methods focus on designing static structures with high accuracy. Here we characterize the structures of three distinct computationally designed protein assemblies that each form multiple unanticipated architectures, and identify flexibility in specific regions of the subunits of each assembly as the source of structural diversity. Cryo-EM single-particle reconstructions and native mass spectrometry showed that only two distinct architectures were observed in two of the three cases, while we obtained six cryo-EM reconstructions that likely represent a subset of the architectures present in solution in the third case. Structural modeling and molecular dynamics simulations indicated that the surprising observation of a defined range of architectures, instead of non-specific aggregation, can be explained by constrained flexibility within the building blocks. Our results suggest that deliberate use of structural flexibility as a design principle will allow exploration of previously inaccessible structural and functional space in designed protein assemblies.

Structure based virtual screening and molecular simulation study of FDA-approved drugs to inhibit human HDAC6 and VISTA as dual cancer immunotherapy

Cancer immunotherapy has significantly contributed to the treatment of various types of cancers mainly by targeting immune checkpoint inhibitors (ICI). Among them, V-domain immunoglobulin suppressor of T cell activation (VISTA) has been explored as a promising therapeutic target. Besides, histone deacetylase 6 (HDAC6) has been demonstrated to be efficacious target for several cancers. The current theoretical work was performed to explore the virtual repurposing of the FDA-approved drugs as inhibitors against these two (VISTA and HDAC6) cancers therapeutic targets. The crystal structure of the two proteins were downloaded from PDB and subjected to virtual screening by DrugRep webserver while using FDA-approved drugs library as ligands database. Our study revealed that Oxymorphone and Bexarotene are the top-ranked inhibitors of VISTA and HDAC6, respectively. The docking score of Bexarotene was predicted as - 10 kcal/mol while the docking score of Oxymorphone was predicted as - 6.2 kcal/mol. Furthermore, a total of 100 ns MD simulation revealed that the two drugs Oxymorphone and Bexarotene formed stable complexes with VISTA and HDAC6 drug targets. As compared to the standard drug the two drugs Oxymorphone and Bexarotene revealed great stability during the whole 100 ns MD simulation. The binding free energy calculation further supported the Root Mean Square Deviation (RMSD) result which stated that as compared to the ref/HDAC6 (- 18.0253 ± 2.6218) the binding free energy score of the Bexarotene/HDAC6 was good (- 51.9698 ± 3.1572 kcal/mol). The binding free energy score of Oxymorphone/VISTA and Ref/VISTA were calculated as - 36.8323 ± 3.4565, and - 21.5611 ± 4.8581 respectively. In conclusion, the two drugs deserve further consideration as cancer treatment option.

Designing Efficient Enzymes: Eight Predicted Mutations Convert a Hydroxynitrile Lyase into an Efficient Esterase

Hydroxynitrile lyase from rubber tree (HbHNL) shares 45% identical amino acid residues with the homologous esterase from tobacco, SABP2, but the two enzymes catalyze different reactions. The x-ray structures reveal a serine-histidine-aspartate catalytic triad in both enzymes along with several differing amino acid residues within the active site. Previous exchange of three amino acid residues in the active site of HbHNL with the corresponding amino acid residue in SABP2 (T11G-E79H-K236M) created variant HNL3, which showed low esterase activity toward p-nitrophenyl acetate. Further structure comparison reveals additional differences surrounding the active site. HbHNL contains an improperly positioned oxyanion hole residue and differing solvation of the catalytic aspartate. We hypothesized that correcting these structural differences would impart good esterase activity on the corresponding HbHNL variant. To predict the amino acid substitutions needed to correct the structure, we calculated shortest path maps for both HbHNL and SABP2, which reveal correlated movements of amino acids in the two enzymes. Replacing four amino acid residues (C81L-N104T-V106F-G176S) whose movements are connected to the movements of the catalytic residues yielded variant HNL7TV (stabilizing substitution H103V was also added), which showed an esterase catalytic efficiency comparable to that of SABP2. The x-ray structure of an intermediate variant, HNL6V, showed an altered solvation of the catalytic aspartate and a partially corrected oxyanion hole. This dramatic increase in catalytic efficiency demonstrates the ability of shortest path maps to predict which residues outside the active site contribute to catalytic activity.

De novo design of energy transfer proteins housing excitonically coupled chlorophyll special pairs

Natural photosystems couple light harvesting to charge separation using a "special pair" of chlorophyll molecules that accepts excitation energy from the antenna and initiates an electron-transfer cascade. To investigate the photophysics of special pairs independent of complexities of native photosynthetic proteins, and as a first step towards synthetic photosystems for new energy conversion technologies, we designed C2-symmetric proteins that precisely position chlorophyll dimers. X-ray crystallography shows that one designed protein binds two chlorophylls in a binding orientation matching native special pairs, while a second positions them in a previously unseen geometry. Spectroscopy reveals excitonic coupling, and fluorescence lifetime imaging demonstrates energy transfer. We designed special pair proteins to assemble into 24-chlorophyll octahedral nanocages; the design model and cryo-EM structure are nearly identical. The design accuracy and energy transfer function of these special pair proteins suggest that de novo design of artificial photosynthetic systems is within reach of current computational methods.

AmberMDrun: A Scripting Tool for Running Amber MD in an Easy Way

MD simulations have been widely applied and become a powerful tool in the field of biomacromolecule simulations and computer-aided drug design, etc., which can estimate binding free energy between receptor and ligand. However, the inputs and force field preparation for performing Amber MD is somewhat complicated, and challenging for beginners. To address this issue, we have developed a script for automatically preparing Amber MD input files, balancing the system, performing Amber MD for production, and predicting receptor-ligand binding free energy. This script is open-source, extensible and can support customization. The core code is written in C++ and has a Python interface, providing both efficient performance and convenience.

Isolation, characterization and in silico study of propenamide alkaloids from Hymenoepmecis bicolor poison against active μ-opioid receptor

Some insects produce venoms to defend against predators and directly interact with opioid receptors. In the present study, it was identified two alkaloids in the wasp venom species Hymenoepimecis bicolor. It was demonstrated that these could act as potential inhibitors of opioid receptors through their robust affinity to the receptors. The interaction profile was given to opioid receptors (μOR), with 60% of targets similar to alkaloid 1, with 0.25 probability, and 46.7% of targets similar to alkaloid 2, with a probability 0.17 of affinity as a target, which is considered signaling macromolecules and can mediate the most potent analgesic and addictive properties of opiate alkaloids. Notably, both alkaloids showed -7.6 kcal/mol affinity to the morphine agonies through six residues, Gly124, Asp147, Trp293, Ile296, Ile322, and Tyr326. These observations suggest further research on opioid receptors using in vitro studies of possible therapeutic applications.Communicated by Ramaswamy H. Sarma.

A multipoint guidance mechanism for β-barrel folding on the SAM complex

Mitochondrial β-barrel proteins are essential for the transport of metabolites, ions and proteins. The sorting and assembly machinery (SAM) mediates their folding and membrane insertion. We report the cryo-electron microscopy structure of the yeast SAM complex carrying an early eukaryotic β-barrel folding intermediate. The lateral gate of Sam50 is wide open and pairs with the last β-strand (β-signal) of the substrate-the 19-β-stranded Tom40 precursor-to form a hybrid barrel in the membrane plane. The Tom40 barrel grows and curves, guided by an extended bridge with Sam50. Tom40's first β-segment (β1) penetrates into the nascent barrel, interacting with its inner wall. The Tom40 amino-terminal segment then displaces β1 to promote its pairing with Tom40's last β-strand to complete barrel formation with the assistance of Sam37's dynamic α-protrusion. Our study thus reveals a multipoint guidance mechanism for mitochondrial β-barrel folding.

Contributions from ClpS surface residues in modulating N-terminal peptide binding and their implications for NAAB development

Numerous technologies are currently in development for use in next-generation protein sequencing platforms. A notable published approach employs fluorescently-tagged binding proteins to identity the N-terminus of immobilized peptides, in-between rounds of digestion. This approach makes use of N-terminal amino acid binder (NAAB) proteins, which would identify amino acids by chemical and shape complementarity. One source of NAABs is the ClpS protein family, which serve to recruit proteins to bacterial proteosomes based on the identity of the N-terminal amino acid. In this study, a Thermosynechococcus vestitus (also known as Thermosynechococcus elongatus) ClpS2 protein was used as the starting point for direct evolution of an NAAB with affinity and specificity for N-terminal leucine. Enriched variants were analyzed and shown to improve the interaction between the ClpS surface and the peptide chain, without increasing promiscuity. Interestingly, interactions were found that were unanticipated which favor different charged residues located at position 5 from the N-terminus of a target peptide.

Naphthoquinones biflorin and bis-biflorin (Capraria biflora) as possible inhibitors of the fungus Candida auris polymerase: molecular docking, molecular dynamics, MM/GBSA calculations and in silico drug-likeness study

A new worldwide concern has emerged with the recent emergence of infections caused by Candida auris. This reflects its comparative ease of transmission, substantial mortality, and the increasing level of resistance seen in the three major classes of antifungal drugs. Efforts to create a better design for structure-based drugs that described numerous modifications and the search for secondary metabolic structures derived from plant species are likely to reduce the virulence of several fungal pathogens. In this context, the present work aimed to evaluate in silico two naphthoquinones isolated from the roots of Capraria biflora, biflorin, and its dimmer, bis-biflorin, as potential inhibitors of Candida auris polymerase. Based on the simulation performed with the two naphthoquinones, biflorin and bis-biflorin, it can be stated that bis-biflorin showed the best interactions with Candida auris polymerase. Still, biflorin also demonstrated favorable coupling energy. Predictive pharmacokinetic assays suggest that biflorin has high oral bioavailability and more excellent metabolic stability compared to the bis-biflorin analogue. constituting a promising pharmacological tool.Communicated by Ramaswamy H. Sarma.

Mechanisms of Electron Transfer Rate Modulations in Cytochrome P450 BM3

Bacterial cytochromes P450 BM3 (CYP450 BM3) catalyze reactions of industrial importance. Despite many successful biotransformations, robust (re)design for novel applications remains challenging. Rational design and evolutionary approaches are not always successful, highlighting a lack of complete understanding of the mechanisms of electron transfer (ET) modulations. Thus, the full potential of CYP450 reactions remains under-exploited. In this work, we report the first molecular dynamics (MD)-based explicit prediction of BM3 ET parameters (reorganization energies; λ and ET free energies; ΔG°), and log ET rates (log k_ET) using the Marcus theory. Overall, the calculated ET rates for the BM3 wild-type (WT), mutants (F393 and L86), ligand-bound state, and ion concentrations agree well with experimental data. In ligand-free (LF) BM3, mutations modulate k_ET via ET ΔG°. Simulations show that the experimental ET rate enhancement is due to increased driving force (more negative ΔG°) upon ligation. This increase is related to the protein reorganization required to accommodate the ligand in the binding pocket rather than binding interactions with the ligand. Our methodology (CYPWare 1.0) automates all the stages of the MD simulation step-up, energy calculations, and estimation of ET parameters. CYPWare 1.0 and this work thus represent an important advancement in the CYP450 ET rate predictions, which has the potential to guide the redesign of ET enzymes. This program and a Web tool are available on GitHub for academic research.

Quantifying the effects of lossy compression on energies calculated from molecular dynamics trajectories

Molecular dynamics (MD) simulations are now able to routinely reach timescales of microseconds and beyond. This has led to a corresponding increase in the amount of MD trajectory data that needs to be stored, particularly when those trajectories contain explicit solvent molecules. As such, it is desirable to be able to compress trajectory data while still retaining as much of the original information as possible. In this work, we describe compressing MD trajectory data using the NetCDF4/HDF5 file format, making use of quantization of the original positions to achieve better compression ratios. We also analyze the affect this has on both the resulting positions and the energies calculated from post-processing these trajectories, and recommend an optimal level of quantization. Overall we find the NetCDF4/HDF5 format to be an excellent choice for storing MD trajectory data in terms of speed, compressibility, and versatility.

Efficient fluorescent recognition of ATP/GTP by a water-soluble bisquinolinium pyridine-2,6-dicarboxamide compound. Crystal structures, spectroscopic studies and interaction mode with DNA

The new dicationic pyridine-2,6-dicarboxamide-based compound 1 bearing two N-alkylquinolinium units was synthesized, structurally determined by single-crystal X-ray diffraction, and studied in-depth as a fluorescent receptor for nucleotides and inorganic phosphorylated anions in pure water. The addition of nucleotides to 1 at pH = 7.0 quenches its blue emission with a selective affinity towards adenosine 5'-triphosphate (ATP) and guanosine 5'-tripohosphate (GTP) over other nucleotides such CTP, UTP, ADP, AMP, dicarboxylates and inorganic anions. On the basis of the spectroscopic tools (1H, 31P NMR, UV-vis, fluorescence), MS measurements and DFT calculations, receptor 1 binds ATP with high affinity (log K = 5.04) through the simultaneous formation of strong hydrogen bonds and π-π interactions between the adenosine fragment and quinolinium ring with binding energy calculated in 8.7 kcal mol-1. High affinity for ATP/GTP is attributed to the high acidity of amides and preorganized rigid structure of 1. Receptor 1 is an order of magnitude more selective for ATP than GTP. An efficient photoinduced electron transfer quenching mechanism with simultaneous receptor-ATP complexation in both the excited and ground states is proposed. Additionally, multiple spectroscopic studies and molecular dynamics simulations showed that 1 can intercalate into DNA base pairs.

BioExcel Building Blocks Workflows (BioBB-Wfs), an integrated web-based platform for biomolecular simulations

We present BioExcel Building Blocks Workflows, a web-based graphical user interface (GUI) offering access to a collection of transversal pre-configured biomolecular simulation workflows assembled with the BioExcel Building Blocks library. Available workflows include Molecular Dynamics setup, protein-ligand docking, trajectory analyses and small molecule parameterization. Workflows can be launched in the platform or downloaded to be run in the users’ own premises. Remote launching of long executions to user's available High-Performance computers is possible, only requiring configuration of the appropriate access credentials. The web-based graphical user interface offers a high level of interactivity, with integration with the NGL viewer to visualize and check 3D structures, MDsrv to visualize trajectories, and Plotly to explore 2D plots. The server requires no login but is recommended to store the users’ projects and manage sensitive information such as remote credentials. Private projects can be made public and shared with colleagues with a simple URL. The tool will help biomolecular simulation users with the most common and repetitive processes by means of a very intuitive and interactive graphical user interface. The server is accessible at https://mmb.irbbarcelona.org/biobb-wfs.

Conformational Dynamics of the Hepatitis B Virus Pre-genomic RNA on Multiple Time Scales: Implications for Viral Replication

Human hepatitis B virus (HBV) replication is initiated by the binding of the viral polymerase (P) to epsilon (ε), an ≈85-nucleotide (nt) cis-acting regulatory stem-loop RNA located at the 5'-end of the pre-genomic RNA (pgRNA). This interaction triggers P and pgRNA packaging and protein-primed reverse transcription and is therefore an attractive therapeutic target. Our recent nuclear magnetic resonance (NMR) structure of ε provides a useful starting point toward a detailed understanding of HBV replication, and hints at the functional importance of ε dynamics. Here, we present a detailed description of ε motions on the ps to ns and μs to ms time scales by NMR spin relaxation and relaxation dispersion, respectively. We also carried out molecular dynamics simulations to provide additional insight into ε conformational dynamics. These data outline a series of complex motions on multiple time scales within ε. Moreover, these motions occur in mostly conserved nucleotides from structural regions (i.e., priming loop, pseudo-triloop, and U43 bulge) that biochemical and mutational studies have shown to be essential for P binding, P-pgRNA packaging, protein-priming, and DNA synthesis. Taken together, our work implicates RNA dynamics as an integral feature that governs HBV replication.

Computational design of mechanically coupled axle-rotor protein assemblies

Natural molecular machines contain protein components that undergo motion relative to each other. Designing such mechanically constrained nanoscale protein architectures with internal degrees of freedom is an outstanding challenge for computational protein design. Here we explore the de novo construction of protein machinery from designed axle and rotor components with internal cyclic or dihedral symmetry. We find that the axle-rotor systems assemble in vitro and in vivo as designed. Using cryo-electron microscopy, we find that these systems populate conformationally variable relative orientations reflecting the symmetry of the coupled components and the computationally designed interface energy landscape. These mechanical systems with internal degrees of freedom are a step toward the design of genetically encodable nanomachines.

HDX-MS and MD Simulations Provide Evidence for Stabilization of the IgG1-FcγRIa (CD64a) Immune Complex Through Intermolecular Glycoprotein Bonds

Previous reports present different models for the stabilization of the Fc-FcγRI immune complex. Although accord exists on the importance of L235 in IgG1 and some hydrophobic contacts for complex stabilization, discord exists regarding the existence of stabilizing glycoprotein contacts between glycans of IgG1 and a conserved FG-loop (¹⁷¹MGKHRY¹⁷⁶) of FcγRIa. Complexes formed from the FcγRIa receptor and IgG1s containing biantennary glycans with N-acetylglucosamine, galactose, and α2,6-N-acetylneuraminic terminations were measured by hydrogen-deuterium exchange mass spectrometry (HDX-MS), classified for dissimilarity with Welch's ANOVA and Games-Howell post hoc procedures, and modeled with molecular dynamics (MD) simulations. For each glycoform of the IgG1-FcγRIa complex peptic peptides of Fab, Fc and FcγRIa report distinct H/D exchange rates. MD simulations corroborate the differences in the peptide deuterium content through calculation of the percent of time that transient glycan-peptide bonds exist. These results indicate that stability of IgG1-FcγRIa complexes correlate with the presence of intermolecular glycoprotein interactions between the IgG1 glycans and the ¹⁷³KHR¹⁷⁵ motif within the FG-loop of FcγRIa. The results also indicate that intramolecular glycan-protein bonds stabilize the Fc region in isolated and complexed IgG1. Moreover, HDX-MS data evince that the Fab domain has glycan-protein binding contacts within the IgG1-FcγRI complex.

Methemoglobin formation in mutant hemoglobin α chains: electron transfer parameters and rates

Hemoglobin-mediated transport of dioxygen (O₂) critically depends on the stability of the reduced (Fe²⁺) form of the heme cofactors. Some protein mutations stabilize the oxidized (Fe³⁺) state (methemoglobin, Hb M), causing methemoglobinemia, and can be lethal above 30%. The majority of the analyses of factors influencing Hb oxidation are retrospective and give insights only for inner-sphere mutations of heme (His58, His87). Herein, we report the first all-atom molecular dynamics simulations on both redox states and calculations of the Marcus electron transfer (ET) parameters for the α chain Hb oxidation and reduction rates for Hb M. The Hb wild-type (WT) and most of the studied α chain variants maintain globin structure except the Hb M Iwate (H87Y). The mutants forming Hb M tend to have lower redox potentials and thus stabilize the oxidized (Fe³⁺) state (in particular, the Hb Miyagi variant with K61E mutation). Solvent reorganization (λ_solv 73-96%) makes major contributions to reorganization free energy, whereas protein reorganization (λ_prot) accounts for 27-30% except for the Miyagi and J-Buda variants (λ_prot ∼4%). Analysis of heme-solvent H-bonding interactions among variants provide insights into the role of Lys61 residue in stabilizing the Fe²⁺ state. Semiclassical Marcus ET theory-based calculations predict experimental k_ET for the Cyt b5-Hb complex and provide insights into relative reduction rates for Hb M in Hb variants. Thus, our methodology provides a rationale for the effect of mutations on the structure, stability, and Hb oxidation reduction rates and has potential for identification of mutations that result in methemoglobinemia.

Improving the speed of volumetric density map generation via cubic spline interpolation

Visualizing data generated from molecular dynamics simulations can be difficult, particularly when there can be thousands to millions of trajectory frames. The creation of a 3D grid of atomic density (i.e. a volumetric map) is one way to easily view the long-time average behavior of a system. One way to generate volumetric maps is by approximating each atom with a Gaussian function centered on that atom and spread over neighboring grid cells. However the calculation of the Gaussian function requires evaluation of the exponential function, which is computationally costly. Here we report on speeding up the calculation of volumetric maps from molecular dynamics trajectory data by replacing the expensive exponential function evaluation with an approximation using interpolating cubic splines. We also discuss the errors involved in this approximation, and recommend settings for volumetric map creation based on this.

Using Molecular Dynamics in the Study of Molecularly Imprinted Polymers

Molecular dynamics (MD) simulations of prepolymerization mixtures can provide detailed insights concerning the molecular-level mechanisms underlying the performance of molecularly imprinted polymers (MIPs) and can be used for the in silico screening of candidate polymer systems. Here, we describe the use of MD simulations of all-atom, all-component MIP prepolymerization mixtures and procedures for the evaluation of the simulation data using the Amber simulation software suite.