Structures prediction and replica exchange molecular dynamics simulations of α-synuclein: A case study for intrinsically disordered proteins.

In recent years, a variety of three-dimensional structure prediction tools, including AlphaFold2, AlphaFold3, I-TASSER, C-I-TASSER, Phyre2, ESMFold, and RoseTTAFold, have been employed in the investigation of intrinsically disordered proteins. However, a comprehensive validation of these tools specifically for intrinsically disordered proteins has yet to be conducted. In this study, we utilize AlphaFold2, AlphaFold3, I-TASSER, C-I-TASSER, Phyre2, ESMFold, and RoseTTAFold to predict the structure of a model intrinsically disordered α-synuclein protein. Additionally, extensive replica exchange molecular dynamics simulations of the intrinsically disordered protein are conducted. The resulting structures from both structure prediction tools and replica exchange molecular dynamics simulations are analyzed for radius of gyration, secondary and tertiary structure properties, as well as Cα and Hα chemical shift values. A comparison of the obtained results with experimental data reveals that replica exchange molecular dynamics simulations provide results in excellent agreement with experimental observations. However, none of the structure prediction tools utilized in this study can fully capture the structural characteristics of the model intrinsically disordered protein. This study shows that a cluster of ensembles are required for intrinsically disordered proteins. Artificial-intelligence based structure prediction tools such as AlphaFold3 and C-I-TASSER could benefit from stochastic sampling or Monte Carlo simulations for generating an ensemble of structures for intrinsically disordered proteins.

Bioinformatics approach for prediction and analysis of the Non-Structural Protein 4B (NSP4B) of the Zika virus.

BACKGROUND: The Nonstructural Protein (NSP) 4B of Zika virus of 251 amino acids from (ZIKV/Human/POLG_ZIKVF) with accession number (A0A024B7W1), Induces the production of Endoplasmic Reticulum ER-derived membrane vesicles, which are the sites of viral replication. To understand the physical basis of how proteins fold in nature and to solve the challenge of protein structure prediction, Ab-initio and comparative modeling are crucial tools. RESULTS: The systematic in silico technique, ThreaDom, had only predicted one domain (4 - 190) of NSP4B. I-TASSER, and Alphafold were ranked as the best servers for full-length 3-D protein structure predictions of NSP4B, where the predicted models were evaluated quantitatively using benchmarked metrics including C-score (-3.43), TM-score (0.77949), RMSD (2.73), and Z-score (1.561). The functional and protein binding motifs were realized using motif databases, secondary and surface accessibility predictions combined with Post-Translational Modification Sites (PTMs) prediction. Two highly conserved protein-binding motifs (Flavi NS4B and Bacillus papRprotein), together with three (PTMs) (Casein Kinase II, Myristyl site, and ASN-Glycosylation site) were predicted utilizing the Motif scan and Scanprosite servers. These patterns and PTMs were associated with NSP4B's role in triggering the development of the viral replication complex and its participation in the localization of NS3 and NS5 on the membrane. Only one hit from Structural Classification of Protein (SCOP) matched the protein sequence at positions 10 to 397 and was categorized six-hairpin glycosidases superfamily according to CATH (Class, Architecture, Topology, and Homology). Integrating this NSP4B information with the templates' SCOP and CATH annotations achieves it easier to attribute structure-function/evolution links to both previously known and recently discovered protein structures.

Comparison, Analysis, and Molecular Dynamics Simulations of Structures of a Viral Protein Modeled Using Various Computational Tools.

The structural analysis of proteins is a major domain of biomedical research. Such analysis requires resolved three-dimensional structures of proteins. Advancements in computer technology have led to progress in biomedical research. In silico prediction and modeling approaches have facilitated the construction of protein structures, with or without structural templates. In this study, we used three neural network-based de novo modeling approaches-AlphaFold2 (AF2), Robetta-RoseTTAFold (Robetta), and transform-restrained Rosetta (trRosetta)-and two template-based tools-the Molecular Operating Environment (MOE) and iterative threading assembly refinement (I-TASSER)-to construct the structure of a viral capsid protein, hepatitis C virus core protein (HCVcp), whose structure have not been fully resolved by laboratory techniques. Templates with sufficient sequence identity for the homology modeling of complete HCVcp are currently unavailable. Therefore, we performed domain-based homology modeling for MOE simulations. The templates for each domain were obtained through sequence-based searches on NCBI and the Protein Data Bank. Then, the modeled domains were assembled to construct the complete structure of HCVcp. The full-length structure and two truncated forms modeled using various computational tools were compared. Molecular dynamics (MD) simulations were performed to refine the structures. The root mean square deviation of backbone atoms, root mean square fluctuation of Cα atoms, and radius of gyration were calculated to monitor structural changes and convergence in the simulations. The model quality was evaluated through ERRAT and phi-psi plot analysis. In terms of the initial prediction for protein modeling, Robetta and trRosetta outperformed AF2. Regarding template-based tools, MOE outperformed I-TASSER. MD simulations resulted in compactly folded protein structures, which were of good quality and theoretically accurate. Thus, the predicted structures of certain proteins must be refined to obtain reliable structural models. MD simulation is a promising tool for this purpose.

