In recent years, the awareness that pesticides can have other effects apart from generic toxicity is growing. In particular, several pieces of evidence highlight their influence on human fertility. In this study, we investigated, by a virtual screening approach, the binding between pesticides and proteins present in human gametes or associated with reproduction, in order to identify new interactions that could affect human fertility. To this aim, we prepared ligand (pesticides) and receptor (proteins) 3D structure datasets from online structural databases (such as PubChem and RCSB), and performed a virtual screening analysis using Autodock Vina. In the comparison of the predicted interactions, we found that famoxadone was predicted to bind Cellular Retinol Binding Protein-III in the retinol-binding site with a better minimum energy value of -10.4 Kcal/mol and an RMSD of 3.77 with respect to retinol (-7.1 Kcal/mol). In addition to a similar network of interactions, famoxadone binding is more stabilized by additional hydrophobic patches including L20, V29, A33, F57, L117, and L118 amino acid residues and hydrogen bonds with Y19 and K40. These results support a possible competitive effect of famoxadone on retinol binding with impacts on the ability of developing the cardiac tissue, in accordance with the literature data on zebrafish embryos. Moreover, famoxadone binds, with a minimum energy value between -8.3 and -8.0 Kcal/mol, to the IZUMO Sperm-Egg Fusion Protein, interacting with a network of polar and hydrophobic amino acid residues in the cavity between the 4HB and Ig-like domains. This binding is more stabilized by a predicted hydrogen bond with the N185 residue of the protein. A hindrance in this position can probably affect the conformational change for JUNO binding, avoiding the gamete membrane fusion to form the zygote. This work opens new interesting perspectives of study on the effects of pesticides on fertility, extending the knowledge to other typologies of interaction which can affect different steps of the reproductive process.

It is nowadays clear that RNA molecules can play active roles in several biological processes. As a result, an increasing number of RNAs are gradually being identified as potentially druggable targets. In particular, noncoding RNAs can adopt highly organized conformations that are suitable for drug binding. However, RNAs are still considered challenging targets due to their complex structural dynamics and high charge density. Thus, elucidating relevant features of drug-RNA binding is fundamental for advancing drug discovery. Here, by using Molecular Dynamics simulations, we compare key features of ligand binding to proteins with those observed in RNA. Specifically, we explore similarities and differences in terms of (i) conformational flexibility of the target, (ii) electrostatic contribution to binding free energy, and (iii) water and ligand dynamics. As a test case, we examine binding of the same ligand, namely riboflavin, to protein and RNA targets, specifically the riboflavin (RF) kinase and flavin mononucleotide (FMN) riboswitch. The FMN riboswitch exhibited enhanced fluctuations and explored a wider conformational space, compared to the protein target, underscoring the importance of RNA flexibility in ligand binding. Conversely, a similar electrostatic contribution to the binding free energy of riboflavin was found. Finally, greater stability of water molecules was observed in the FMN riboswitch compared to the RF kinase, possibly due to the different shape and polarity of the pockets.

Molecular Dynamics (MD) is a computational technique widely used to evaluate a molecular system\'s thermodynamic properties and conformational behavior over time. In particular, the energy analysis of a protein conformation ensemble produced though MD simulations plays a crucial role in explaining the relationship between protein dynamics and its mechanism of action. In this research work, the HINT (Hydropathic INTeractions) LogP-based scoring function was first used to handle MD trajectories and investigate the molecular basis behind the intricate PPARγ mechanism of activation. The Peroxisome Proliferator-Activated Receptor γ (PPARγ) is an emblematic example of a highly flexible protein due to the extended ω-loop delimiting the active site, and it is responsible for the receptor\'s ability to bind chemically different compounds. In this work, we focused on the PPARγ complex with Rosiglitazone, a common anti-diabetic compound and analyzed the molecular basis of the flexible ω-loop stabilization effect produced by the Oleic Acid co-binding. The HINT-based analysis of the produced MD trajectories allowed us to account for all of the energetic contributions involved in interconverting between conformational states and describe the intramolecular interactions between the flexible ω-loop and the helix H3 triggered by the allosteric binding mechanism.

Cellular behavior, cell differentiation and ontogenetic development in eukaryotes result from complex interactions between epigenetic and classic molecular genetic mechanisms, with many of these interactions still to be elucidated. Histone deacetylase enzymes (HDACs) promote the interaction of histones with DNA by compacting the nucleosome, thus causing transcriptional repression. MADS-domain transcription factors are highly conserved in eukaryotes and participate in controlling diverse developmental processes in animals and plants, as well as regulating stress responses in plants. In this work, we focused on finding out putative interactions of Arabidopsis thaliana HDACs and MADS-domain proteins using an evolutionary perspective combined with bioinformatics analyses and testing the more promising predicted interactions through classic molecular biology tools. Through bioinformatic analyses, we found similarities between HDACs proteins from different organisms, which allowed us to predict a putative protein-protein interaction between the Arabidopsis thaliana deacetylase HDA15 and the MADS-domain protein XAANTAL1 (XAL1). The results of two-hybrid and Bimolecular Fluorescence Complementation analysis demonstrated in vitro and in vivo HDA15-XAL1 interaction in the nucleus. Likely, this interaction might regulate developmental processes in plants as is the case for this type of interaction in animals.

