Conformational response to ligand binding of TMPRSS2, a protease involved in SARS-CoV-2 infection: Insights through computational modeling.

Thanks to the considerable research which has been undertaken in the last few years to improve our understanding of the biology and mechanism of action of SARS-CoV-2, we know how the virus uses its surface spike protein to infect host cells. The transmembrane prosthesis, serine 2 (TMPRSS2) protein, located on the surface of human cells, recognizes the cleavage site in the spike protein, leading to the release of the fusion peptide and entry of the virus into the host cells. Because of its role, TMPRSS2 has been proposed as a drug target to prevent infection by the virus. In this study, we aim to increase our understanding of TMPRSS2 using long scale microsecond atomistic molecular dynamics simulations, focusing on the conformational changes over time. The comparison between simulations conducted on the protein in the native (apo) and inhibited form (holo), has shown that in the holo form the inhibitor stabilizes the catalytic site and induces rearrangements in the extracellular domain of the protein. In turn, it leads to the formation of a new cavity in the vicinity of the ligand binding pocket that is stable in the microsecond time scale. Given the low specificity of known protease inhibitors, these findings suggest a new potential drug target site that can be used to improve TMPRSS2 specific recognition by newly designed inhibitors.

Computational Studies of SARS-CoV-2 3CLpro: Insights from MD Simulations.

Given the enormous social and health impact of the pandemic triggered by severe acute respiratory syndrome 2 (SARS-CoV-2), the scientific community made a huge effort to provide an immediate response to the challenges posed by Coronavirus disease 2019 (COVID-19). One of the most important proteins of the virus is an enzyme, called 3CLpro or main protease, already identified as an important pharmacological target also in SARS and Middle East respiratory syndrome virus (MERS) viruses. This protein triggers the production of a whole series of enzymes necessary for the virus to carry out its replicating and infectious activities. Therefore, it is crucial to gain a deeper understanding of 3CLpro structure and function in order to effectively target this enzyme. All-atoms molecular dynamics (MD) simulations were performed to examine the different conformational behaviors of the monomeric and dimeric form of SARS-CoV-2 3CLpro apo structure, as revealed by microsecond time scale MD simulations. Our results also shed light on the conformational dynamics of the loop regions at the entry of the catalytic site. Studying, at atomic level, the characteristics of the active site and obtaining information on how the protein can interact with its substrates will allow the design of molecules able to block the enzymatic function crucial for the virus.

Key Players in I-DmoI Endonuclease Catalysis Revealed from Structure and Dynamics.

Homing endonucleases, such as I-DmoI, specifically recognize and cleave long DNA target sequences (∼20 bp) and are potentially powerful tools for genome manipulation. However, inefficient and off-target DNA cleavage seriously limits specific editing in complex genomes. One approach to overcome these limitations is to unambiguously identify the key structural players involved in catalysis. Here, we report the E117A I-DmoI mutant crystal structure at 2.2 Šresolution that, together with the wt and Q42A/K120M constructs, is combined with computational approaches to shed light on protein cleavage activity. The cleavage mechanism was related both to key structural effects, such as the position of water molecules and ions participating in the cleavage reaction, and to dynamical effects related to protein behavior. In particular, we found that the protein perturbation pattern significantly changes between cleaved and noncleaved DNA strands when the ions and water molecules are correctly positioned for the nucleophilic attack that initiates the cleavage reaction, in line with experimental enzymatic activity. The proposed approach paves the way for an effective, general, and reliable procedure to analyze the enzymatic activity of endonucleases from a very limited data set, i.e., structure and dynamics.

Kinetics and mechanistic study of competitive inhibition of thymidine phosphorylase by 5-fluoruracil derivatives.

In a previous investigation, cationic liposomes formulated with new 5-FU derivatives, differing for the length of the polyoxyethylenic spacer that links the N(3) position of 5-FU to an alkyl chain of 12 carbon atoms, showed a higher cytotoxicity compared to free 5-FU, the cytotoxic effect being directly related to the length of the spacer. To better understand the correlation of the spacer length with toxicity, we carried out initial rate studies to determine inhibition, equilibrium and kinetic constants (KI, KM, kcat), and get inside inhibition activity of the 5-FU derivatives and their mechanism of action, a crucial information to design structural variations for improving the anticancer activity. The experimental investigation was supported by docking simulations based on the X-ray structure of thymidine phosphorylase (TP) from Escherichia coli complexed with 3'-azido-2'-fluoro-dideoxyuridin. Theoretical and experimental results showed that all the derivatives exert the same inhibition activity of 5-FU either as monomer and when embedded in lipid bilayer.