A Case Study on PPM1D and 9 Other Shared Germline Alterations in a Family.

BACKGROUND: The use of high-throughput genotyping techniques has enabled us to identify the rare germline genetic variants with different pathogenicity and penetrance, and understand their role in cancer predisposition. We report here a familial cancer case, a study from Western Indian. METHODS: NGS-WES was carried out in a lung cancer patient who has a family history of multiple cancers across generations, including tongue, lung, brain, cervical, urothelial, and esophageal cancer. The results were validated by data mining from available data bases. I-TASSER, RasMol and PyMol were used for protein structure modelling. RESULTS: The sequencing by NGS-WES revealed PPM1D c.1654C>T (p.Arg552Ter) mutation in hotspot region exon 6 leading to sudden protein truncation and loss of the C-terminal, due to the substitution of C>T. This mutation was classified as a variant of uncertain significance (VUS), due to limited data on lung cancer, The three unaffected siblings of proband did not show any pathogenic variants and comparative analysis of the four siblings indicate 9 shared genetic variants, classified as benign as per ClinVar. CONCLUSION: PPM1D constitutional genetic alterations are rare and uncommon in different ethnic populations. This gene encodes a phosphatase playing role in regulating the P53 tumor suppressor pathway and DNA damage response. Genetic alterations in the PPM1D gene maybe linked to history of gliomas, breast cancer, and ovarian cancer onset in the proband's family.
.

Molecular modeling, simulation and docking of Rv1250 protein from Mycobacterium tuberculosis.

Computational prediction and protein structure modeling have come to the aid of various biological problems in determining the structure of proteins. These technologies have revolutionized the biological world of research, allowing scientists and researchers to gain insights into their biological questions and design experimental research much more efficiently. Pathogenic Mycobacterium spp. is known to stay alive within the macrophages of its host. Mycobacterium tuberculosis is an acid-fast bacterium that is the most common cause of tuberculosis and is considered to be the main cause of resistance of tuberculosis as a leading health issue. The genome of Mycobacterium tuberculosis contains more than 4,000 genes, of which the majority are of unknown function. An attempt has been made to computationally model and dock one of its proteins, Rv1250 (MTV006.22), which is considered as an apparent drug-transporter, integral membrane protein, and member of major facilitator superfamily (MFS). The most widely used techniques, i.e., homology modeling, molecular docking, and molecular dynamics (MD) simulation in the field of structural bioinformatics, have been used in the present work to study the behavior of Rv1250 protein from M. tuberculosis. The structure of unknown TB protein, i.e., Rv1250 was retrived using homology modeling with the help of I-TASSER server. Further, one of the sites responsible for infection was identified and docking was done by using the specific Isoniazid ligand which is an inhibitor of this protein. Finally, the stability of protein model and analysis of stable and static interaction between protein and ligand molecular dynamic simulation was performed at 100 ns The designing of novel Rv1250 enzyme inhibitors is likely achievable with the use of proposed predicted model, which could be helpful in preventing the pathogenesis caused by M. tuberculosis. Finally, the MD simulation was done to evaluate the stability of the ligand for the specific protein.

Comparative In Silico Analysis and Functional Characterization of TANK-Binding Kinase 1-Binding Protein 1.

