Search

Cancer is one of the leading causes of death in the U.S., and tumorous cancers such as cervical, lung, breast, and ovarian cancers are the most common types. APOBEC3B is a nonessential cytidine deaminase found in humans and theorized to defend against viral infection. However, overexpression of APOBEC3B is linked to cancer in humans, which makes APOBEC3B a potential cancer treatment target through competitive inhibition for several tumorous cancers. Computational studies can help reveal a small molecule inhibitor using high-throughput virtual screening of millions of candidates with relatively little cost. This study aims to narrow the field of potential APOBEC3B inhibition candidates for future in vitro assays and provide an effective scaffold for drug design studies. Another goal of this project is to provide critical amino acid targets in the active site for future drug design studies. This study simulated 7.8 million drug candidates using high-throughput virtual screening and further processed the top scoring 241 molecules from AutoDock Vina, DOCK 6, and de novo design. Using virtual screening, de novo design, and molecular dynamics simulations, a competitive inhibitor candidate was discovered with an average binding free energy score of -46.03 kcal/mol, more than 10 kcal/mol better than the substrate control (dCMP). These results indicate that this molecule (or a structural derivative) may be an effective inhibitor of APOBEC3B and prevent host genome mutagenesis resulting from protein overexpression. Another important finding is the confirmation of essential amino acid targets, such as Tyr250 and Gln213 within the active site of APOBEC3B. Therefore, study used novel computational methods to provide a theoretical scaffold for future drug design studies that may prove useful as a treatment for epithelial cancers.Communicated by Ramaswamy H. Sarma.

The COVID-19 pandemic was declared due to the spread of the novel coronavirus, SARS-CoV-2. Viral infection is caused by the interaction between the SARS-CoV-2 receptor binding domain (RBD) and the human ACE2 receptor (hACE2). Previous computational studies have identified repurposed small molecules that target the RBD, but very few have screened drugs in the RBD-hACE2 interface. When studies focus solely on the binding affinity between the drug and the RBD, they ignore the effect of hACE2, resulting in an incomplete analysis. We screened ACE inhibitors and previously identified SARS-CoV-2 inhibitors for binding to the RBD-hACE2 interface, and then conducted 500 ns of unrestrained molecular dynamics (MD) simulations of fosinopril, fosinoprilat, lisinopril, emodin, diquafosol, and physcion bound to the interface to assess the binding characteristics of these ligands. Based on MM-GBSA analysis, all six ligands bind favorably in the interface and inhibit the RBD-hACE2 interaction. However, when we repeat our simulation by first binding the drug to the RBD before interacting with hACE2, we find that fosinopril, fosinoprilat, and lisinopril result in a strongly interacting trimeric complex (RBD-drug-hACE2). Hydrogen bonding and pairwise decomposition analyses further suggest that fosinopril is the best RBD inhibitor. However, when lisinopril is bound, it stabilizes the trimeric complex and, therefore, is not an ideal potential drug candidate. Overall, these results reveal important atomistic interactions critical to the binding of the RBD to hACE2 and highlight the significance of including all protein partners in the evaluation of a potential drug candidate., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Carter, Airas, Gladden, Miller and Parish.)

Background: Microtubule-associated proteins regulate the dynamics, organization, and function of microtubules, impacting a number of vital cellular processes. CRMPs have been shown to control microtubule assembly and axon outgrowth during neuronal differentiation. While many microtubule-associated proteins have been linked to roles in cell division and neuronal development, it is still unclear the complement that control the formation of parallel microtubule arrays in epithelial cells., Results: Here we show through time-lapse DIC microscopy that Caenorhabditis elegans embryos homozygous for the weak loss-of-function allele unc-33(e204) progress more slowly through epidermal morphogenesis, while animals homozygous for strong loss-of-function alleles exhibit more embryonic lethality. Identification of two novel missense mutations in unc-33(e572), Val476Gly, and Ser731Thr, lead to computational approaches to determine the potential effects of these changes on UNC-33/CRMP structure. Molecular dynamics simulations show that for Asp389Asn and Arg502His, two other known missense mutations, local changes in protein-protein hydrogen bonding affect the stability of the protein. However, the Val476Gly/Ser731Thr combination does not alter the structure or energetics of UNC-33 drastically when compared to the wild-type protein., Conclusions: These results support a novel role for UNC-33/CRMP in C. elegans epidermal development and shed light on how individual amino acid changes cause a loss-of-function in UNC-33., (© 2022 American Association for Anatomy.)

