JED: a Java Essential Dynamics Program for comparative analysis of protein trajectories.

BACKGROUND: Essential Dynamics (ED) is a common application of principal component analysis (PCA) to extract biologically relevant motions from atomic trajectories of proteins. Covariance and correlation based PCA are two common approaches to determine PCA modes (eigenvectors) and their eigenvalues. Protein dynamics can be characterized in terms of Cartesian coordinates or internal distance pairs. In understanding protein dynamics, a comparison of trajectories taken from a set of proteins for similarity assessment provides insight into conserved mechanisms. Comprehensive software is needed to facilitate comparative-analysis with user-friendly features that are rooted in best practices from multivariate statistics. RESULTS: We developed a Java based Essential Dynamics toolkit called JED to compare the ED from multiple protein trajectories. Trajectories from different simulations and different proteins can be pooled for comparative studies. JED implements Cartesian-based coordinates (cPCA) and internal distance pair coordinates (dpPCA) as options to construct covariance (Q) or correlation (R) matrices. Statistical methods are implemented for treating outliers, benchmarking sampling adequacy, characterizing the precision of Q and R, and reporting partial correlations. JED output results as text files that include transformed coordinates for aligned structures, several metrics that quantify protein mobility, PCA modes with their eigenvalues, and displacement vector (DV) projections onto the top principal modes. Pymol scripts together with PDB files allow movies of individual Q- and R-cPCA modes to be visualized, and the essential dynamics occurring within user-selected time scales. Subspaces defined by the top eigenvectors are compared using several statistical metrics to quantify similarity/overlap of high dimensional vector spaces. Free energy landscapes can be generated for both cPCA and dpPCA. CONCLUSIONS: JED offers a convenient toolkit that encourages best practices in applying multivariate statistics methods to perform comparative studies of essential dynamics over multiple proteins. For each protein, Cartesian coordinates or internal distance pairs can be employed over the entire structure or user-selected parts to quantify similarity/differences in mobility and correlations in dynamics to develop insight into protein structure/function relationships.

Investigation on interaction of tannic acid with type I collagen and its effect on thermal, enzymatic, and conformational stability for tissue engineering applications.

Collagen is an essential component of tissues, which is the most abundant component in extracellular matrix and highly conserved across the animal kingdom. It can assemble into fiber and play an essential role in cell adhesion and growth and could be extremely useful in tissue engineering. In this study, the effect of tannic acid (TA) on the thermal, enzymatic and conformational stability of type I collagen has been investigated for the development of collagen-based biomaterials. Interaction of TA with collagen demonstrates the role of hydrogen bonding and hydrophobic interaction in providing the thermal and enzymatic stability. Thermal analysis studies reveal that, hydrothermal stability of collagen increases as well as inhibits the breakdown of collagenase by formation of hydrogen bonds and hydrophobic interactions. TA binds to the collagen with high affinity because the structural flexibility of the collagen compensates for the structural rigidity of the phenolics. Increase in concentration of TA induces significant change in the conformation of triple helix. The free binding energy of TA with collagen-like peptide was determined to be in the range of -9.4 to -11.2 kcal mol(-1), which was calculated by using Autodock Vina software and showed numerous hydrophobic and hydrogen bond interactions. We anticipate that these collagen-based biomaterials hold great potential for biomedical applications.

In Vitro characterization of the endocrine disrupting effects of per- and poly-fluoroalkyl substances (PFASs) on the human androgen receptor.