Protein modelling plays a vital role in the drug discovery process. TANK-binding kinase 1-binding protein 1 is also called an adapter protein, which is encoded by gene TBK1 present in Homo sapiens. It is found in lungs, small intestine, leukocytes, heart, placenta, muscle, kidney, lower level of thymus, and brain. It has a number of protein-binding sites, to which TBK1 and IKBKE bind and perform different functions as immunomodulatory, antiproliferative, and antiviral innate immunity which release different types of interferons. Our study predicts the comparative model of 3-dimensional (3D) structure through different bioinformatics tools that will be helpful for further studies in future. The reactivity and stability of these proteins were evaluated physicochemically and through domain determination and prediction of secondary structure using bioinformatics methods such as ProtParam, Pfam, and SOPMA, respectively. Robetta, an ab initio approach, I-TASSER, and AlphaFold was used for 3D structure prediction, and the models were validated using the SAVESv6.0 (PROCHECK) server. Conclusively, the best 3D structure of TBK1-binding protein 1 was predicted using Robetta software. After unveiling the 3D structure of the novel protein, we concluded that this structure will help us to find out its role other than in antiviral innate immunity and by producing torsion in its 3D structure researchers will be able to detect either this protein is involved in any disease or not because according to previous studies it was not associated with any disease.

Simulated Docking Predicts Putative Channels for the Transport of Long-Chain Fatty Acids in Vibrio cholerae.

Fatty acids (FA) play an important role in biological functions, such as membrane homeostasis, metabolism, and as signaling molecules. FadL is the only known protein that uptakes long-chain fatty acids in Gram-negative bacteria, and this uptake has traditionally been thought to be limited to fatty acids up to 18 carbon atoms in length. Recently however, it was found Vibrio cholerae has the ability to uptake fatty acids greater than 18 carbon atoms and this uptake corresponds to bacterial survivability. Using E. coli's FadL as a template, V. cholerae FadL homologs vc1042, vc1043, and vca0862 have been computationally folded, simulated on an atomistic level using Molecular Dynamics, and docked in silico to analyze the FadL transport channels. For the vc1042 and vc1043 homologs, these transport channels have more structural accommodations for the many rigid unsaturated bonds of long-chain polyunsaturated fatty acids, while the vca0862 homolog was found to lack transport channels within the signature beta barrel of FadL proteins.

In-silico analysis of nonsynonymous genomic variants within CCM2 gene reaffirm the existence of dual cores within typical PTB domain.

PURPOSE: The objective of this study is to validate the existence of dual cores within the typical phosphotyrosine binding (PTB) domain and to identify potentially damaging and pathogenic nonsynonymous coding single nuclear polymorphisms (nsSNPs) in the canonical PTB domain of the CCM2 gene that causes cerebral cavernous malformations (CCMs). METHODS: The nsSNPs within the coding sequence for PTB domain of human CCM2 gene, retrieved from exclusive database searches, were analyzed for their functional and structural impact using a series of bioinformatic tools. The effects of mutations on the tertiary structure of the PTB domain in human CCM2 protein were predicted to examine the effect of nsSNPs on the tertiary structure of PTB Cores. RESULTS: Our mutation analysis, through alignment of protein structures between wildtype CCM2 and mutant, predicted that the structural impacts of pathogenic nsSNPs is biophysically limited to only the spatially adjacent substituted amino acid site with minimal structural influence on the adjacent core of the PTB domain, suggesting both cores are independently functional and essential for proper CCM2 PTB function. CONCLUSION: Utilizing a combination of protein conservation and structure-based analysis, we analyzed the structural effects of inherited pathogenic mutations within the CCM2 PTB domain. Our results predicted that the pathogenic amino acid substitutions lead to only subtle changes locally, confined to the surrounding tertiary structure of the PTB core within which it resides, while no structural disturbance to the neighboring PTB core was observed, reaffirming the presence of independently functional dual cores in the CCM2 typical PTB domain.

In silico analysis of AhyI protein and AI-1 inhibition using N-cis-octadec-9z-enoyl-l-homoserine lactone inhibitor in Aeromonas hydrophila.