Photoacoustic (PA) imaging has proved versatile for many biomedical applications from drug delivery tracking to disease diagnostics and postoperative surveillance. It recently emerged as a tool for accurate and real-time heparin monitoring to avoid bleeding complications associated with anticoagulant therapy. However, molecular-dye-based application is limited by high concentration requirements, photostability, and a strong background hemoglobin signal. We developed polydopamine nanocapsules (PNCs) via supramolecular templates and loaded them with molecular dyes for enhanced PA-mediated heparin detection. Depending on surface charge, the dye-loaded PNCs undergo disassembly or aggregation upon heparin recognition: both experiments and simulation have revealed that the increased PA signal mainly results from dye-loaded PNC-heparin aggregation. Importantly, Nile blue (NB)-loaded PNCs generated a 10-fold higher PA signal than free NB dye, and such PNC enabled the direct detection of heparin in a clinically relevant therapeutic window (0-4 U/mL) in whole human blood ( R 2 = 0.91). Furthermore, the PA signal of PNC@NB obtained from 17 patients linearly correlated with ACT values ( R 2 = 0.73) and cumulative heparin ( R 2 = 0.83). This PNC-based strategy for functional nanocapsules offers a versatile engineering platform for robust biomedical contrast agents and nanocarriers.

Deposition of amyloid β (Aβ) fibrils in the brain is a key pathologic hallmark of Alzheimer's disease. A class of polyphenolic biflavonoids is known to have anti-amyloidogenic effects by inhibiting aggregation of Aβ and promoting disaggregation of Aβ fibrils. In the present study, we further sought to investigate the structural basis of the Aβ disaggregating activity of biflavonoids and their interactions at the atomic level. A thioflavin T (ThT) fluorescence assay revealed that amentoflavone-type biflavonoids promote disaggregation of Aβ fibrils with varying potency due to specific structural differences. The computational analysis herein provides the first atomistic details for the mechanism of Aβ disaggregation by biflavonoids. Molecular docking analysis showed that biflavonoids preferentially bind to the aromatic-rich, partially ordered N-termini of Aβ fibril via the π-π interactions. Moreover, docking scores correlate well with the ThT EC 50 values. Molecular dynamic simulations revealed that biflavonoids decrease the content of β-sheet in Aβ fibril in a structure-dependent manner. Hydrogen bond analysis further supported that the substitution of hydroxyl groups capable of hydrogen bond formation at two positions on the biflavonoid scaffold leads to significantly disaggregation of Aβ fibrils. Taken together, our data indicate that biflavonoids promote disaggregation of Aβ fibrils due to their ability to disrupt the fibril structure, suggesting biflavonoids as a lead class of compounds to develop a therapeutic agent for Alzheimer's disease.

