Principal component analysis of alpha-helix deformations in transmembrane proteins

Immunoinformatics and Reverse Vaccinology Approach for the Identification of Potential Vaccine Candidates against Vandammella animalimors

Vandammella animalimorsus is a Gram-negative and non-motile bacterium typically transmitted to humans through direct contact with the saliva of infected animals, primarily through biting, scratches, or licks on fractured skin. The absence of a confirmed post-exposure treatment of V. animalimorsus bacterium highlights the imperative for developing an effective vaccine. We intended to determine potential vaccine candidates and paradigm a chimeric vaccine against V. animalimorsus by accessible public data analysis of the strain by utilizing reverse vaccinology. By subtractive genomics, five outer membranes were prioritized as potential vaccine candidates out of 2590 proteins. Based on the instability index and transmembrane helices, a multidrug transporter protein with locus ID A0A2A2AHJ4 was designated as a potential candidate for vaccine construct. Sixteen immunodominant epitopes were retrieved by utilizing the Immune Epitope Database. The epitope encodes the strong binding affinity, nonallergenic properties, non-toxicity, high antigenicity scores, and high solubility revealing the more appropriate vaccine construct. By utilizing appropriate linkers and adjuvants alongside a suitable adjuvant molecule, the epitopes were integrated into a chimeric vaccine to enhance immunogenicity, successfully eliciting both adaptive and innate immune responses. Moreover, the promising physicochemical features, the binding confirmation of the vaccine to the major innate immune receptor TLR-4, and molecular dynamics simulations of the designed vaccine have revealed the promising potential of the selected candidate. The integration of computational methods and omics data has demonstrated significant advantages in discovering novel vaccine targets and mitigating vaccine failure rates during clinical trials in recent years.

Designing of a multi-epitopes based vaccine against Haemophilius parainfluenzae and its validation through integrated computational approaches

Haemophilus parainfluenzae is a Gram-negative opportunist pathogen within the mucus of the nose and mouth without significant symptoms and has an ability to cause various infections ranging from ear, eye, and sinus to pneumonia. A concerning development is the increasing resistance of H. parainfluenzae to beta-lactam antibiotics, with the potential to cause dental infections or abscesses. The principal objective of this investigation is to utilize bioinformatics and immuno-informatic methodologies in the development of a candidate multi-epitope Vaccine. The investigation focuses on identifying potential epitopes for both B cells (B lymphocytes) and T cells (helper T lymphocytes and cytotoxic T lymphocytes) based on high non-toxic and non-allergenic characteristics. The selection process involves identifying human leukocyte antigen alleles demonstrating strong associations with recognized antigenic and overlapping epitopes. Notably, the chosen alleles aim to provide coverage for 90% of the global population. Multi-epitope constructs were designed by using suitable linker sequences. To enhance the immunological potential, an adjuvant sequence was incorporated using the EAAAK linker. The final vaccine construct, comprising 344 amino acids, was achieved after the addition of adjuvants and linkers. This multi-epitope Vaccine demonstrates notable antigenicity and possesses favorable physiochemical characteristics. The three-dimensional conformation underwent modeling and refinement, validated through in-silico methods. Additionally, a protein-protein molecular docking analysis was conducted to predict effective binding poses between the multi-epitope Vaccine and the Toll-like receptor 4 protein. The Molecular Dynamics (MD) investigation of the docked TLR4-vaccine complex demonstrated consistent stability over the simulation period, primarily attributed to electrostatic energy. The docked complex displayed minimal deformation and enhanced rigidity in the motion of residues during the dynamic simulation. Furthermore, codon translational optimization and computational cloning was performed to ensure the reliability and proper expression of the multi-Epitope Vaccine. It is crucial to emphasize that despite these computational validations, experimental research in the laboratory is imperative to demonstrate the immunogenicity and protective efficacy of the developed vaccine. This would involve practical assessments to ascertain the real-world effectiveness of the multi-epitope Vaccine.

Studying Conformational Properties of Transmembrane Domain of KCNE3 in a Lipid Bilayer Membrane Using Molecular Dynamics Simulations