Molecular Packing in Langmuir Monolayers Composed of a Phosphatidylcholine and a Pyrene Lipid.

Pyrene lipids are useful tools to investigate membrane organization and intracellular lipid trafficking. The molecular interactions controlling the organization of lipid monolayers composed of a cationic amphiphile tagged with a pyrene residue and a saturated or unsaturated phospholipid, namely, 1,2-dimyristoyl-sn-glycero-3-phosphocholine and 1,2-dioleoyl-sn-glycero-3-phosphocholine, were investigated by Langmuir trough isotherms to understand how the molecular structure of the components and their relative amount affect the physicochemical properties of lipid monolayers. The obtained results show that the cationic headgroups and unsaturation of hydrophobic chains strongly affect the organization of the lipid monolayer as a function of the amount of components. On the other hand, the presence of the pyrene moiety does not seem to have a marked influence on the interaction within lipid assembly.

Modeling conformational transitions in kinases by molecular dynamics simulations: achievements, difficulties, and open challenges.

Protein kinases work because their flexibility allows to continuously switch from inactive to active form. Despite the large number of structures experimentally determined in such states, the mechanism of their conformational transitions as well as the transition pathways are not easily to capture. In this regard, computational methods can help to shed light on such an issue. However, due to the intrinsic sampling limitations, much efforts have been done to model in a realistic way the conformational changes occurring in protein kinases. In this review we will address the principal biological achievements and structural aspects in studying kinases conformational transitions and will focus on the main challenges related to computational approaches such as molecular modeling and MD simulations.

Molecular mechanisms of activation in CDK2.

Cyclin-dependent kinases (CDKs) are enzymes involved in crucial cellular processes. Their biological activity is directly linked to their high conformational variability, which involves large protein conformational rearrangements. We present here the application of an enhancing sampling technique to the study of conformational transitions between the open and closed state of CDKs. The analysis of the conformational intermediates supports the idea that the process is regulated by two important protein regions, which sequentially rearrange in order to allow the protein to reach its final conformation. Furthermore, the two paths involve additional (minor) protein rearrangements which are specific to the paths. Our results show that our procedure can provide reasonable transition pathways between the two protein forms at a very reduced computational cost. The robustness and the simplicity of our approach make it of general application to describe virtually any macromolecular conformational transitions.

Using metadynamics and path collective variables to study ligand binding and induced conformational transitions.

Large-scale conformational transitions represent both a challenge and an opportunity for computational drug design. Exploring the conformational space of a druggable target with sufficient detail is computationally demanding. However, if it were possible to fully account for target flexibility, one could exploit this knowledge to rationally design more potent and more selective drug candidates. Here, we discuss how molecular dynamics together with free energy algorithms based on Metadynamics and Path Collective Variables can be used to study both large-scale conformational transitions and ligand binding to flexible targets. We show real-life examples of how these methods have been applied in the case of cyclin-dependent kinases, a family of flexible targets that shows promise in cancer therapy.

Structural insight into peroxisome proliferator-activated receptor gamma binding of two ureidofibrate-like enantiomers by molecular dynamics, cofactor interaction analysis, and site-directed mutagenesis.

Molecular dynamics simulations were performed on two ureidofibrate-like enantiomers to gain insight into their different potency and efficacy against PPARgamma. The partial agonism of the S enantiomer seems to be due to its capability to stabilize different regions of the receptor allowing the interaction with both coactivators and corepressors as shown by fluorescence resonance energy transfer (FRET) assays. The recruitment of the corepressor N-CoR1 by the S enantiomer on two different responsive elements of PPARgamma regulated promoters was confirmed by chromatin immunoprecipitation assays. Cell-based transcription assays show that PPARgamma coactivator 1alpha (PGC-1alpha) and cAMP response element binding protein-binding protein (CBP) enhance the basal and ligand-stimulated receptor activity acting as coactivators of PPARgamma, whereas the receptor interacting protein 140 (RIP140) and the nuclear corepressor 1 (N-CoR1) repress the transcriptional activity of PPARgamma. We also tested the importance of the residue Q286 on the transcriptional activity of the receptor by site-directed mutagenesis and confirmed its key role in the stabilization of helix 12. Molecular modeling studies were performed to provide a molecular explanation for the different behavior of the mutants.

Synthesis, biological evaluation, and docking studies of N-substituted acetamidines as selective inhibitors of inducible nitric oxide synthase.