The protein PCSK9 (proprotein convertase subtilisin/Kexin type 9) negatively regulates the recycling of LDLR (low-density lipoprotein receptor), leading to an elevated plasma level of LDL. Inhibition of PCSK9-LDLR interaction has emerged as a promising therapeutic strategy to manage hypercholesterolemia. However, the large interaction surface area between PCSK9 and LDLR makes it challenging to identify a small molecule competitive inhibitor. An alternative strategy would be to identify distal cryptic sites as targets for allosteric inhibitors that can remotely modulate PCSK9-LDLR interaction. Using several microseconds long molecular dynamics (MD) simulations, we demonstrate that on binding with LDLR, there is a significant conformational change (population shift) in a distal loop (residues 211-222) region of PCSK9. Consistent with the bidirectional nature of allostery, we establish a clear correlation between the loop conformation and the binding affinity with LDLR. Using a thermodynamic argument, we establish that the loop conformations predominantly present in the apo state of PCSK9 would have lower LDLR binding affinity, and they would be potential targets for designing allosteric inhibitors. We elucidate the molecular origin of the allosteric coupling between this loop and the LDLR binding interface in terms of the population shift in a set of salt bridges and hydrogen bonds. Overall, our work provides a general strategy toward identifying allosteric hotspots: compare the conformational ensemble of the receptor between the apo and bound states of the protein and identify distal conformational changes, if any. The inhibitors should be designed to bind and stabilize the apo-specific conformations.

For antimicrobial agents in particular, plasma protein binding (PPB) plays a pivotal role in deciphering key properties of drug candidates. Animal models are generally used in the preclinical development of new drugs to predict their effects in humans using translational pharmacokinetics/pharmacodynamics (PK/PD). Thus, we compared the protein binding (PB) of cefazolin as well as bacterial growth under various conditions in vitro. The PB extent of cefazolin was studied in human, bovine, and rat plasmas at different antibiotic concentrations in buffer and media containing 20-70% plasma or pure plasma using ultrafiltration (UF) and equilibrium dialysis (ED). Moreover, bacterial growth and time-kill assays were performed in Mueller Hinton Broth (MHB) containing various plasma percentages. The pattern for cefazolin binding to plasma proteins was found to be similar for both UF and ED. There was a significant decrease in cefazolin binding to bovine plasma compared to human plasma, whereas the pattern in rat plasma was more consistent with that in human plasma. Our growth curve analysis revealed considerable growth inhibition of Escherichia coli at 70% bovine or rat plasma compared with 70% human plasma or pure MHB. As expected, our experiments with cefazolin at low concentrations showed that E. coli grew slightly better in 20% human and rat plasma compared to MHB, most probably due to cefazolin binding to proteins in the plasma. Based on the example of cefazolin, our study highlights the interspecies differences of PB with potential impact on PK/PD. These findings should be considered before preclinical PK/PD data can be extrapolated to human patients.

A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand\'s unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.

In this study, we utilize Protein Residue Networks (PRNs), constructed using Local Spatial Pattern (LSP) alignment, to explore the dynamic behavior of Catabolite Activator Protein (CAP) upon the sequential binding of cAMP. We employed the Degree Centrality of these PRNs to investigate protein dynamics on a sub-nanosecond time scale, hypothesizing that it would reflect changes in CAP\'s entropy related to its thermal motions. We show that the binding of the first cAMP led to an increase in stability in the Cyclic-Nucleotide Binding Domain A (CNBD-A) and destabilization in CNBD-B, agreeing with previous reports explaining the negative cooperativity of cAMP binding in terms of an entropy-driven allostery. LSP-based PRNs also allow for the study of Betweenness Centrality, another graph-theoretical characteristic of PRNs, providing insights into global residue connectivity within CAP. Using this approach, we were able to correctly identify amino acids that were shown to be critical in mediating allosteric interactions in CAP. The agreement between our studies and previous experimental reports validates our method, particularly with respect to the reliability of Degree Centrality as a proxy for entropy related to protein thermal dynamics. Because LSP-based PRNs can be easily extended to include dynamics of small organic molecules, polynucleotides, or other allosteric proteins, the methods presented here mark a significant advancement in the field, positioning them as vital tools for a fast, cost-effective, and accurate analysis of entropy-driven allostery and identification of allosteric hotspots.

Since its first report in 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posed a grave threat to public health. Virus-specific countermeasures, such as vaccines and therapeutics, have been developed and have contributed to the control of the viral pandemic, which has become endemic. Nonetheless, new variants continue to emerge and could cause a new pandemic. Consequently, it is important to comprehensively understand viral evolution and the roles of mutations in viral infectivity and transmission. SARS-CoV-2 beta variant encode mutations (D614G, N501Y, E484K, and K417N) in the spike which are frequently found in other variants as well. While their individual role in viral infectivity has been elucidated against various therapeutic antibodies, it still remains unclear whether those mutations may act additively or synergistically when combined. Here, we report that N501Y mutation shows differential effect on two therapeutic antibodies tested. Interestingly, the relative importance of E484K and K417N mutations in antibody evasion varies depending on the antibody type. Collectively, these findings suggest that continuous efforts to develop effective antibody therapeutics and combinatorial treatment with multiple antibodies are more rational and effective forms of treatment.

The biological target identification process, a pivotal phase in the drug discovery workflow, becomes particularly challenging when mutations affect proteins\' mechanisms of action. COVID-19 Spike glycoprotein mutations are known to modify the affinity toward the human angiotensin-converting enzyme ACE2 and several antibodies, compromising their neutralizing effect. Predicting new possible mutations would be an efficient way to develop specific and efficacious drugs, vaccines, and antibodies. In this work, we developed and applied a computational procedure, combining constrained logic programming and careful structural analysis based on the Structural Activity Relationship (SAR) approach, to predict and determine the structure and behavior of new future mutants. \"Mutations rules\" that would track statistical and functional types of substitutions for each residue or combination of residues were extracted from the GISAID database and used to define constraints for our software, having control of the process step by step. A careful molecular dynamics analysis of the predicted mutated structures was carried out after an energy evaluation of the intermolecular and intramolecular interactions using the HINT (Hydrophatic INTeraction) force field. Our approach successfully predicted, among others, known Spike mutants.