Structural analysis and inhibition capacity of dioscorin protein yam: Theoretical study by molecular docking and dynamic simulations.

Proteases such as trypsins in the gut of Spodoptera frugiperda are responsible for breaking down dietary proteins into amino acids necessary for insect growth and development. In this study, we characterized the insecticidal potential of dioscorin, the storage protein of yam (Dioscorea alata), using molecular docking and molecular dynamics simulations to determine the interactions between trypsin enzymes and the protein inhibitor dioscorin. To achieve this, we used the three-dimensional structures of the trypsin-like digestive enzymes of S. frugiperda, a pest of corn and cotton, as receptors or target molecules. We performed protein-protein docking using Cluspro software, estimation of the binding free energy, and information on the dynamic and time-dependent behavior of dioscorin-trypsin complexes using the NAMD package. Our computational analysis showed that dioscorin can bind to the digestive trypsins of S. frugiperda, as confirmed by the affinity energy values (-1022.4 to -1236.9), stability of the complexes during the simulation trajectory, and binding free energy values between -57.3 and -66.9 kcal/mol. Additionally, dioscorin uses two reactive sites to bind trypsin, but the largest contribution to the interaction energy is made by amino acid residues between amino acid backbone positions 8-14 by hydrogen bonds, hydrophobic, and Van der Waals (VdW) interactions. VdW is the energy that makes the greatest contribution to the binding energy. Collectively, our findings demonstrate, for the first time, the binding capacity of the yam protein dioscorin to the digestive trypsin of S. frugiperda. These promising results suggest a possible bioinsecticide action of dioscorin.

Atomistic mechanism of leptin and leptin-receptor association.

The leptin-leptin receptor complex is at the very core of energy homeostasis and immune system regulation, among many other functions. In this work, we built homology models of leptin and the leptin binding domain (LBD) of the receptor from humans and mice. Docking analyses were used to obtain the coordinates of the native leptin-LBD complexes and a mixed heterodimer formed by human leptin and mouse LBD. Molecular dynamics (MD) simulations were performed using all models (monomers and heterodimers) as initial coordinates and the GROMACS program. The overall structural and dynamical behaviors are similar for the three complexes. Upon MD simulations, several new interactions appear. In particular, hydrophobic interactions, with more than 90% persistence, seem to be the most relevant for the stability of the dimers, as well as the pair formed by Asp85Lep and Arg468LBD. This in silico analysis provides structural and dynamical information, at the atomistic level, about the mechanism of leptin-LBD complex formation and leptin receptor activation. This knowledge might be used in the rational drug design of therapeutics to modulate leptin signaling.Communicated by Ramaswamy H. Sarma.

Higher infectivity of the SARS-CoV-2 new variants is associated with K417N/T, E484K, and N501Y mutants: An insight from structural data.

The evolution of the SARS-CoV-2 new variants reported to be 70% more contagious than the earlier one is now spreading fast worldwide. There is an instant need to discover how the new variants interact with the host receptor (ACE2). Among the reported mutations in the Spike glycoprotein of the new variants, three are specific to the receptor-binding domain (RBD) and required insightful scrutiny for new therapeutic options. These structural evolutions in the RBD domain may impart a critical role to the unique pathogenicity of the SARS-CoV-2 new variants. Herein, using structural and biophysical approaches, we explored that the specific mutations in the UK (N501Y), South African (K417N-E484K-N501Y), Brazilian (K417T-E484K-N501Y), and hypothetical (N501Y-E484K) variants alter the binding affinity, create new inter-protein contacts and changes the internal structural dynamics thereby increases the binding and eventually the infectivity. Our investigation highlighted that the South African (K417N-E484K-N501Y), Brazilian (K417T-E484K-N501Y) variants are more lethal than the UK variant (N501Y). The behavior of the wild type and N501Y is comparable. Free energy calculations further confirmed that increased binding of the spike RBD to the ACE2 is mainly due to the electrostatic contribution. Further, we find that the unusual virulence of this virus is potentially the consequence of Darwinian selection-driven epistasis in protein evolution. The triple mutants (South African and Brazilian) may pose a serious threat to the efficacy of the already developed vaccine. Our analysis would help to understand the binding and structural dynamics of the new mutations in the RBD domain of the Spike protein and demand further investigation in in vitro and in vivo models to design potential therapeutics against the new variants.

Characterization of the adipogenic protein E4orf1 from adenovirus 36 through an in silico approach.

