<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"\@SimSun";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:10.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Hi Elaine,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Wow, your response is sooooooo thorough and helpful, which is exactly what I was looking for and even extended beyond that. Thank you sooooooo much for your thoughtful help! You have saved my day!!!<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Best regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Zhen<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Elaine Meng <meng@cgl.ucsf.edu><br>
<b>Date: </b>Thursday, July 21, 2022 at 15:15<br>
<b>To: </b>Gong, Zhen <zg2234@cumc.columbia.edu><br>
<b>Cc: </b>Gong, Zhen via ChimeraX-users <chimerax-users@cgl.ucsf.edu><br>
<b>Subject: </b>[EXTERNAL] Re: [chimerax-users] Match/Matchmaker multiple segments of two protein sequences<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:11.0pt">Hi Zhen,<br>
If you already have the lists of the exact residue numbers you want to pair, then you can just use the "align" command. The "matchmaker" command is to make it easy for you when you don't feel like figuring out the residue-number pairing (and the sequences
are similar enough to align).<br>
<br>
E.g. something like command:<br>
<br>
align #1:3-13,23-70,72-81,86-101@ca to #2:7-17,18-65,67-76,77-92@ca cutoffDistance 3
<br>
<br>
You have to give equal numbers of atoms, of course. So I had to change "18-66" in your example to "18-65" to work with the structures I had open, and also add "@ca" to specify only using the CA atom of each residue (since different residue types would have
different numbers of atoms).<br>
<br>
<<a href="https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Falign.html&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=b7%2BmnANKHVHEJgiq6VxvK3bwgQVmNEPB6ZUNnqB6UHc%3D&reserved=0">https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Falign.html&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=b7%2BmnANKHVHEJgiq6VxvK3bwgQVmNEPB6ZUNnqB6UHc%3D&reserved=0</a>><br>
<br>
If you were using "matchmaker" however, there are lots of options to adjust how the sequence alignment is calculated, e.g. using heaver weight of secondary structure vs. residue type and/or using a more evolutionarily distant scoring matrix. See matchmaker
alignment scoring options:<br>
<<a href="https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fmatchmaker.html%23scoring&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=7HUNI%2FR25YqGBKt3uxvap%2BMf3JTnuhU%2B0d22PGVrMu4%3D&reserved=0">https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fmatchmaker.html%23scoring&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=7HUNI%2FR25YqGBKt3uxvap%2BMf3JTnuhU%2B0d22PGVrMu4%3D&reserved=0</a>><br>
<br>
The alignment scoring options can also be adjusted in the Matchmaker GUI:<br>
<<a href="https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Ftools%2Fmatchmaker.html%23alignment&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=%2FmufuBUyyH852XOjX1DRtYwKdwx1Nn5WpnNWRFwejjY%3D&reserved=0">https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Ftools%2Fmatchmaker.html%23alignment&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=%2FmufuBUyyH852XOjX1DRtYwKdwx1Nn5WpnNWRFwejjY%3D&reserved=0</a>><br>
<br>
In Chimera (but not ChimeraX currently) you can open your own sequence alignment from a file and then use it to superimpose the structures. Chimera Multalign Viewer and how to superimpose structures after opening the alignment file:<br>
<<a href="https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.rbvi.ucsf.edu%2Fchimera%2Fdocs%2FContributedSoftware%2Fmultalignviewer%2Fframemav.html&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=tWPRdBd0r5mMmFTtdxp7go%2FUuIwKQtYRlNczXMx5kwg%3D&reserved=0">https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.rbvi.ucsf.edu%2Fchimera%2Fdocs%2FContributedSoftware%2Fmultalignviewer%2Fframemav.html&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=tWPRdBd0r5mMmFTtdxp7go%2FUuIwKQtYRlNczXMx5kwg%3D&reserved=0</a>><br>
<<a href="https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.rbvi.ucsf.edu%2Fchimera%2Fdocs%2FContributedSoftware%2Fmultalignviewer%2Fmultalignviewer.html%23superposition&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=8YmSYFDz%2BtbUhDVUEFJ8P1j7HYe1c2qXJE3PTDvjo4k%3D&reserved=0">https://nam02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.rbvi.ucsf.edu%2Fchimera%2Fdocs%2FContributedSoftware%2Fmultalignviewer%2Fmultalignviewer.html%23superposition&data=05%7C01%7Czg2234%40cumc.columbia.edu%7C1f1199f04f894919e31508da6b4d5001%7Cb0002a9b0017404d97dc3d3bab09be81%7C1%7C0%7C637940277051928371%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=8YmSYFDz%2BtbUhDVUEFJ8P1j7HYe1c2qXJE3PTDvjo4k%3D&reserved=0</a>><br>
<br>
I hope this helps,<br>
Elaine<br>
-----<br>
Elaine C. Meng, Ph.D. <br>
UCSF Chimera(X) team<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<br>
> On Jul 21, 2022, at 10:39 AM, Gong, Zhen via ChimeraX-users <chimerax-users@cgl.ucsf.edu> wrote:<br>
> <br>
> Dear Chimerax developers and users, <br>
> <br>
> I am wondering if chimerax is able to do match/matchmaker multiple segments of two protein sequences, something like:<br>
> match #1:3-13, 23-70, 72-81, 86-101 to #2:7-17, 18-66, 67-76, 77-92 cutoffDistance 3 showAlignment true<br>
> Maybe it can but I did not put the command in correct format. Please help.<br>
> <br>
> A related question would be that the default alignment algorithm in chimerax (Needleman-Wunsch or Smith-Waterman) did not do a good job when the two proteins share high structural similarities with low sequence identity. I would like to do the manual sequence
alignment, save the alignment as .aln file and use that for matchmaker, iterate with cutoff and get final rmsd of aligned atoms. Is it possible to execute such procedures in chimerax?<br>
> <br>
> Many thanks in advance for any helpful suggestion and discussion!<br>
> <br>
> Best regards,<br>
> Zhen<o:p></o:p></span></p>
</div>
</div>
</body>
</html>