AhyI is homologous to the protein LuxI and is conserved throughout bacterial species including Aeromonas hydrophila. A. hydrophila causes opportunistic infections in fish and other aquatic organisms. Furthermore, this pathogennot only poses a great risk for the aquaculture industry, but also for human public health. AhyI (expressing acylhomoserine lactone) is responsible for the biosynthesis of autoinducer-1 (AI-1), commonly referred to as a quorum sensing (QS) signaling molecule, which plays an essential role in bacterial communication. Studying protein structure is essential for understanding molecular mechanisms of pathogenicity in microbes. Here, we have deduced a predicted structure of AhyI protein and characterized its function using in silico methods to aid the development of new treatments for controlling A.hydrophila infections. In addition to modeling AhyI, an appropriate inhibitor molecule was identified via high throughput virtual screening (HTVS) using mcule drug-like databases.The AhyI-inhibitor N-cis-octadec-9Z-enoyl-l-Homoserine lactone was selected withthe best drug score. In order to understand the pocket sites (ligand binding sites) and their interaction with the selected inhibitor, docking (predicted protein binding complex) servers were used and the selected ligand was docked with the predicted AhyI protein model. Remarkably, N-cis-octadec-9Z-enoyl-l-Homoserine lactone established interfaces with the protein via16 residues (V24, R27, F28, R31, W34, V36, D45, M77, F82, T101, R102, L103, 104, V143, S145, and V168), which are involved with regulating mechanisms of inhibition. These proposed predictions suggest that this inhibitor molecule may be used as a novel drug candidate for the inhibition of auto-inducer-1 (AI-1) activity.The N-cis-octadec-9Z-enoyl-l-Homoserine lactone inhibitor molecule was studied on cultured bacteria to validate its potency against AI-1 production. At a concentration of 40 µM, optimal inhibition efficiency of AI-1 was observedin bacterial culture media.These results suggest that the inhibitor molecule N-cis-octadec-9Z-enoyl-l-Homoserine lactone is a competitive inhibitor of AI-1 biosynthesis.

Prokaryotic expression, evaluation, and prediction of the structure and function of the ecarin metalloproteinase domain.

Ecarin is one of the most widely used drug compounds in blood clotting experiments and is used to monitor and treat many diseases such as cancer, liver, lupus, and cardiovascular disease. The metalloproteinase domain is known as the active site of ecarin. In this study, an ecarin metalloproteinase cassette was designed and synthesized in the pUC57 vector. The gene fragment was released and cloned into the pET-28a vector and expressed in Escherichia coli. The recombinant protein was confirmed by western blotting. Enzyme activity was estimated by a laboratory coagulation test, and prothrombin time and tertiary structure were determined by using the Iterative Threading ASSEmbly Refinement (I-TASSER) server. Data from blood clotting tests for the produced ecarin activity were analyzed using an independent t test. As per I-TASSER server prediction, model 1 with the highest confidence score 0.95, template modeling score (0.84 ± 0.08), and root mean square deviation (3.5 ± 2.4 Å) was considered as the best model, and the 2e3xA enzyme was more similar to the target protein. The predictive results helped to better understand the relationship between the structure and function of the ecarin metalloproteinase domain. Also, the production of this active site in the prokaryotic expression system, which is simpler and more cost-effective than the production of the eukaryotic system, showed that this recombinant ecarin could be used as a substitute for the raw snake venom of Echis carinatus because it converts prothrombin into thrombin, and its activity, as estimated using the prothrombin time test, was found to be faster than normal ecarin.

Efficient production of oxidized terpenoids via engineering fusion proteins of terpene synthase and cytochrome P450.

The functionalization of terpenes using cytochrome P450 enzymes is a versatile route to the production of useful derivatives that can be further converted to value-added products. Many terpenes are hydrophobic and volatile making their availability as a substrate for P450 enzymes significantly limited during microbial production. In this study, we developed a strategy to improve the accessibility of terpene molecules for the P450 reaction by linking terpene synthase and P450 together. As a model system, fusion proteins of 1,8-cineole synthase (CS) and P450cin were investigated and it showed an improved hydroxylation of the monoterpenoid 1,8-cineole up to 5.4-fold. Structural analysis of the CS-P450cin fusion proteins by SEC-SAXS indicated a dimer formation with preferred orientations of the active sites of the two domains. We also applied the enzyme fusion strategy to the oxidation of a sesquiterpene epi-isozizaene and the fusion enzymes significantly improved albaflavenol production in engineered E. coli. From the analysis of positive and negative examples of the fusion strategy, we proposed key factors in structure-based prediction and evaluation of fusion enzymes. Developing fusion enzymes for terpene synthase and P450 presents an efficient strategy toward oxidation of hydrophobic terpene compounds. This strategy could be widely applicable to improve the biosynthetic titer of the functionalized products from hydrophobic terpene intermediates.

A Computational Approach for Identifying Potential Phytochemicals Against Non-structural Protein 1 (Nsp1) of SARS-CoV-2.