Recent experiments indicate that the C-Jun amino-terminal kinase-interacting protein 1 (JIP1) binds to and activates the c-Jun N-terminal kinase (JNK) protein. JNK is an integral part of cell apoptosis, and misregulation of this process is a causative factor in diseases such as Alzheimer's disease (AD), obesity, and cancer. It has also been shown that JIP1 may increase the phosphorylation of tau by facilitating the interaction between the tau protein and JNK, which could also be a causative factor in AD. Very little is known about the structure and dynamics of JIP1; however, the amino acid composition of the first 350 residues suggests that it contains an intrinsically disordered region. Molecular dynamics (MD) simulations using AMBER 14 were used to study the structure and dynamics of a functionally active JIP1 10mer fragment to better understand the solution behavior of the fragment. Two microseconds of unbiased MD was performed on the JIP1 10mer fragment in 10 different seeds for a total of 20 μs of simulation time, and from this, seven structurally stable conformations of the 10mer fragment were identified via classical clustering. The 10mer ensemble was also used to build a Markov state model (MSM) that identified four metastable states that encompassed six of the seven conformational families identified by classical dimensional reduction. Based on this MSM, conformational interconversions between the four states occur via two dominant pathways with probability fluxes of 55 and 44% for each individual pathway. Transitions between the initial and final states occur with mean first passage times of 31 (forward) and 16 (reverse) μs.

Shigella is the causative agent of bacillary dysentery and is responsible for an estimated 165 million infections and 600,000 deaths annually. Like many Gram-negative pathogens, Shigella relies on a type three secretion system (T3SS) to initiate and sustain infection by directly injecting effector proteins into host cells. Protein secretion through the needle-like injectisome and overall Shigella virulence rely on the T3SS ATPase Spa47, making it a likely means for T3SS regulation and an attractive target for therapeutic small molecule inhibitors. Here, we utilize a recently solved 2.15 Å crystal structure of Spa47 to computationally screen 7.6 million drug-like compounds for candidates which avoid the highly conserved active site by targeting a distal, but critical, interface between adjacent protomers of the Spa47 homohexamer. Ten of the top inhibitor candidates were characterized, identifying novel Spa47 inhibitors that reduce in vitro ATPase activity by as much as 87.9 ± 10.5% with IC 50 's as low as 25 ± 20 μM and reduce in vivo Shigella T3SS protein secretion by as much as 94.7 ± 3.0%. Kinetic analyses show that the inhibitors operate through a noncompetitive mechanism that likely supports the inhibitors' low cytotoxicity, as they avoid off-target ATPases involved in either Shigella or mammalian cell metabolism. Interestingly, the inhibitors display nearly identical inhibition profiles for Spa47 and the T3SS ATPases EscN from E. coli and FliI from Salmonella. Together, the results of this study provide much-needed insight into T3SS ATPase inhibition mechanisms and a strong platform for developing broadly effective cross-pathogen T3SS ATPase inhibitors.

We recently reported a real-time method to measure heparin in human whole blood based on the photoacoustic change of methylene blue (MB). Intriguingly, the MB behaved unlike other "turn on" photoacoustic probes-the absorbance decreased as the photoacoustic signal increased. The underlying mechanism was not clear and motivated this study. We studied the binding mechanism of MB and heparin in water and phosphate buffer saline (PBS) with both experimental and computational methods. We found that the photoacoustic enhancement of the MB-heparin mixture was a result of MB-heparin aggregation due to charge neutralization and resulting sequestration of MB in these aggregates. The sequestration of MB in the MB-heparin aggregates led to decreased absorbance-there was simply less free dye in solution to absorb light. The highest photoacoustic signal and aggregation occurred when the number of negatively charged sulfate groups on heparin was approximately equal to the number of positively charged MB molecule. The MB-heparin aggregates dissociated when there were more sulfated groups from heparin than MB molecules because of the electrostatic repulsion between negatively charged sulfate groups. PBS facilitated MB dimer formation regardless of heparin concentration and reprecipitated free MB in aggregates due to ionic strength and ionic shielding. Further molecular dynamics experiments found that binding of heparin occurred at the sulfates and glucosamines in heparin. Phosphate ions could interact with the heparin via sodium ions to impair the MB-heparin binding. Finally, our model found 3.7-fold more MB dimerization upon addition of heparin in MB solution confirming that heparin facilitates MB aggregation. We conclude that the addition of heparin in MB decreases the absorbance of the sample because of MB-heparin aggregation leading to fewer MB molecules in solution; however, the aggregation also increases the PA intensity because the MB molecules in the MB-heparin aggregate have reduced degrees of freedom and poor heat transfer to solvent.