Per- and poly-fluoroalkyl substances (PFASs) are used extensively in a broad range of industrial applications and consumer products. While a few legacy PFASs have been voluntarily phased out, over 5000 PFASs have been produced as replacements for their predecessors. The potential endocrine disrupting hazards of most emerging PFASs have not been comprehensively investigated. In silico molecular docking to the human androgen receptor (hAR) combined with machine learning techniques were previously applied to 5206 PFASs and predicted 23 PFASs bind the hAR. Herein, the in silico results were validated in vitro for the five candidate AR ligands that were commercially available. Three manufactured PFASs namely (9-(nonafluorobutyl)- 2,3,6,7-tetrahydro-1 H,5 H,11 H-pyrano[2,3-f]pyrido[3,2,1-ij]quinolin-11-one (NON), 2-(heptafluoropropyl)- 3-phenylquinoxaline (HEP), and 2,2,3,3,4,4,5,5,5-nonafluoro-N-(4-nitrophenyl)pentanamide (NNN) elicited significant antiandrogenic effects at relatively low concentrations. We further investigated the mechanism of AR inhibition and found that all three PFASs inhibited AR transactivation induced by testosterone through a competitive binding mechanism. We then examined the antiandrogenic effects of these PFASs on AR expression and its responsive genes. Consistently, these PFASs significantly decreased the expression of PSA and FKBP5 and increased the expression of AR, similar to the effects elicited by a known competitive AR inhibitor, hydroxyflutamide. This suggests they are competitive antagonists of AR activity and western blot analysis revealed these PFASs decreased intracellular AR protein in androgen sensitive human prostate cancer cells. Hence, the findings presented here corroborate our published in silico approach and indicate these emerging PFASs may adversely affect the human endocrine system.

Adsorption of collagen onto single walled carbon nanotubes: a molecular dynamics investigation.

Classical molecular dynamics (MD) simulation has been carried out to understand the adsorption of collagen like peptides onto single walled carbon nanotubes (CNT) in an aqueous environment. It is observed that the triple helical structure of all the model collagen like peptides (CPs) has been unaltered upon adsorption onto CNT. The model CPs do not wrap around the CNT, however, the axis of the triple helix subtends a cross angle with respect to the axis of the CNT. The interaction between the CPs and CNT as well as that between the CPs and water molecules was observed by MD simulation snapshots. The inherent nature of the interaction of CPs with CNT facilitates the penetration of CPs into the water/CNT interface. During this process, water molecules trapped between the CPs and CNT are appreciably displaced. Although, hydrophobic-hydrophobic interaction is crucial for the interaction, the role of πR (R = OH and NH(2)) interactions are also observed from the geometrical parameters. The sequence specific interaction of CPs with CNT is evident from the results. It is found that the length of the CNT, curvature of the CNT and length of the CPs do not significantly influence interaction between the two systems. Overall the findings provide important information for the development of nanocomposite materials from collagen and CNT.

Identification of an Electrostatic Ruler Motif for Sequence-Specific Binding of Collagenase to Collagen.

Sequence-specific cleavage of collagen by mammalian collagenase plays a pivotal role in cell function. Collagenases are matrix metalloproteinases that cleave the peptide bond at a specific position on fibrillar collagen. The collagenase Hemopexin-like (HPX) domain has been proposed to be responsible for substrate recognition, but the mechanism by which collagenases identify the cleavage site on fibrillar collagen is not clearly understood. In this study, Brownian dynamics simulations coupled with atomic-detail and coarse-grained molecular dynamics simulations were performed to dock matrix metalloproteinase-1 (MMP-1) on a collagen IIIα1 triple helical peptide. We find that the HPX domain recognizes the collagen triple helix at a conserved R-X11-R motif C-terminal to the cleavage site to which the HPX domain of collagen is guided electrostatically. The binding of the HPX domain between the two arginine residues is energetically stabilized by hydrophobic contacts with collagen. From the simulations and analysis of the sequences and structural flexibility of collagen and collagenase, a mechanistic scheme by which MMP-1 can recognize and bind collagen for proteolysis is proposed.

Molecular dynamic simulation studies on the effect of one residue chain staggering on the structure and stability of heterotrimeric collagen-like peptides with interruption.

A systematic molecular dynamics (MD) simulation has been carried out on collagen-like peptides with different combinations of interruptions in the Gly-X(AA) -Y(AA) repeats. Although experimental studies have been carried out to elucidate the structural consequences of homotrimeric collagen-like peptides, this is the first report on the structural effect on the heterotrimeric models with G4G and G1G breaks present simultaneously in the constituent chains with difference in one residue chain staggering. The results reveal that the axial registry of the interrupted region changes significantly from that of conventional triple helical peptide without interruption. Further, results from MD simulations show the formation of a kink in the interrupted region of the triple-helical peptides. The conformational analysis reveals that the interruption in the Gly-X(AA) -Y(AA) pattern in these peptides induces ß-strand conformation in triple helical peptides. The conventional hydrogen bonds in the interrupted triad are affected and new nonconventional H-bonds are formed in the triple helical structure, and as a result interrupted region becomes locally fragile. MM-PBSA calculations on the different systems clearly suggest that the binding affinity varies marginally due to one residue staggering. However, it is found from the structural parameters that hydrogen-bonding pattern differs significantly due to the difference in the staggering of chains.

3D-QSAR and docking studies on the HEPT derivatives of HIV-1 reverse transcriptase.

Three-dimensional Quantitative Structure Activity Relationship (3D-QSAR) has been derived for a set of HEPT derivatives of HIV-1 reverse transcriptase (RT) using Comparative Molecular Field Analysis (CoMFA). The CoMFA models have been developed using two different alignment procedures such as common substructure and bioactive conformation. The CoMFA model I is derived from a common substructure procedure that includes steric and electrostatic fields with the cross-validated q(2) and the non-cross-validated r(2) value of 0.86 and 0.97, respectively. The same for the CoMFA model II that is derived based on the bioactive conformation are 0.19 and 0.77, respectively. It is evident from the results that the common substructure-based alignment model has good statistical significance when compared with that of bioactive conformation for the selected systems in this study. The docking study revealed that the conformational flexibility observed at the R3 position favors different orientations of the substitution at the active site of HIV-1 RT and thereby leads to inconsistency in the CoMFA alignment based on bioactive conformation.