AIM AND OBJECTIVE: A recent study has revealed that non-structural protein 1 (Nsp1) of the SARS-CoV-2 is one of the novel targets for developing new antiviral drugs. To date, there is no significant exact medication available to treat Covid-19. As a result, both the death toll and the number of people affecting by this disease are increasing with each passing day. 35 phytochemicals having antiviral properties were taken to get the best compounds against Nsp1 Materials and Methods: As no PDB structure of this protein is available, homology modeling was done to predict the probable structure. After homology modeling, the best model was taken according to C-score and TM- score and then validated using different web servers. After validation, docking of these compounds was done using AutoDock vina, vega zz, and PyRx, and consensus docking score was considered to select molecules after docking. Finally, the orbitals energy calculation of these compounds was done to check their activity and the binding interactions of these molecules also analyzed. RESULTS: Molecules having a consensus score of -8kcal/mol or more negative were kept for further study and it was seen that 16 molecules had the given criteria. Then, drug-likeness filtration was done according to Lipinski's rule of five and 11 molecules remained. Out of these 11 molecules, 5 molecules had satisfactory ADMET properties. Calculation of orbital energy revealed their activity. CONCLUSION: It is expected that this research might be helpful for the development of new antiviral drugs active against SARS-CoV-2 targeting Nsp1.

Identification of Yak's TLR4 Alternative Spliceosomes and Bioinformatic Analysis of TLR4 Protein Structure and Function.

In this study, the yak's TLR4 gene alternative spliceosomes were investigated using PCR amplification and cloning to improve disease-resistance in yak and promote efficient utilization of yak's resources. qRT-PCR was used to determine the expression levels of two alternatively spliced transcripts of the TLR4 gene in seven distinct tissues. To predict the function of proteins expressed by each TLR4 spliceosome, bioinformatic analysis of yak's TLR4 protein structure and function was performed, which led to the identification of two alternative spliceosomes of yak's TLR4 gene. The TLR4-X1 sequence length was 2526 bp, and it encoded full-length TLR4 protein (841 amino acids). The sequence length of the exon-2 deleted TLR4-X2 sequence was 1926 bp, and it encoded truncated TLR4 protein (641 amino acids). TLR4-X2 sequence was consistent with the predicted sequence of the TLR4 gene in GenBank. Each tissue showed significantly different expression levels of these two alternative spliceosomes. As per the bioinformatic analysis of the structure and function of TLR4 protein, deletion of exon-2 in the TLR4 gene resulted in frameshift mutations of the reading frame in the corresponding protein, which altered its ligand-binding and active sites. Besides, biological property such as substrate specificity of truncated TLR4 protein was also altered, leading to altered protein function. This study has laid a theoretical foundation for exploring the role of two variants of the TLR4 gene in yak's disease resistance. Besides, this study's data could be analyzed further to explore the molecular mechanism associated with disease-resistance in the yak.

Bioinformatic Prediction of Gene Ontology Terms of Uncharacterized Proteins from Chromosome 11.