Optimizing DNA polymerases for a broad range of tasks requires an understanding of the factors influencing polymerase fidelity, but many details of polymerase behavior remain unknown, especially in the presence of mismatched nascent base pairs. Using molecular dynamics, the large fragment of Bacillus stearothermophilus DNA polymerase I is simulated in the presence of all 16 possible standard nucleoside triphosphate-template (dNTP-dN) pairs, including four Watson-Crick pairs and 12 mismatches. The precatalytic steps of nucleotide addition from nucleotide insertion to immediately preceding catalysis are explored using three starting structures representing different stages of nucleotide addition. From these simulations, interactions between dNTPs and the DNA-protein complex formed by the polymerase are elucidated. Patterns of large-scale conformational shifts, classification of nucleotide pairs based on composition, and investigation of the roles of residues interacting with dNTPs are completed on 50+ μs of simulation. The role of molecular dynamics in studies of polymerase behavior is discussed.

The accuracy of computational models for simulating biomolecules under specific solution pH conditions is critical for properly representing the effect of pH in biological processes. Constant pH (CpH) simulations involving implicit solvent using the AMBER software often incorrectly estimate pK a values of aspartate and glutamate residues due to large effective radii stemming from the presence of dummy protons. These inaccuracies stem from problems in the sampled ensembles of titratable residues that can influence other observable pH-dependent behavior, such as conformational change. We investigate new radii assignments for atoms in titratable residues with carboxylate groups to mitigate the systematic overestimation in the current method. We find that decreased carboxylate radii correspond with increased agreement with experimentally derived pK a values for residues in hen egg-white lysozyme and Δ+PHS variants of staphylococcal nuclease (SNase) and improved conformation state sampling compared to experimentally described expectations of native-like structure. Our CpH simulations suggest that decreasing the effective radii of these carboxylate groups is essential for eliminating a significant source of systematic error that hurts the accuracy of both conformational and protonation state sampling with implicit solvent.