KCNE3 is a single-pass integral membrane protein that regulates numerous voltage-gated potassium channel functions such as KCNQ1. Previous solution NMR studies suggested a moderate degree of curved α-helical structure in the transmembrane domain (TMD) of KCNE3 in lyso-myristoylphosphatidylcholine (LMPC) micelles and isotropic bicelles with the residues T71, S74 and G78 situated along the concave face of the curved helix. During the interaction of KCNE3 and KCNQ1, KCNE3 pushes its transmembrane domain against KCNQ1 to lock the voltage sensor in its depolarized conformation. A cryo-EM study of KCNE3 complexed with KCNQ1 in nanodiscs suggested a deviation of the KCNE3 structure from its independent structure in isotropic bicelles. Despite the biological significance of KCNE3 TMD, the conformational properties of KCNE3 are poorly understood. Here, all atom molecular dynamics (MD) simulations were utilized to investigate the conformational dynamics of the transmembrane domain of KCNE3 in a lipid bilayer containing a mixture of POPC and POPG lipids (3:1). Further, the effect of the interaction impairing mutations (V72A, I76A and F68A) on the conformational properties of the KCNE3 TMD in lipid bilayers was investigated. Our MD simulation results suggest that the KCNE3 TMD adopts a nearly linear α helical structural conformation in POPC-POPG lipid bilayers. Additionally, the results showed no significant change in the nearly linear α-helical conformation of KCNE3 TMD in the presence of interaction impairing mutations within the sampled time frame. The KCNE3 TMD is more stable with lower flexibility in comparison to the N-terminal and C-terminal of KCNE3 in lipid bilayers. The overall conformational flexibility of KCNE3 also varies in the presence of the interaction-impairing mutations. The MD simulation data further suggest that the membrane bilayer width is similar for wild-type KCNE3 and KCNE3 containing mutations. The Z-distance measurement data revealed that the TMD residue site A69 is close to the lipid bilayer center, and residue sites S57 and S82 are close to the surfaces of the lipid bilayer membrane for wild-type KCNE3 and KCNE3 containing interaction-impairing mutations. These results agree with earlier KCNE3 biophysical studies. The results of these MD simulations will provide complementary data to the experimental outcomes of KCNE3 to help understand its conformational dynamic properties in a more native lipid bilayer environment.

Non-canonical helical transitions and conformational switching are associated with characteristic flexibility and disorder indices in TRP and Kv channels

Structural evidence and much experimental data have demonstrated the presence of non-canonical helical substructures (π and 3₁₀) in regions of great functional relevance both in TRP as in Kv channels. Through an exhaustive compositional analysis of the sequences underlying these substructures, we find that each of them is associated with characteristic local flexibility profiles, which in turn are implicated in significant conformational rearrangements and interactions with specific ligands. We found that α-to-π helical transitions are associated with patterns of local rigidity whereas α-to-3₁₀ transitions are mainly leagued with high local flexibility profiles. We also study the relationship between flexibility and protein disorder in the transmembrane domain of these proteins. By contrasting these two parameters, we located regions showing a sort of structural discrepancy between these similar but not identical protein attributes. Notably, these regions are presumably implicated in important conformational rearrangements during the gating in those channels. In that sense, finding these regions where flexibility and disorder are not proportional allows us to detect regions with potential functional dynamism. From this point of view, we highlighted some conformational rearrangements that occur during ligand binding events, the compaction, and refolding of the outer pore loops in several TRP channels, as well as the well-known S4 motion in Kv channels.

Computational exploration of natural compounds targeting Staphylococcus aureus: inhibiting AgrA promoter binding for antimicrobial intervention

Staphylococcus aureus is a highly virulent nosocomial pathogen that poses a significant threat to individuals exposed to healthcare settings. Due to its sophisticated machinery for producing virulence factors, S. aureus can cause severe and potentially fatal infections in humans. This study focuses on the response regulator AgrA, which plays a crucial role in regulating the production of virulence factors in S. aureus. The objective is to identify natural compounds that can inhibit the binding of AgrA to its promoter site, thus inhibiting the expression of virulence genes. To achieve this, a pharmacophore model was generated using known drugs and applied to screen the ZINC natural product database. The resulting compounds were subjected to molecular docking-based virtual screening against the C-terminal DNA binding domain of AgrA. Three compounds, namely ZINC000077269178, ZINC000051012304, and ZINC000004266026, were shortlisted based on their strong affinity for key residues involved in DNA binding and transcription initiation. Subsequently, the unbound and ligand-bound complexes were subjected to a 200 ns molecular dynamics simulation to assess their conformational stability. Various analyses, including RMSD, RMSF, Rg, SASA, Principal Component Analysis, and Gibbs free energy landscape, were conducted on the simulation trajectory. The RMSD profile indicated similar fluctuations in both bound and unbound structures, while the Rg profile demonstrated the compactness of the protein without any unfolding during the simulation. Furthermore, Principal component analysis revealed that ligand binding reduced the overall atomic motion of the protein whereas free energy landscape suggested the energy variations obtained in complexes.Communicated by Ramaswamy H. Sarma.

Length-dependent collective vibrational dynamics in alpha-helices

Functions of protein molecules in nature are closely associated with their well-defined three-dimensional structures and dynamics in body fluid. So far, many efforts have been made to reveal the relation of protein structure, dynamics, and function, but the structural origin of protein dynamics, especially for secondary structures as building blocks of conformation transition, is still ambiguous. Here we theoretically uncover the collective vibrations of elastic poly-alanine α-helices and find vibration patterns that are distinctively different over residue numbers ranging from 20 to 80. Contrary to the decreasing vibration magnitude from ends to the middle region for short helices, the vibration magnitude for long helices takes the minimum at approximately 1/5 of helix length from ends but reaches a peak at the center. Further analysis indicates that major vibrational modes of helical structures strongly depend on their residue numbers, where the twist mode dominates in the vibrations of short helices while the bend mode dominates the long ones analogous to an elastic Euler beam. The helix-coil transition pathway is also affected by the alternation of the first-order mode in helices with different lengths. The dynamic properties of the helical polypeptides are promising to be harnessed for de novo design of protein-based materials and artificial biomolecules in clinical treatments.