Adenovirus 36 (Ad-36) is related to human obesity due to its adipogenic activity mediated by the early 4 open reading frame 1 (E4orf1) protein. Mechanisms underlying the adipogenic effect of E4orf1 are not completely understood; however, the proliferation and differentiation of fat cells are increased through the activation of the phosphatidyl inositol 3 kinase pathway by binding proteins containing PDZ domain. This study characterized E4orf1 tridimensional structure and analyzed its interactions with PDZ domain-containing proteins in order to provide new information about the behavior of this viral protein and its targets, which could provide an interesting druggable target for obesity-related cardiometabolic alterations. In silico strategies such as homology modeling, docking, and molecular dynamics (MD) were used to study the interaction of E4orf1 with five PDZ domains of disk large homolog 1 (PDZ-1 and PDZ-2), membrane-associated guanylate kinase 1 (PDZ-3), and multi-PDZ domain protein 1 (PDZ-7 and PDZ-10). Mutagenesis analysis of selected residues was performed to evaluate their effects on the stabilization of E4orf1:PDZ complexes. MD simulations showed that the E4orf1:PDZ10 complex was more stable than the others ones. The highly hydrophobic residues at the C-terminal region (114-125) of the E4orf1 are essential in the initial phase stabilization of the complexes. Moreover, the residues 80-85 in the core region contribute to longer stabilization of the E4orf1:PDZ10 complex, a result that was confirmed by in silico mutagenesis. In conclusion, E4orf1 forms a stable complex with PDZ10 domain, and the residues 80-85 are of particular importance. The characterization of E4orf1 interactions with PDZ domains provides an initial approach to discover druggable targets for Ad-36-induced obesity.

Novel multi-epitope protein containing conserved epitopes from different Leishmania species as potential vaccine candidate: Integrated immunoinformatics and molecular dynamics approach.

Leishmaniosis, caused by intracellular parasites of the genus Leishmania, has become a serious public health problem around the world, and for which there are currently extensive limitations. In this work, a theoretical model was proposed for the development of a multi-epitope vaccine. The protein GP63 of the parasite was selected for epitopes prediction, due to its important biological role for the infection process and abundance. IEDB tools were used to determine epitopes B and T in Leishmania braziliensis; besides, other conserved epitopes in three species were selected. To improve immunogenicity, 50S ribosomal protein L7 / L12 (ID: P9WHE3) was used as a domain of adjuvant in the assembly process. The folding arrangement of the vaccine was obtained through homologous modeling multi-template with MODELLER v9.21, and a Ramachandran plot analysis was done. Furthermore, physicochemical properties were described with the ProtParam tool and secondary structure prediction combining GOR-IV and SOPMA tools. Finally, a molecular dynamics simulation (50 ns) was performed to establish flexibility and conformational changes. The analysis of the results indicates high conservancy in the epitopes predicted among the four species. Moreover, Ramachandran plot, physicochemical parameters, and secondary structure prediction suggest a stable conformation of the vaccine, after a minimum conformational change that was evaluated with the free energy landscape. The conformational change does not drive any substantial change for epitope exposition on the surface. The vaccine proposed could be tested experimentally to guide new approaches in the development of pan-vaccines; vaccines with regions conserved in multiple species.

Molecular Modeling of Structures and Interaction of Human Corticotropin-Releasing Factor (CRF) Binding Protein and CRF Type-2 Receptor.

The corticotropin-releasing factor (CRF) system is a key mediator of the stress response and addictive behavior. The CRF system includes four peptides: The CRF system includes four peptides: CRF, urocortins I-III, CRF binding protein (CRF-BP) that binds CRF with high affinity, and two class B G-protein coupled receptors CRF1R and CRF2R. CRF-BP is a secreted protein without significant sequence homology to CRF receptors or to any other known class of protein. Recently, it has been described a potentiation role of CRF-BP over CRF signaling through CRF2R in addictive-related neuronal plasticity and behavior. In addition, it has been described that CRF-BP is capable to physically interact specifically with the α isoform of CRF2R and acts like an escort protein increasing the amount of the receptor in the plasma membrane. At present, there are no available structures for CRF-BP or for full-length CRFR. Knowing and studying the structure of these proteins could be beneficial in order to characterize the CRF-BP/CRF2αR interaction. In this work, we report the modeling of CRF-BP and of full-length CRF2αR and CRF2ßR based on the recently solved crystal structures of the transmembrane domains of the human glucagon receptor and human CRF1R, in addition with the resolved N-terminal extracellular domain of CRFRs. These models were further studied using molecular dynamics simulations and protein-protein docking. The results predicted a higher possibility of interaction of CRF-BP with CRF2αR than CRF2ßR and yielded the possible residues conforming the interacting interface. Thus, the present study provides a framework for further investigation of the CRF-BP/CRF2αR interaction.

Heterodimerization of the Entamoeba histolytica EhCPADH virulence complex through molecular dynamics and protein-protein docking.

