Journal of Molecular Biology
Regular ArticleComparative Protein Modelling by Satisfaction of Spatial Restraints
Abstract
We describe a comparative protein modelling method designed to find the most probable structure for a sequence given its alignment with related structures. The three dimensional (3D) model is obtained by optimally satisfying spatial restraints derived from the alignment and expressed as probability density functions (pdfs) for the features restrained. For example, the probabilities for main-chain conformations of a modelled residue may be restrained by its residue type, main-chain conformation of an equivalent residue in a related protein, and the local similarity between the two sequences. Several such pdfs are obtained from the correlations between structural features in 17 families of homologous proteins which have been aligned on the basis of their 3D structures. The pdfs restrain Cα-Cα distances, main-chain N-O distances, main-chain and side-chain dihedral angles. A smoothing procedure is used in the derivation of these relationships to minimize the problem of a sparse database. The 3D model of a protein is obtained by optimization of the molecular pdf such that the model violates the input restraints as little as possible. The molecular pdf is derived as a combination of pdfs restraining individual spatial features of the whole molecule. The optimization procedure is a variable target function method that applies the conjugate gradients algorithm to positions of all non hydrogen atoms. The method is automated and is illustrated by the modelling of trypsin from two other serine proteinases.
References (0)
Cited by (11069)
Supreme glutathione-dependent ketosteroid isomerase in the yellow-fever transmitting mosquito Aedes aegypti
2024, Biochemical and Biophysical Research CommunicationsThe steroid hormone ecdysone is essential for the reproduction and survival of insects. The hormone is synthesized from dietary sterols such as cholesterol, yielding ecdysone in a series of consecutive enzymatic reactions. In the insect orders Lepidoptera and Diptera a glutathione transferase called Noppera-bo (Nobo) plays an essential, but biochemically uncharacterized, role in ecdysteroid biosynthesis. The Nobo enzyme is consequently a possible target in harmful dipterans, such as disease-carrying mosquitoes. Flavonoid compounds inhibit Nobo and have larvicidal effects in the yellow-fever transmitting mosquito Aedes aegypti, but the enzyme is functionally incompletely characterized. We here report that within a set of glutathione transferase substrates the double-bond isomerase activity with 5-androsten-3,17-dione stands out with an extraordinary specific activity of 4000 μmol min−1 mg−1. We suggest that the authentic function of Nobo is catalysis of a chemically analogous ketosteroid isomerization in ecdysone biosynthesis.
The role of glycosylation in non-productive adsorption of cellulase to lignin isolated from pretreated corn stover
2024, International Journal of Biological MacromoleculesGlycosylation, a general post-translational modification for fungal cellulase, has been shown to affect cellulase binding to its substrate. However, the exact impact of glycosylation on cellulase-lignin interaction remain unclear. Here, we demonstrated that the lignin isolated from tetrahydrofuran-pretreated corn stover exhibits strong adsorption capability to cellulase due to its negatively charged and porous structure. For the cellulases with varying glycosylation levels, the less-glycosylated protein showed high adsorption capability to lignin, and that trend was observed for the main cellulase components secreted by Penicillium oxilicum, including endoglucanase PoCel5B, cellobiohydrolase PoCel7A-2, and β-glucosidase PoBgl1. Additionally, N-glycan sites and motifs were examined using mass spectrometry, and protein structures with N-glycans were constructed, where PoBgl1 and PoCel7A-2 contained 13 and 1 glycosylated sites respectively. The results of molecular dynamics simulations indicated that the N-glycans impacted on the solvent-accessible surface area and secondary structure of protein, and the binding conformation of lignin fragment on cellulase, resulting in a decrease in binding energy (14 kcal/mol for PoBgl1 and 13 kcal/mol for PoCel7A-2), particularly for van der Waals and electrostatic interaction. Those findings suggested that glycosylation negatively impacted the lignin-cellulase interaction, providing a theoretical basis for the rational engineering of enzymes to reduce lignin-enzyme interaction.
Insights into nucleoside hydrolase from Leishmania donovani inhibition: A new bioaffinity chromatography-based screening assay and docking studies
2024, Bioorganic ChemistryLeishmaniasis, a group of neglected infectious diseases, encompasses a serious health concern, particularly with visceral leishmaniasis exhibiting potentially fatal outcomes. Nucleoside hydrolase (NH) has a fundamental role in the purine salvage pathway, crucial for Leishmania donovani survival, and presents a promising target for developing new drugs for visceral leishmaniasis treatment. In this study, LdNH was immobilized into fused silica capillaries, resulting in immobilized enzyme reactors (IMERs). The LdNH-IMER activity was monitored on-flow in a multidimensional liquid chromatography system, with the IMER in the first dimension. A C18 analytical column in the second dimension furnished the rapid separation of the substrate (inosine) and product (hypoxanthine), enabling direct enzyme activity monitoring through product quantification. LdNH-IMER exhibited high stability and was characterized by determining the Michaelis-Menten constant. A known inhibitor (1-(β-d-Ribofuranosyl)-4-quinolone derivative) was used as a model to validate the established method in inhibitor recognition. Screening of three additional derivatives of 1-(β-d-Ribofuranosyl)-4-quinolone led to the discovery of novel inhibitors, with compound 2a exhibiting superior inhibitory activity (Ki = 23.37 ± 3.64 µmol/L) compared to the employed model inhibitor. Docking and Molecular Dynamics studies provided crucial insights into inhibitor interactions at the enzyme active site, offering valuable information for developing new LdNH inhibitors. Therefore, this study presents a novel screening assay and contributes to the development of potent LdNH inhibitors.
Peptide translocation across asymmetric phospholipid membranes
2024, Biophysical JournalThe transport of molecules across cell membranes is vital for proper cell function and effective drug delivery. While most cell membranes naturally possess an asymmetric lipid composition, research on membrane transport predominantly uses symmetric lipid membranes. The permeation through the asymmetric membrane is then calculated as a sum of the inverse permeabilities of leaflets from symmetric bilayers. In this study, we examined how two types of amphiphilic molecules translocate across both asymmetric and symmetric membranes. Using computer simulations with both coarse-grained and atomistic force fields, we calculated the free energy profiles for the passage of model amphiphilic peptides and a lipid across various membranes. Our results consistently demonstrate that while the free energy profiles for asymmetric membranes with a small differential stress concur with symmetric ones in the region of lipid headgroups, the profiles differ around the center of the membrane. In this region, the free energy for the asymmetric membrane transitions between the profiles for two symmetric membranes. In addition, we show that peptide permeability through an asymmetric membrane cannot always be predicted from the permeabilities of the symmetric membranes. This indicates that using symmetric membranes falls short in providing an accurate depiction of peptide translocation across asymmetric membranes.
A multiscale approach reveals the molecular architecture of the autoinhibited kinesin KIF5A
2024, Journal of Biological ChemistryKinesin-1 is a microtubule motor that transports cellular cargo along microtubules. KIF5A is one of three kinesin-1 isoforms in humans, all of which are autoinhibited by an interaction between the motor and an IAK motif in the proximal region of the C-terminal tail. The C-terminal tail of KIF5A is ∼80 residues longer than the other two kinesin-1 isoforms (KIF5B and KIF5C) and it is unclear if it contributes to autoinhibition. Mutations in KIF5A cause neuronal diseases and could affect autoinhibition, as reported for a mutation that skips exon 27, altering its C-terminal sequence. Here, we combined negative-stain electron microscopy, crosslinking mass spectrometry (XL-MS) and AlphaFold2 structure prediction to determine the molecular architecture of the full-length autoinhibited KIF5A homodimer, in the absence of light chains. We show that KIF5A forms a compact, bent conformation, through a bend between coiled-coils 2 and 3, around P687. XL-MS of WT KIF5A revealed extensive interactions between residues in the motor, between coiled-coil 1 and the motor, between coiled-coils 1 and 2, with coiled-coils 3 and 4, and the proximal region of the C-terminal tail and the motor in the autoinhibited state, but not between the distal C-terminal region and the rest of the molecule. While negative-stain electron microscopy of exon-27 KIF5A splice mutant showed the presence of autoinhibited molecules, XL-MS analysis suggested that its autoinhibited state is more labile. Our model offers a conceptual framework for understanding how mutations within the motor and stalk domain may affect motor activity.
Tumor-acquired somatic mutation affects conformation to abolish ABCG2-mediated drug resistance
2024, Drug Resistance UpdatesABCG2 is an important ATP-binding cassette transporter impacting the absorption and distribution of over 200 chemical toxins and drugs. ABCG2 also reduces the cellular accumulation of diverse chemotherapeutic agents. Acquired somatic mutations in the phylogenetically conserved amino acids of ABCG2 might provide unique insights into its molecular mechanisms of transport. Here, we identify a tumor-derived somatic mutation (Q393K) that occurs in a highly conserved amino acid across mammalian species. This ABCG2 mutant seems incapable of providing ABCG2-mediated drug resistance. This was perplexing because it is localized properly and retained interaction with substrates and nucleotides. Using a conformationally sensitive antibody, we show that this mutant appears “locked” in a non-functional conformation. Structural modeling and molecular dynamics simulations based on ABCG2 cryo-EM structures suggested that the Q393K interacts with the E446 to create a strong salt bridge. The salt bridge is proposed to stabilize the inward-facing conformation, resulting in an impaired transporter that lacks the flexibility to readily change conformation, thereby disrupting the necessary communication between substrate binding and transport.