The human immunodeficiency virus (HIV) infects healthy human cells by binding to the glycoprotein cluster of differentiation 4 receptors on the surface of helper T-cells, along with either of two chemokine receptors, CC chemokine receptor type 5 (CCR5) or C-X-C chemokine receptor type 4 (CXCR4). Recently, a pyrazolo-piperdine ligand was synthesized and the corresponding biological data showed good binding to both chemokine receptors, effectively blocking HIV-1 entry. Here, we exhaustively assess the atomistic binding interactions of this compound with both CCR5 and CXCR4, and we find that binding is driven by π-stacking interactions between aromatic rings on the ligand and receptor residues, as well as electrostatic interactions involving the protonated piperidine nitrogen. However, these favorable binding interactions were partially offset by unfavorable desolvation of active site glutamates and aspartates, prompting our proposal of a new, synthetically-accessible derivative designed to increase the electrostatic interactions without compromising the π-stacking features., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Dynemicin A has the ability to undergo the Bergman cyclization, forming a para-benzyne moiety with the ability to induce DNA strand scission. This property of dynemicin A makes it a promising anti-tumor agent. Past research has shown conclusively that dynemicin A binds to and abstracts a hydrogen atom (H5') from the DNA backbone, but the molecular mechanism of the binding event is not fully understood. We have used AMBER Molecular Dynamics simulations to investigate the dynamics associated with the reaction mechanisms. Previously, two binding mechanisms have been proposed, of which the second is more supported: (1) dynemicin A intercalates between two base pairs and directly abstracts a hydrogen atom from DNA and (2) dynemicin A inserts into the minor groove and directly abstracts a hydrogen atom from DNA. We propose a third mechanism, where dynemicin A intercalates, then undergoes a proximate, intramolecular hydrogen atom abstraction (internal abstraction). While not studied here, the resulting radical would then subsequently abstract a hydrogen atom from DNA., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Mutations in the amyloid precursor protein (APP) are responsible for the formation of amyloid-β peptides. These peptides play a role in Alzheimer's and other dementia-related diseases. The cargo binding domain of the kinesin-1 light chain motor protein (KLC1) may be responsible for transporting APP either directly or via interaction with C-jun N-terminal kinase-interacting protein 1 (JIP1). However, to date there has been no direct experimental or computational assessment of such binding at the atomistic level. We used molecular dynamics and free energy estimations to gauge the affinity for the binary complexes of KLC1, APP, and JIP1. We find that all binary complexes (KLC1:APP, KLC1:JIP1, and APP:JIP1) contain conformations with favorable binding free energies. For KLC1:APP the inclusion of approximate entropies reduces the favorability. This is likely due to the flexibility of the 42-residue APP protein. In all cases we analyze atomistic/residue driving forces for favorable interactions. Proteins 2017; 85:221-234. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

DNA polymerases must quickly and accurately distinguish between similar nucleic acids to form Watson-Crick base pairs and avoid DNA replication errors. Deoxynucleoside triphosphate (dNTP) binding to the DNA polymerase active site induces a large conformational change that is difficult to characterize experimentally on an atomic level. Here, we report an X-ray crystal structure of DNA polymerase I bound to DNA in the open conformation with a dNTP present in the active site. We use this structure to computationally simulate the open to closed transition of DNA polymerase in the presence of a Watson-Crick base pair. Our microsecond simulations allowed us to characterize the key steps involved in active site assembly, and propose the sequence of events involved in the prechemistry steps of DNA polymerase catalysis. They also reveal new features of the polymerase mechanism, such as a conserved histidine as a potential proton acceptor from the primer 3'-hydroxyl., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Assembly of polymerase chain reactions at room temperature can sometimes lead to low yields or unintentional products due to mispriming. Mutation of isoleucine 707 to leucine in DNA polymerase I from Thermus aquaticus substantially decreases its activity at room temperature without compromising its ability to amplify DNA. To understand why a conservative change to the enzyme over 20 Å from the active site can have a large impact on its activity at low temperature, we solved the X-ray crystal structure of the large (5'-to-3' exonuclease-deleted) fragment of Taq DNA polymerase containing the cold-sensitive mutation in the ternary (E-DNA-ddNTP) and binary (E-DNA) complexes. The I707L KlenTaq1 ternary complex was identical to the wild-type in the closed conformation except for the mutation and a rotamer change in nearby phenylalanine 749, suggesting that the enzyme should remain active. However, soaking out of the nucleotide substrate at low temperature results in an altered binary complex made possible by the rotamer change at F749 near the tip of the polymerase O-helix. Surprisingly, two adenosines in the 5'-template overhang fill the vacated active site by stacking with the primer strand, thereby blocking the active site at low temperature. Replacement of the two overhanging adenosines with pyrimidines substantially increased activity at room temperature by keeping the template overhang out of the active site, confirming the importance of base stacking. These results explain the cold-sensitive phenotype of the I707L mutation in KlenTaq1 and serve as an example of a large conformational change affected by a conservative mutation.

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme:DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.

Chagas disease (or American trypanosomiasis) is a deadly tropical disease that affects millions of people worldwide, primarily in rural regions of South America. Trypanosoma cruzi, the parasitic cause of Chagas disease, possesses a membrane-anchored trans-sialidase enzyme that transfers sialic acids from the host cell surface to the parasitic cell surface, allowing T. cruzi to effectively evade the host's immune system. This enzyme has no analogous human counterpart and thus has become an interesting drug target to combat the parasite. Recent computational efforts have improved our knowledge about the enzyme's structure, dynamics and catalyzed reaction. Many compounds have been tested against trans-sialidase activity, but no strong inhibitors have been identified yet. The current lack of drugs for Chagas disease necessitates more R&D into the design and discovery of strong inhibitors of T. cruzi trans-sialidase.

Chagas' is a fatal disease that affects millions of people worldwide. The lack of safe and effective treatments for Chagas' highlights the need for the discovery of new drugs to fight the disease. Trypanosoma cruzi, the parasitic cause of Chagas' disease, synthesizes a trans-sialidase (TcTS) enzyme responsible for the transfer of sialic acids from the host cell surface to glycoconjugates on the parasitic cell surface. TcTS has no human analogs and is vital to the life cycle of T. cruzi, making TcTS an important enzyme for drug design against Chagas' disease. We use fragment docking to generate various e-pharmacophore hypotheses depicting protein residues important for ligand binding. Virtual screening of the ZINC Clean Leads database with more than 4 million compounds using the e-pharmacophore models found 82 potential inhibitors of TcTS. Molecular dynamics and free energy of binding calculations were used to rank the compounds based on their affinity for TcTS. Two compounds-ZINC13359679 and ZINC02576132-were found to be the most promising lead candidates for TcTS inhibition, and their binding modes are analyzed in detail., (Copyright © 2013 Elsevier Inc. All rights reserved.)

MM-PBSA is a post-processing end-state method to calculate free energies of molecules in solution. MMPBSA.py is a program written in Python for streamlining end-state free energy calculations using ensembles derived from molecular dynamics (MD) or Monte Carlo (MC) simulations. Several implicit solvation models are available with MMPBSA.py, including the Poisson-Boltzmann Model, the Generalized Born Model, and the Reference Interaction Site Model. Vibrational frequencies may be calculated using normal mode or quasi-harmonic analysis to approximate the solute entropy. Specific interactions can also be dissected using free energy decomposition or alanine scanning. A parallel implementation significantly speeds up the calculation by dividing frames evenly across available processors. MMPBSA.py is an efficient, user-friendly program with the flexibility to accommodate the needs of users performing end-state free energy calculations. The source code can be downloaded at http://ambermd.org/ with AmberTools, released under the GNU General Public License.

Substrate ingress and product egress from the active site of urease is tightly controlled by an active-site flap. Molecular dynamics simulations of urease have revealed a previously unobserved wide-open flap state that, unlike the well-characterized closed and open states, allows ready access to the metal cluster in the active site. This state is easily reached from the open state via low free energy barriers. Additionally, we have found that even when the flap is closed, a region of the binding pocket is solvent-exposed, leading to the hypothesis that it may act as a substrate/product reservoir. The newly identified wide-open state offers further opportunities for small-molecule drug discovery by defining a more extensive active-site pocket than has been previously described.

Recognition of the human immunodeficiency virus Rev-responsive element (RRE) RNA by the Rev protein is an essential step in the viral life cycle. Formation of the Rev-RRE complex signals nucleocytoplasmic export of unspliced and partially spliced viral RNA. Essential components of the complex have been localized to a minimal arginine-rich Rev peptide and stem IIB of RRE. In vitro selection studies have identified a synthetic peptide known as RSG 1.2 that binds with better specificity and affinity to RRE than the Rev peptide. NMR structures of both peptide-RNA complexes of Rev and RSG 1.2 bound to RRE stem IIB have been solved and reveal gross structural differences between the two bound complexes. Molecular dynamics simulations of the Rev and RSG 1.2 peptides in complex with RRE stem IIB have been simulated to better understand on an atomic level how two arginine-rich peptides of similar length recognize the same sequence of RNA with such different structural motifs. While the Rev peptide employs some base-specific hydrogen bonding for recognition of RRE, shape recognition, through contact with the sugar-phosphate backbone, and cation-pi interactions are also important. Molecular dynamics simulations suggest that RSG 1.2 binds more tightly to the RRE sequence than Rev by forming more base-specific contacts, using water to mediate peptide-RNA contacts, and is held in place by a strong salt bridge network spanning the major groove of the RNA.