EhCPADH is a protein complex involved in the virulence of Entamoeba histolytica, the protozoan responsible for human amebiasis. It is formed by the EhCP112 cysteine protease and the EhADH adhesin. To explore the molecular basis of the complex formation, three-dimensional models were built for both proteins and molecular dynamics simulations (MDS) and docking calculations were performed. Results predicted that the pEhCP112 proenzyme and the mEhCP112 mature enzyme were globular and peripheral membrane proteins. Interestingly, in pEhCP112, the propeptide appeared hiding the catalytic site (C167, H329, N348); while in mEhCP112, this site was exposed and its residues were found structurally closer than in pEhCP112. EhADH emerged as an extended peripheral membrane protein with high fluctuation in Bro1 and V shape domains. 500 ns-long MDS and protein-protein docking predictions evidenced different heterodimeric complexes with the lowest free energy. pEhCP112 interacted with EhADH by the propeptide and C-terminal regions and mEhCP112 by the C-terminal through hydrogen bonds. In contrast, EhADH bound to mEhCP112 by 442-479 residues, adjacent to the target cell-adherence region (480-600 residues), and by the Bro1 domain (9-349 residues). Calculations of the effective binding free energy and per residue free energy decomposition showed that EhADH binds to mEhCP112 with a higher binding energy than to pEhCP112, mainly through van der Waals interactions and the nonpolar part of solvation energy. The EhADH and EhCP112 structural relationship was validated in trophozoites by immunofluorescence, TEM, and immunoprecipitation assays. Experimental findings fair agreed with in silico results.

Estudio y análisis comparativo de interacciones entre la proteina integrina con fragmentos de la proteína fibrilina-1 y fragmentos mutados de esta utilizando la metodología de docking molecular / Study and comparative analysis of interactions protein between the integrina with fragments fibrillin-1 and mutated fragments of this by using molecular docking

Resumen Objetivo: Evaluar las interacciones proteína-proteína que pueden generarse entre fragmentos de la proteína fibrilina-1, cuyas mutaciones causan el síndrome de Marfan (SM). Materiales y métodos: Se realizó una serie de cálculos docking proteína-proteína entre las macromoléculas de interés; se empleó el programa Molsoft ICM; se utilizaron las estructuras cristalinas de la proteína integrina αVβ3 y los fragmentos de la proteína fibrilina-1; además se generó una sucesión de mutaciones en la fibrilina-1, las cuales son características de pacientes con síndrome de Marfan, y posteriormente se realizó el acoplamiento molecular. Adicionalmente se determinó los aminoácidos que con mayor frecuencia estaban presentes en el sitio de interacción y su hidrofobicidad. Resultados: Se cuantificó la cantidad de aminoácidos hidrófobos presentes en las zonas de interacción producidas por los acoplamientos, teniendo en cuenta la energía del sistema; esta ponderación estuvo entre el 40 y 50 % de los aminoácidos de la zona de interacción, con un porcentaje mayor con respecto a aminoácidos neutros o cargados. En los resultados obtenidos utilizando las mutaciones realizadas sobre los fragmentos cbEGF22-TB4-cbEGF23 y cbEGF9-hyb2-cbEGF10 de la fibrilina-1 se encontró que no se ubicaron en zonas cercanas al sitio de interacción en la mayoría de los casos. Conclusiones: Las interacciones entre los fragmentos de fibrilina-1, y estos con respecto a la integrina, mostraron en sus zonas de interacción la presencia mayoritariamente de aminoácidos hidrofóbicos, que es lo esperado normalmente.

Abstract Objective: To assess protein-protein interactions between fragments of fibrillin-1 protein, whose mutations cause Marfan syndrome (MS). Materials and Methods: We performed a series of protein-protein docking calculations between the macromolecules of interest for this purpose was used Molsoft ICM program. We used the crystal structures of αVβ3 Integrin protein and fragments of fibrillin-1 protein also were generated mutations in the fibrillin-1, which are characteristic in patients with Marfan syndrome and subsequently to the molecular docking. We determined the amino acids most often present at the site of interaction and its hydrophobicity. Results: The amount of hydrophobic amino acids present in the areas of interaction given by the couplings was quantified. Given the energy of the system, was between 40 and 50% of the amino acids of the interaction zone, with a higher proportion relative to charged or neutral amino acids. In the results obtained using the mutations performed on fragments cbEGF23 cbEGF22-TB4-and-cbEGF10 cbEGF9-HYB2 of fibrillin-1, was found they were not placed in areas near the site of interaction in most cases. Conclusion: The interaction between fragments of fibrillin-1, and those with respect to their integrin showed the presence interaction zones mostly hydrophobic amino acids, which are normally expected.