New acetamidines structurally related to N-(3-(aminomethyl)benzyl)acetamidine (1, W1400) were designed as inhibitors of inducible nitric oxide synthase (iNOS). Six compounds were found to be selective for iNOS over endothelial nitric oxide synthase (eNOS), and among them, the most active and selective compound was the N-benzylacetamidine 2. A docking study was also performed to shed light on the effects of the structural modifications on the interaction of the designed inhibitors with the NOS.

Aquation of the ruthenium-based anticancer drug NAMI-A: a density functional study.

We carried out density functional theory (DFT) calculations to investigate the thermodynamics and the kinetics of the double aquation reaction of the anticancer drug NAMI-A. Three explicit water molecules were included in the calculations to improve the PB solvation energies. Our calculations show that the chloride substitution reactions on the considered Ru(III) octahedral complex follow a dissociative interchange mechanism, I(d), passing through a loose heptacoordinate transition state. We calculated an activation enthalpy and free energy for the first aquation step of 101.5 and 103.7 kJ mol(-1), respectively, values that are in good agreement with the available experimental results. The activation enthalpy and free energy for the second aquation step were found significantly higher, 118.7 and 125.0 kJ mol(-1), again in agreement with the experimental evidence indicating a slower rate for the second aquation.

Binding of antitumor ruthenium complexes to DNA and proteins: a theoretical approach.

The thermodynamics of the binding of the antitumor ammine, amine, and immine complexes of ruthenium(II) and ruthenium(III) to DNA and peptides was studied computationally using model molecules. We performed density functional calculations on several monofunctional ruthenium complexes of the formula [Ru(NH3)5B]z+, where B is an adenine, guanine, or cytosine nucleobase or an 4-methylimidazole, a dimethylthioether, or a dimethylphosphate anion and z = 2 and 3. The pentammineruthenium fragment has been intensively studied and also constitutes a good model for a wide class of antitumor ammine, amine, and imine complexes of Ru(II) and Ru(III), while the considered bases/ligands have been chosen as models for the main binding sites of DNA, nucleobases, and phosphate backbone and proteins, histidyl, and sulfur-containing residue such as methionine or cysteine. Bond dissociation enthalpies and free energies have been calculated for all the considered metal binding sites both in the gas phase and in solution and allow building a binding affinity order for the considered nucleic acid or protein binding sites. The binding of guanine to some bifunctional complexes, [Ru(NH3)(4)Cl2], [cis-RuCl(2)(bpy)2], and [cis-RuCl(2)(azpy)2], has also been considered to evaluate the effect of a second labile chloro or aquo ligand and more realistic polypyridyl and arylazopyridine ligands.

Stereoselectivity by enantiomeric inhibitors of matrix metalloproteinase-8: new insights from molecular dynamics simulations.

Molecular Dynamics simulations in aqueous solution were performed for the matrix metalloproteinase-8 (MMP-8) free catalytic domain and for its complexes with the (R)- and (S)-[1-(4'-methoxybiphenyl-4-sulfonylamino)-2-methylpropyl] phosphonate. The 144-155 loop of the enzyme undergoes a drastic decrease of mobility once complexed with both enantiomers. The two enantiomers induce a different decrease of conformational entropy upon complexation. The higher affinity of the R-enantiomer can be related to the lower loss of conformational entropy accompanying its binding. The differences in the dynamical behavior of the protein induced by the two enantiomers are discussed at molecular level and the mode of binding of the simulated complexes is compared with that previously determined by X-ray crystallography.

Theoretical models of possible compact nucleosome structures.

Chromatin structure seems related to the DNA linker length. This paper presents a systematic search of the possible chromatin structure as a function of the linker lengths, starting from three different low-resolution molecular models of the nucleosome. Gay-Berne potential was used to evaluate the relative nucleosome packing energy. Results suggest that linker DNAs, which bridges and orientate nucleosomes, affect both the geometry and the rigidity of the global chromatin structure.

Systematic search for compact structures of telomeric nucleosomes.

Telomeres are structures functionally and structurally distinct from bulk chromatin. They are constituted of highly conserved 5-7 bp tandemly repeated units, organized into nucleosomes with short linkers, whereas the knowledge of the linker histone role in telomeric chromatin is still fragmentary. Experimental evidence suggests the structural organization of telomeric nucleosomes is different from that of the bulk chromatin. This work presents a systematic search of the telomeric nucleosome arrangements. A low-resolution molecular model was used to evaluate the relative nucleosome packing energy. Structures with favorable energy were found, reducing the possible telomeric chromatin conformations to two different three-dimensional folds.