In chromosome 11, 71 out of its 1254 proteins remain functionally uncharacterized on the basis of their existence evidence (uPE1s) following the latest version of neXtProt (release 2020-01-17). Because in vivo and in vitro experimental strategies are often time-consuming and labor-intensive, there is a need for a bioinformatics tool to predict the function annotation. Here, we used I-TASSER/COFACTOR provided on the neXtProt web site, which predicts gene ontology (GO) terms based on the 3D structure of the protein. I-TASSER/COFACTOR predicted 2413 GO terms with a benchmark dataset of the 22 proteins belonging to PE1 of chromosome 11. In this study, we developed a filtering algorithm in order to select specific GO terms using the GO map generated by I-TASSER/COFACTOR. As a result, 187 specific GO terms showed a higher average precision-recall score at the least cellular component term compared to 2413 predicted GO terms. Next, we applied 65 proteins belonging to uPE1s of chromosome 11, and then 409 out of 6684 GO terms survived, where 103 and 142 GO terms of molecular function and biological process, respectively, were included. Representatively, the cellular component GO terms of CCDC90B, C11orf52, and the SMAP were predicted and validated using the overexpression system into 293T cells and immunofluorescence staining. We will further study their biological and molecular functions toward the goal of the neXt-CP50 project as a part of C-HPP. We shared all results and programs in Github (https://github.com/heeyounh/I-TASSER-COFACTOR-filtering.git).

Structural modeling, mutation analysis, and in vitro expression of usherin, a major protein in inherited retinal degeneration and hearing loss.

Usherin is the most common causative protein associated with autosomal recessive retinitis pigmentosa (RP) and Usher syndrome (USH), which are characterized by retinal degeneration alone and in combination with hearing loss, respectively. Usherin is essential for photoreceptor survival and hair cell bundle integrity. However, the molecular mechanism underlying usherin function in normal and disease conditions is unclear. In this study, we investigated structural models of usherin domains and localization of usherin pathogenic small in-frame mutations, mainly homozygous missense mutations. We found that usherin fibronectin III (FN3) domains and most laminin-related domains have a ß-sandwich structure. Some FN3 domains are predicted to interact with each other and with laminin-related domains. The usherin protein may bend at some FN3 linker regions. RP- and USH-associated small in-frame mutations are differentially located in usherin domains. Most of them are located at the periphery of ß-sandwiches, with some at the interface between interacting domains. The usherin laminin epidermal growth factor repeats adopt a rod-shaped structure, which is maintained by disulfide bonds. Most missense mutations and deletion of exon 13 in this region disrupt the disulfide bonds and may affect local protein folding. Despite low expression of the recombinant entire protein and protein fragments in mammalian cell culture, usherin FN3 fragments are more robustly expressed and secreted than its laminin-related fragments. Our findings provide new insights into the usherin structure and the disease mechanisms caused by pathogenic small in-frame mutations, which will help inform future experimental research on diagnosis, disease mechanisms, and therapeutic approaches.

Structural and posttranslational analysis of human calcium-binding protein, spermatid-associated 1.

Recently, we detected a novel biomarker in human saliva called calcium-binding protein, spermatid-associated 1 (CABS1). CABS1 protein had previously been described only in testis, and little was known of its characteristics other than it was considered a structurally disordered protein. Levels of human CABS1 (hCABS1) in saliva correlate with stress, whereas smaller sized forms of hCABS1 in saliva are associated with resilience to stress. Interestingly, hCABS1 also has an anti-inflammatory peptide sequence near its carboxyl terminus, similar to that of a rat prohormone, submandibular rat 1. We performed phylogenetic and sequence analysis of hCABS1. We found that from 72 CABS1 sequences currently annotated in the National Center for Biotechnology Information protein database, only 14 contain the anti-inflammatory domain "TxIFELL," all of which are primates. We performed structural unfoldability analysis using PONDER and FoldIndex and discovered three domains that are highly disordered. Predictions of three-dimensional structure of hCABS1 using RaptorX, IonCom, and I-TASSER software agreed with these findings. Predicted neutrophil elastase cleavage density also correlated with hCABS1 regions of high structural disorder. Ligand binding prediction identified Ca2+ , Mg2+ , Zn2+ , leucine, and thiamine pyrophosphate, a pattern observed in enzymes associated with energy metabolism and mitochondrial localization. These new observations on hCABS1 raise intriguing questions about the interconnection between the autonomic nervous system, stress, and the immune system. However, the precise molecular mechanisms involved in the complex biology of hCABS1 remain unclear. We provide a detailed in silico analysis of relevant aspects of the structure and function of hCABS1 and postulate extracellular and intracellular roles.

Molecular Cloning of an Amino Acid Permease Gene and Structural Characterization of the Protein in Common Bean (Phaseolus vulgaris L.).

Plants synthesize amino acids by collateral metabolic pathways using primary elements carbon and oxygen from air, hydrogen from water in soil and nitrogen from soil. Following synthesis, amino acids are immediately used for metabolism, transient storage or transported to the phloem. Different families of transporters have been identified for import of amino acids into plant cells. The first identified amino acid transporter, amino acid permease 1 (AAP1) in Arabidopsis belongs to a family of eight members and transports acidic, neutral, and basic amino acids. Legumes fix atmospheric nitrogen through a symbiotic relationship with root nodules bacteria. Following fixation, nitrogen is reduced to amino acids and is exported via different amino acid transporters. However, information is lacking about the structure of these important classes of amino acid transporter proteins in plant. We have amplified AAP from Phaseolus vulgaris, an economically important leguminous plant grown all over the world, and sequenced. The sequence has been characterized in silico and a three-dimensional structure of AAP has been predicted and validated. The information obtained not only enhances the knowledge about the structure of an amino acid permease gene in P. vulgaris, but will also help in designing protein-ligand studies using this protein as well.

Molecular Dynamics Simulation Reveals Exposed Residues in the Ligand-Binding Domain of the Low-Density Lipoprotein Receptor that Interacts with Vesicular Stomatitis Virus-G Envelope.

Familial hypercholesterolemia (FH) is an autosomal dominant disease most often caused by mutations in the low-density lipoprotein receptor (LDLR) gene, which consists of 18 exons spanning 45 kb and codes for a precursor protein of 860 amino acids. Mutations in the LDLR gene lead to a reduced hepatic clearance of LDL as well as a high risk of coronary artery disease (CAD) and sudden cardiac death (SCD). Recently, LDLR transgenes have generated interest as potential therapeutic agents. However, LDLR packaging using a lentiviral vector (LVV) system pseudotyped with a vesicular stomatitis virus (VSV)-G envelope is not efficient. In this study, we modified the LVV system to improve transduction efficiency and investigated the LDLR regions responsible for transduction inhibition. Transduction efficiency of 293T cells with a 5'-LDLReGFP-3' fusion construct was only 1.55% compared to 42.32% for the eGFP construct. Moreover, co-expression of LDLR affected eGFP packaging. To determine the specific region of the LDLR protein responsible for packaging inhibition, we designed constructs with mutations or sequential deletions at the 3' and 5' ends of LDLR cDNA. All constructs except one without the ligand-binding domain (LBD) (pWoLBD-eGFP) resulted in low transduction efficiency, despite successful packaging of viral RNA in the VSV envelope, as confirmed through RT-PCR. When we evaluated a direct interaction between LDLR and the VSV envelope glycoprotein using MD simulation and protein-protein interactions, we uncovered Val119, Thr120, Thr67, and Thr118 as exposed residues in the LDLR receptor that interact with the VSV protein. Together, our results suggest that the LBD of LDLR interacts with the VSV-G protein during viral packaging, which significantly reduces transduction efficiency.

Blinded Testing of Function Annotation for uPE1 Proteins by I-TASSER/COFACTOR Pipeline Using the 2018-2019 Additions to neXtProt and the CAFA3 Challenge.

In 2018, we reported a hybrid pipeline that predicts protein structures with I-TASSER and function with COFACTOR. I-TASSER/COFACTOR achieved Gene Ontology (GO) high prediction accuracies of Fmax = 0.69 and 0.57 for molecular function (MF) and biological process (BP), respectively, on 100 comprehensively annotated proteins. Now we report blinded analyses of newly annotated proteins in the critical assessment of function annotation (CAFA) three function prediction challenge and in neXtProt. For CAFA3 results released in May 2019, our predictions on 267 and 912 human proteins with newly annotated MF and BP terms achieved Fmax = 0.50 and 0.42, respectively, on "No Knowledge" proteins, and 0.51 and 0.74, respectively, on "Limited Knowledge" proteins. While COFACTOR consistently outperforms simple homology-based analysis, its accuracy still depends on template availability. Meanwhile, in neXtProt 2019-01, 25 proteins acquired new function annotation through literature curation at UniProt/Swiss-Prot. Before the release of these curated results, we submitted to neXtProt blinded predictions of free-text function annotation based on predicted GO terms. For 10 of the 25, a good match of free-text or GO term annotation was obtained. These blind tests represent rigorous assessments of I-TASSER/COFACTOR. neXtProt now provides links to precomputed I-TASSER/COFACTOR predictions for proteins without function annotation to facilitate experimental planning on "dark proteins".

I-TASSER gateway: A protein structure and function prediction server powered by XSEDE.

There is an increasing gap between the number of known protein sequences and the number of proteins with experimentally characterized structure and function. To alleviate this issue, we have developed the I-TASSER gateway, an online server for automated and reliable protein structure and function prediction. For a given sequence, I-TASSER starts with template recognition from a known structure library, followed by full-length atomic model construction by iterative assembly simulations of the continuous structural fragments excised from the template alignments. Functional insights are then derived from comparative matching of the predicted model with a library of proteins with known function. The I-TASSER pipeline has been recently integrated with the XSEDE Gateway system to accommodate pressing demand from the user community and increasing computing costs. This report summarizes the configuration of the I-TASSER Gateway with the XSEDE-Comet supercomputer cluster, together with an overview of the I-TASSER method and milestones of its development.